Perguntas com a marcação «r»

Eu estou querendo saber se existem algoritmos (talvez algoritmos genéticos) em R para construção de recursos (derivados de preditores candidatos de preditores existentes)? Estou pensando em uma rotina para testar potências, interações, proporções e combinações lineares de ordem superior e funções não lineares de variáveis ​​existentes (sin, cos, atan etc.). …

8 r  machine-learning  feature-construction  genetic-algorithms 

Acabei de descobrir como executar um script R no R Console no Windows. source("arrrFile.R") O problema é que esse comando é executado "silenciosamente". Como posso executar o arquivo comando por comando, como eu digitaria no console?

8 r 

Eu estava brincando com o ggplot2 usando os seguintes comandos para ajustar uma linha aos meus dados: ggplot(data=datNorm, aes(x=Num, y=Val)) + geom_point() + stat_summary(fun.data = "mean_cl_boot", geom="errorbar", colour="red", width=0.8) + stat_sum_single(median) + stat_sum_single(mean, colour="blue") + geom_smooth(level = 0.95, aes(group=1), method="lm") Os pontos vermelhos são valores medianos, o azul é a …

8 r  modeling  dataset  data-mining  ggplot2 

Eu tenho um experimento de curso de tempo que segue 8 grupos de tratamento de 12 peixes por 24 horas com observações feitas em intervalos de 5 segundos. Entre as medições feitas está o quão longe cada peixe viaja (em mm) entre as observações. As 24 horas são divididas em …

8 r  time-series  pattern-recognition 

Vou dar meus exemplos com chamadas R. Primeiro, um exemplo simples de regressão linear com uma variável dependente 'vida útil' e duas variáveis ​​explicativas contínuas. data.frame(height=runif(4000,160,200))->human.life human.life$weight=runif(4000,50,120) human.life$lifespan=sample(45:90,4000,replace=TRUE) summary(lm(lifespan~1+height+weight,data=human.life)) Call: lm(formula = lifespan ~ 1 + height + weight, data = human.life) Residuals: Min 1Q Median 3Q Max -23.0257 -11.9124 …

8 r  multiple-regression  interpretation  regression-coefficients 

Fechadas. Esta questão está fora de tópico . No momento, não está aceitando respostas. Deseja melhorar esta pergunta? Atualize a pergunta para que ela esteja no tópico de Validação cruzada. Fechado há 8 meses . Teste A / B: http://20bits.com/articles/statistical-analysis-and-ab-testing/ http://elem.com/~btilly/effective-ab-testing/ Eu não estou muito familiarizado com o teste A …

8 r  python  ab-test 

Eu gostaria de fazer algo em R que o SAS possa fazer usando o proc misto do SAS (existe alguma maneira de se fazer isso no STATA), ou seja, ajustar o chamado modelo bivariado de Reitsma et al (2005). Este modelo é um modelo misto especial em que a variação …

8 r  optimization 

Eu tenho alguns dados que estão claramente truncados à esquerda. Desejo ajustá-lo com uma estimativa de densidade que o manipule de alguma maneira, em vez de tentar suavizá-lo. Quais métodos conhecidos (como de costume, em R) podem resolver isso? Código de amostra: set.seed(1341) x <- c(runif(30, 0, 0.01), rnorm(100,3)) hist(x, …

8 r  pdf  histogram  kernel-smoothing 

Eu tenho um aplicativo para o qual eu preciso de uma aproximação à soma lognormal pdf para usar como parte de uma função de probabilidade. A distribuição da soma lognormal não tem forma fechada, e há vários papéis em periódicos de processamento de sinais sobre diferentes aproximações. Eu tenho usado …

8 r  lognormal 

Estou interessado em montar uma ANOVA Bayesiana de Dois Fatores em ERROS ou utilizando algum pacote R. Infelizmente, estou tendo dificuldade em encontrar recursos sobre este tópico. Alguma sugestão? Mesmo um artigo descrevendo a abordagem seria útil.

8 r  bayesian  anova  bugs 

Eu calculei a proporção de pintos nascidos do número de ovos nascidos em cada ano usando prop.test()em R. Eu vejo que isso me dá a proporção nascida, mas também o intervalo de confiança de 95%, que é o que eu estou procurando. Depois de ler as excelentes informações de outra …

8 r  confidence-interval  binomial 

Pesquisei este online por horas, mas nenhuma postagem online é o que estou procurando. Minha pergunta é muito fácil de implementar no procedimento misto SAS Proc, mas não sei como fazê-lo nos pacotes lme e / ou lmer. Suponha, eu tenho um modelo, , em que é fixo, mas ey=μ+α+β+αβ+ey=μ+α+β+αβ+ey …

8 r  mixed-model  cross-validation  biostatistics 

Suponha que uma variável de nível dependente seja medida em um nível de unidade (nível 1) aninhado em unidades do tipo (nível 2 ) e que unidades do tipo A estejam aninhadas em níveis do tipo B (nível 3 ).yyyAAA222AAABBB333 Suponha que eu ajuste a seguinte fórmula: y ~ "FIXED …

8 r  lme4-nlme  random-effects-model 

Costumo enfrentar a questão de ter que escolher um número de clusters. A partição que acabo escolhendo é mais frequentemente baseada em preocupações visuais e teóricas do que em critérios de qualidade. Eu tenho duas perguntas principais. O primeiro diz respeito à idéia geral de qualidade dos clusters. Pelo que …

8 r  machine-learning  classification  clustering  hierarchical-clustering 

Eu já vi ótimos posts explicando o PCA e por que, sob essa abordagem, os autovetores de uma matriz de correlação (simétrica) são ortogonais. Também entendo as maneiras de mostrar que esses vetores são ortogonais entre si (por exemplo, pegar os produtos cruzados da matriz desses vetores próprios resultará em …

8 r  correlation  pca  orthogonal