Eu tenho uma distribuição de amostras com um pequeno número de valores em cada um (menos de ). Eu calculei a mediana para cada amostra, que quero comparar com um modelo e obter a diferença entre o modelo e a mediana de cada amostra. Para ter um resultado consistente, preciso de um erro nessa...
Estou com problemas para obter resultados equivalentes entre um aovmodelo de medidas repetidas entre dentro e um lmermodelo misto. Meus dados e script têm a seguinte aparência data=read.csv("https://www.dropbox.com/s/zgle45tpyv5t781/fitness.csv?dl=1") data$id=factor(data$id) data id FITNESS TEST...
Quero calcular / avaliar a convolução g(x)=∫Df(x−t)ϕ(t)dt,g(x)=∫Df(x−t)ϕ(t)dt,g(x)=\int_D f(x-t) \phi(t) dt, onde é um densidade e é uma função suave com suporte compacto . A convolução não está disponível em formato fechado e preciso integrá-la numericamente. Minha pergunta é: Existe uma maneira...
Eu usei a função 'polr' no pacote MASS para executar uma regressão logística ordinal para uma variável de resposta categórica ordinal com 15 variáveis explicativas contínuas. Utilizei o código (mostrado abaixo) para verificar se meu modelo atende à premissa proporcional de chances, seguindo as...
A versão curta fora do contexto Seja uma variável aleatória com CDF yyyF(⋅)≡{θθ+(1−θ)×CDFlog-normal(⋅;μ,σ) y = 0 y > 0F(⋅)≡{θ y = 0 θ+(1−θ)×CDFlog-normal(⋅;μ,σ) y > 0 F(\cdot) \equiv \cases{\theta & y = 0 \\ \theta + (1-\theta) \times \text{CDF}_{\text{log-normal}}(\cdot; \mu, \sigma) & y >...
Em um modelo misto linear generalizado logístico (família = binomial), não sei como interpretar a variação de efeitos aleatórios: Random effects: Groups Name Variance Std.Dev. HOSPITAL (Intercept) 0.4295 0.6554 Number of obs: 2275, groups: HOSPITAL, 14 Como interpreto esse resultado...
Eu tenho três variáveis macroeconômicas (ICS - sentimento do consumidor, ER - taxa de emprego, DGO - pedido de bens duráveis) e executei testes de causalidade de Granger em R sobre eles. Eu realmente não sei como interpretar os resultados de um teste de Granger. Alguém poderia me ajudar com a...
Estou tentando estimar o efeito médio do tratamento a partir de dados observacionais usando a ponderação do escore de propensão (especificamente IPTW). Acho que estou calculando o ATE corretamente, mas não sei como calcular o intervalo de confiança do ATE, considerando os pesos da pontuação de...
Tenho dados do seguinte desenho experimental: minhas observações são contagens do número de sucessos ( K) do número correspondente de tentativas ( N), medidos para dois grupos, cada um composto por Iindivíduos, a partir de Ttratamentos, nos quais em cada combinação de fatores existem Rrepetições ....
Eu tenho o seguinte conjunto de dados simples com duas variáveis contínuas; ou seja: d = data.frame(x=runif(100,0,100),y = runif(100,0,100)) plot(d$x,d$y) abline(lm(y~x,d), col="red") cor(d$x,d$y) # = 0.2135273 Preciso reorganizar os dados para que a correlação entre as variáveis seja ~...
Eu estava tentando reaprender as estatísticas bayesianas (toda vez que achava que finalmente as recebia, surge outra coisa que não havia considerado antes ...), mas não estava claro (para mim) o que o processo de geração de dados na estrutura bayesiana, na verdade é. A estrutura frequentista é...
Ao selecionar um número apropriado de nós para um GAM, convém levar em consideração o número de dados e incrementos no eixo x. E se tivermos 100 incrementos no eixo x com 1000 pontos de dados em cada incremento. A informação aqui diz: Se eles não forem fornecidos, os nós do spline serão...
Ao estimar parâmetros com MAP, por que está escrito que estamos estimando o "modo"? Eu pensei que seria o meio da distribuição
Estou procurando um pacote para me ajudar a resolver alguns problemas de otimização quadrática e vejo que há pelo menos meia dúzia de pacotes diferentes. De acordo com esta página: QP (programação quadrática, 90C20): cplexAPI , kernlab , limSolve , LowRankQP , quadprog , Rcplex , Rmosek Alguns...
Meus dados brutos consistem em uma série temporal de 60 dias com uma tendência de queda. Os dados são semanais, portanto a frequência é definida como 7. Eu calculei a diferença dos dados que se parece com isso Quando executo gráficos de ACF e PACF sobre a diferença, pareço obter resultados...
Uso a decomposefunção Re crio os três componentes da minha série temporal mensal (tendência, sazonal e aleatória). Se plotar o gráfico ou olhar para a tabela, posso ver claramente que a série temporal é afetada pela sazonalidade. No entanto, quando regredo a série temporal para as 11 variáveis...
Suponha que, em uma análise de regressão em R, eu tenha uma variável independente do tipo de fator com 3 níveis no meu conjunto de dados de trem. Mas no conjunto de dados de teste essa mesma variável de fator tem 5 níveis. Portanto, não posso prever os valores de resposta para o conjunto de dados...
Estou tentando entender a saída da análise de componentes principais executada da seguinte maneira: > head(iris) Sepal.Length Sepal.Width Petal.Length Petal.Width Species 1 5.1 3.5 1.4 0.2 setosa 2 4.9 3.0 1.4 0.2 setosa 3 4.7 3.2 1.3 0.2 setosa 4 4.6 3.1 1.5 0.2 setosa 5 5.0 3.6 1.4 0.2...
É sabido que uma combinação linear de 2 variáveis normais aleatórias também é uma variável normal aleatória. Existem famílias de distribuição não normais comuns (por exemplo, Weibull) que também compartilham essa propriedade? Parece haver muitos contra-exemplos. Por exemplo, uma combinação linear...
Eu tenho um conjunto de dados com variável contínua e uma variável de destino binário (0 e 1). Preciso discretizar as variáveis contínuas (para regressão logística) com relação à variável alvo e com a restrição de que a frequência de observação em cada intervalo seja equilibrada. Tentei...