Perguntas com a marcação «r»

Digamos que eu tenha um experimento fatorial 2x2 simples no qual eu quero fazer a ANOVA. Assim, por exemplo: d <- data.frame(a=factor(sample(c('a1','a2'), 100, rep=T)), b=factor(sample(c('b1','b2'), 100, rep=T))); d$y <- as.numeric(d$a)*rnorm(100, mean=.75, sd=1) + as.numeric(d$b)*rnorm(100, mean=1.2, sd=1) + as.numeric(d$a)*as.numeric(d$b)*rnorm(100, mean=.5, sd=1) + rnorm(100); Na ausência de uma interação significativa, por padrão …

15 r  anova  contrasts 

Estou tentando usar a função ' densidade ' em R para fazer estimativas de densidade do kernel. Estou tendo alguma dificuldade em interpretar os resultados e comparar vários conjuntos de dados, pois parece que a área sob a curva não é necessariamente 1. Para qualquer função de densidade de probabilidade …

15 r  estimation  pdf  kernel-smoothing  auc 

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á 2 anos . Alguém poderia inventar o código R para plotar uma elipse a partir dos valores próprios …

15 r  multivariate-analysis  matrix  matrix-decomposition 

Como você obteria meios horários para várias colunas de dados, por um período diário, e mostraria resultados para doze "hosts" no mesmo gráfico? Ou seja, eu gostaria de representar graficamente como é um período de 24 horas, por uma semana de dados. O objetivo final seria comparar dois conjuntos desses …

15 r  time-series  aggregation 

Atualmente, estou coletando dados para um experimento sobre características psicossociais associadas à experiência da dor. Como parte disso, estou coletando eletronicamente medidas de GSR e BP de meus participantes, juntamente com várias medidas de auto-relato e implícitas. Tenho formação psicológica e me sinto à vontade com análise fatorial, modelos lineares …

15 r  time-series  references 

Estou usando o GNU R em um PC Ubuntu-Lucid que possui 4 CPUs. Para usar todas as 4 CPUs, instalei o pacote "r-cran-multicore". Como o manual do pacote carece de exemplos práticos que eu entendo, preciso de conselhos sobre como otimizar meu script para fazer uso de todas as 4 …

15 r 

Fechadas. Esta questão está fora de tópico . No momento, não está aceitando respostas. Fechado há 4 anos . Bloqueado . Esta pergunta e suas respostas estão bloqueadas porque a questão está fora do tópico, mas tem um significado histórico. No momento, não está aceitando novas respostas ou interações. Existem …

15 r 

Dadas duas matrizes x e y, ambas de comprimento n, ajustei um modelo y = a + b * x e quero calcular um intervalo de confiança de 95% para a inclinação. Este é (b - delta, b + delta) onde b é encontrado da maneira usual e delta = …

15 r  regression  bayesian  confidence-interval  frequentist 

Ao programar em R, usei o pacote multicore algumas vezes. No entanto, nunca vi uma declaração sobre como ele lida com seus números aleatórios. Quando uso o openMP com C, tenho o cuidado de usar um RNG paralelo adequado, mas com o R suponho que algo sensato aconteça. Alguém pode …

15 r  random-generation  parallel-computing  multicore 

Eu tenho um experimento que tentarei abstrair aqui. Imagine que eu jogo três pedras brancas na sua frente e peço que faça um julgamento sobre a posição delas. Registro uma variedade de propriedades das pedras e sua resposta. Eu faço isso em vários assuntos. Eu gero dois modelos. Uma é …

15 r  mixed-model  model-selection 

Estou trabalhando em um projeto que envolve 14 variáveis ​​e 345.000 observações para dados de habitação (itens como ano construído, metragem quadrada, preço vendido, município de residência, etc.). Estou preocupado em tentar encontrar boas técnicas gráficas e bibliotecas R que contenham boas técnicas de plotagem. Já estou vendo o que …

15 r  data-visualization  large-data  eda 

Estou tentando entender a diferença básica entre regressão stepwise e backward em R usando a função step. Para regressão gradual, usei o seguinte comando step(lm(mpg~wt+drat+disp+qsec,data=mtcars),direction="both") Eu tenho a saída abaixo para o código acima. Para seleção de variável para trás, usei o seguinte comando step(lm(mpg~wt+drat+disp+qsec,data=mtcars),direction="backward") E eu tenho a saída …

15 r  regression 

Estou tentando duplicar os resultados da sklearnbiblioteca de regressão logística usando o glmnetpacote em R. A partir da documentação dasklearn regressão logística , ele está tentando minimizar a função de custo sob pena de l2 min w , c 1minw , c12WTw + C∑i = 1Nregistro( exp( - yEu( XTEuw …

15 r  logistic  python  scikit-learn  glmnet 

Uma distribuição Tweedie pode modelar dados assimétricos com uma massa de pontos em zero quando o parâmetro (expoente na relação média-variância) estiver entre 1 e 2.ppp Da mesma forma, um modelo inflado a zero (seja ele contínuo ou discreto) pode ter um grande número de zeros. Estou tendo problemas para …

15 r  generalized-linear-model  prediction  zero-inflation  tweedie-distribution 

Atualmente, estou ajustando florestas aleatórias para um problema de classificação usando o randomForestpacote em R e não tenho certeza sobre como relatar erros de treinamento para esses modelos. Meu erro de treinamento é próximo de 0% quando o computo usando as previsões que recebo com o comando: predict(model, data=X_train) Onde …

15 r  machine-learning  classification  random-forest  overfitting