Perguntas com a marcação «r»

O bildpacote parece ser um excelente pacote para respostas binárias em série. Mas é por tempo discreto. Gostaria de especificar uma função suave do tempo para a conexão do odds ratio da resposta atual Y com respostas binárias medidas em tempos anteriores, ou pelo menos uma versão Markov de primeira …

13 r  repeated-measures  binary-data  panel-data 

Estou procurando um pacote R geral, limpo e rápido (ou seja, usando rotinas C ++) para simular caminhos de uma difusão não-linear não homogênea como (1) usando o esquema de Euler-Maruyama, o esquema de Milstein (ou qualquer outro). Isso está destinado a ser incorporado a um código de estimativa maior …

13 r  simulation  stochastic-processes  markov-process 

Quero prever o Tree Heights em uma determinada área usando algumas variáveis ​​obtidas por sensoriamento remoto. Como a biomassa aproximada, etc. Quero primeiro usar uma regressão linear (sei que não é a melhor ideia, mas é um passo obrigatório para o meu projeto). Eu queria saber o quanto a autocorrelação …

13 r  multiple-regression  spatial  autocorrelation 

O "Aprendizado de Máquina para Hackers" de O'Reilly diz que cada componente principal representa uma porcentagem da variação. Citei a parte relevante da página abaixo (capítulo 8, p.207). Falando com outro especialista, eles concordaram que é a porcentagem. No entanto, os 24 componentes somam 133,2095%. Como pode ser? Tendo nos …

13 r  pca 

Essas duas funções existem no R, mas não conheço suas diferenças. Parece que eles retornam apenas os mesmos valores p ao chamar wilcox.testcom correct=FALSE, e wilcox_test(no pacote de moedas) com distribution="aymptotic". Para outros valores, eles retornam valores p diferentes. Também wilcox.testestá sempre retornando W = 0 para meu conjunto de …

13 r  hypothesis-testing  p-value  wilcoxon-mann-whitney 

Estou tentando executar um logit bayesiano nos dados aqui . Estou usando bayesglm()no armpacote em R. A codificação é direta o suficiente: df = read.csv("http://dl.dropbox.com/u/1791181/bayesglm.csv", header=T) library(arm) model = bayesglm(PASS ~ SEX + HIGH, family=binomial(link="logit"), data=df) summary(model) fornece a seguinte saída: Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) 0.10381 …

13 r  bayesian  p-value 

Eu adoraria realizar um teste post-hoc de TukeyHSD após meu Anova de duas vias com R, obtendo uma tabela contendo os pares classificados agrupados por diferença significativa. (Desculpe pela redação, ainda sou novo em estatísticas.) Eu gostaria de ter algo parecido com isto: Então, agrupados com estrelas ou letras. Qualquer …

13 r  anova  multiple-comparisons  post-hoc  tukey-hsd 

Estou tentando modelar alguns dados usando o glmnetpacote em R. Digamos que tenho os seguintes dados training_x <- data.frame(variable1 = c(1, 2, 3, 2, 3), variable2 = c(1, 2, 3, 4, 5)) y <- c(1, 2, 3, 4, 5) (Isso é uma simplificação; meus dados são muito mais complicados.) Então …

13 r  glmnet 

Estou tentando fazer previsões usando um modelo aleatório de floresta em R. No entanto, recebo erros, pois alguns fatores têm valores diferentes no conjunto de testes e no conjunto de treinamento. Por exemplo, um fator Cat_2possui valores 34, 68, 76, etc., no conjunto de testes que não aparecem no conjunto …

13 r  categorical-data  random-forest 

Eu queria saber como você geraria dados de uma equação de regressão de Poisson em R? Estou meio confuso sobre como abordar o problema. Portanto, se eu assumir que temos dois preditores e X 2 que são distribuídos . E a interceptação é 0 e ambos os coeficientes são iguais …

13 r  regression  poisson-distribution  simulation 

Em resposta a uma pergunta sobre seleção de modelos na presença de multicolinearidade , Frank Harrell sugeriu : Coloque todas as variáveis ​​no modelo, mas não teste o efeito de uma variável ajustada para os efeitos das variáveis ​​concorrentes ... Os testes de chunk de variáveis ​​concorrentes são poderosos porque …

13 r  model-selection  multicollinearity 

Essa também é uma questão que se relaciona muito com o Python como bancada de estatística e o Excel como bancada de estatística . Eu sei que há uma enorme discussão sobre Ruby versus Python, mas este não é o ponto nesta questão. Eu pensei que Ruby ser mais rápido …

13 r  python  software  ruby 

Estou jogando com um randomForest e descobri que geralmente aumentar o sampSize leva a um melhor desempenho. Existe uma regra / fórmula / etc que sugere qual deve ser o sampSize ideal ou é uma tentativa e erro? Eu acho que outra maneira de expressar isso; Quais são os meus …

13 r  random-forest 

Eu tenho um conjunto de dados contendo várias proporções que somam 1. Estou interessado na alteração dessas proporções ao longo de um gradiente (veja abaixo, por exemplo, dados). gradient <- 1:99 A1 <- gradient * 0.005 A2 <- gradient * 0.004 A3 <- 1 - (A1 + A2) df <- …

13 r  multinomial 

Alguém pode me explicar meu modelo de Cox em inglês simples? Eu ajustei o seguinte modelo de regressão de Cox em todos os meus dados usando a cphfunção Meus dados são salvos em um objeto chamado Data. As variáveis w, xe ysão contínuos; zé um fator de dois níveis. O …

13 r  cross-validation  survival  cox-model