Uma regra de pontuação adequada é uma regra que é maximizada por um modelo 'verdadeiro' e não permite 'hedging' ou jogo do sistema (relatando deliberadamente resultados diferentes, como é a crença verdadeira do modelo para melhorar a pontuação). A pontuação Brier é adequada, a precisão (proporção...
Digamos que temos um modelo mod <- Y ~ X*Condition + (X*Condition|subject) # Y = logit variable # X = continuous variable # Condition = values A and B, dummy coded; the design is repeated # so all participants go through both Conditions # subject = random effects for different subjects...
Pelo que entendi, as escolas do Reino Unido ensinam que o desvio padrão é encontrado usando: Considerando que as escolas dos EUA ensinam: (em um nível básico de qualquer maneira). Isso causou vários problemas aos meus alunos no passado, enquanto eles procuravam na Internet, mas encontraram a...
Fico impressionado com a idéia do encolhimento de James-Stein (ou seja, que uma função não-linear de uma única observação de um vetor de normais possivelmente independentes pode ser um melhor estimador das médias das variáveis aleatórias, onde 'melhor' é medido por erro ao quadrado ) No entanto,...
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...
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 é que a pedra...
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 confirmar que algo...
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á...
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 = qt(0.975,df=n-2)*se.slope e se.slope é o erro padrão na inclinação....
Eu fiz um modelo de regressão logística usando glm em R. Eu tenho duas variáveis independentes. Como posso plotar o limite de decisão do meu modelo no gráfico de dispersão das duas variáveis. Por exemplo, como plotar uma figura como:
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...
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 CPUs. Meu conjunto de...
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...
Eu brinquei com alguns testes de raiz unitária em R e não sei bem o que fazer com o parâmetro k lag. Usei o teste Dickey Fuller aumentado e o teste Philipps Perron do pacote tseries . Obviamente, o parâmetro padrão (para o ) depende apenas do comprimento da série. Se eu escolher valores k...
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...
Eu estaria interessado em encontrar maneiras em R para atualizar eficientemente um modelo linear quando uma observação ou um preditor é adicionado. O biglm possui um recurso de atualização ao adicionar observações, mas meus dados são pequenos o suficiente para residir na memória (embora eu tenha um...
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) +
Ao fazer um GLM e você obtém o erro "não definido por causa de singularidades" na saída anova, como combater esse erro? Alguns sugeriram que é devido à colinearidade entre covariáveis ou que um dos níveis não está presente no conjunto de dados (consulte: interpretando "não definido por causa de...
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...
Bloqueado . Esta pergunta e suas respostas estão bloqueadas porque a questão está fora de tópico, mas tem um significado histórico. No momento, não está aceitando novas respostas ou interações. Eu tenho o gráfico assim: O código R para gerá-lo é: DF <-...