Perguntas com a marcação «r»

Tenho algumas perguntas básicas sobre PCA (análise de componentes principais) e LDA (análise discriminante linear): No PCA, existe uma maneira de calcular a proporção de variação explicada. Também é possível para o LDA? Se sim, como? A saída “Proporção de rastreio” da ldafunção (na biblioteca R MASS) é equivalente à …

8 r  variance  pca  discriminant-analysis 

Meu design é o seguinte. yyy é a resposta de Bernoulli x1x1x_1 é uma variável contínua x2x2x_2 é uma variável categórica (fator) com dois níveis O experimento é completamente dentro dos sujeitos. Ou seja, cada sujeito recebe cada combinação de e .x 2x1x1x_1x2x2x_2 Esta é uma configuração de regressão logística …

8 r  repeated-measures  glmm  lme4-nlme 

Eu tenho um conjunto de dados meteorológicos diários, que surpreendentemente tem um efeito sazonal muito forte. Adaptei um modelo ARIMA a esse conjunto de dados usando a função auto.arima do pacote de previsão. Para minha surpresa, a função não aplica operações sazonais - diferencial sazonal, componentes sazonais ar ou ma. …

8 r  time-series  forecasting  arima  arma 

Fiz um agrupamento de pontos de coordenadas (longitude, latitude) e encontrei resultados adversos surpreendentes dos critérios de agrupamento para o número ideal de agrupamentos. Os critérios são retirados do clusterCrit()pacote. Os pontos que estou tentando agrupar em um gráfico (as características geográficas do conjunto de dados são claramente visíveis): O …

8 r  clustering 

Alguém conhece um método para comparar duas variáveis ​​com um teste do qui quadrado se as variáveis ​​são de pesquisas diferentes com svydesign()declarações diferentes ? Estou procurando testar uma diferença em uma distribuição variável em duas ondas de uma pesquisa, mas a svychisq()declaração é limitada a um objeto de design. …

8 r  chi-squared  survey 

Quero incluir pesos de amostra no meu modelo de regressão quantil, mas não sei como fazer isso. Já defini meu peso, que são pesos replicados já fornecidos no conjunto de dados da pesquisa (calculado no pacote de pesquisa): w<-svrepdesign(variables=data[,1:10],repweights=data[,11:30],type="BRR", combined.weights=TRUE, weights=r.weights, rho=0.5,dbname="") e meu modelo rq é: rq(y~x,tau=c(.1,.2,.3,.4,.5,.6,.7,.8,.9),data=my.data)) Eu tentei …

8 r  quantile-regression 

Eu preciso de um único conjunto de dados imputados (por exemplo, para criar um manequim de grupo de países a partir dos dados de renda per capita do país imputado). O R oferece pacotes de pacotes para criar vários dados imputados (por exemplo, Amelia) e combinar resultados de vários conjuntos …

8 r  data-imputation 

Estou tentando executar uma regressão de Cox em um exemplo de conjunto de dados de 2.000.000 de linhas da seguinte maneira, usando apenas R. Esta é uma tradução direta de um PHREG no SAS. A amostra é representativa da estrutura do conjunto de dados original. ## library(survival) ### Replace 100000 …

8 r  logistic  survival  sas  cox-model 

Eu tenho 2 estudos analisando a resposta do paciente ao mesmo medicamento. O estudo 1 encontrou 10.000 genes expressos acima do plano de fundo e 500 deles são diferencialmente expressos e referidos como a assinatura da resposta ao medicamento. O estudo 2 encontrou 1.000 genes representando a assinatura da resposta …

8 r  multiple-comparisons  genetics  bonferroni  hypergeometric 

Por favor, permita-me fazer uma pergunta básica. Entendo a mecânica de Naive Bayes para variáveis ​​discretas e posso refazer os cálculos "à mão". (código de HouseVotes84 todo o caminho abaixo). No entanto - estou lutando para ver como a mecânica funciona para variáveis ​​contínuas (código de exemplo abaixo). Como o …

8 r  naive-bayes  bayes 

Eu gostaria de fazer análises de poder para testes de hipóteses de (não) igualdade de proporções nas quais as proporções são muito pequenas. Eu gostaria de fazer isso sem usar aproximações normais (ou Poisson) da distribuição binomial. Existem vários tipos gerais de questões de poder que eu gostaria de resolver. …

8 r  binomial  power-analysis 

Quando executo este código: require(nlme) a <- matrix(c(1,3,5,7,4,5,6,4,7,8,9)) b <- matrix(c(3,5,6,2,4,6,7,8,7,8,9)) res <- lm(a ~ b) print(summary(res)) res_gls <- gls(a ~ b) print(summary(res_gls)) Eu obtenho os mesmos coeficientes e a mesma significância estatística nos coeficientes: Loading required package: nlme Call: lm(formula = a ~ b) Residuals: Min 1Q Median 3Q …

8 r  least-squares  generalized-least-squares 

Eu tenho alguns conjuntos de dados EEG que estou testando em duas classes. Posso obter uma taxa de erro decente do LDA (as distribuições condicionais de classe não são gaussianas, mas têm caudas semelhantes e separação suficientemente boa) e, portanto, quero plotar o ROC do preditor LDA contra conjuntos de …

8 r  roc  discriminant-analysis 

Preciso calcular o CDF Dirichlet , mas só consigo encontrar implementações do PDF . Vocês conhecem alguma biblioteca (de preferência a R) implementando-a?

8 r  matlab  pdf  cdf  dirichlet-distribution 

Esta é a pergunta de um iniciante em um exercício na “Computação Bayesiana com R” de Jim Albert. Observe que, embora isso possa ser um dever de casa, no meu caso não é, pois estou aprendendo métodos bayesianos em R porque acho que posso usá-lo em minhas análises futuras. De …

8 r  bayesian  simulation