Entendo que o princípio básico por trás do algoritmo para LLE consiste em três etapas. Localizando a vizinhança de cada ponto de dados por alguma métrica como k-nn. Encontre pesos para cada vizinho que denotem o efeito que o vizinho tem no ponto de dados. Construa a incorporação de baixa dimensão...
Estou usando a biblioteca R 'multcomp' ( http://cran.r-project.org/web/packages/multcomp/ ) para calcular o teste de Dunnett. Estou usando o script abaixo: Group <- factor(c("A","A","B","B","B","C","C","C","D","D","D","E","E","F","F","F")) Value <-
Existe um conjugado antes da distribuição de Laplace ? Caso contrário, existe uma expressão de forma fechada conhecida que se aproxima do posterior para os parâmetros da distribuição de Laplace? Eu pesquisei bastante por aí, sem sucesso, então meu palpite atual é "não" nas perguntas acima...
Eu escrevi algum código que pode fazer a filtragem Kalman (usando vários filtros do tipo Kalman [Information Filter et al.]) Para a Análise Linear Gaussiana de Espaço de Estado para um vetor de estado n-dimensional. Os filtros funcionam muito bem e estou obtendo uma saída agradável. No entanto, a...
Eu li sobre o teste t de Student, mas parece funcionar quando podemos assumir que as distribuições originais são normalmente distribuídas. No meu caso, eles definitivamente não são. Além disso, se eu tiver 13 distribuições, preciso fazer
Entendo o conceito de escalar a matriz de dados para usar em um modelo de regressão linear. Por exemplo, em R você pode usar: scaled.data <- scale(data, scale=TRUE) Minha única pergunta é, para novas observações para as quais quero prever os valores de saída, como elas são dimensionadas...
O erro padrão da intercepção termo ( β 0 ) em Y = β 1 x + β 0 + ε é dado por S E ( β 0 ) 2 = σ 2 [ 1β^0β^0\hat{\beta}_0y=β1x+β0+εy=β1x+β0+εy=\beta_1x+\beta_0+\varepsilonSE(β^0)2=σ2[1n+x¯2∑ni=1(xi−x¯)2]SE(β^0)2=σ2[1n+x¯2∑i=1n(xi−x¯)2]SE(\hat{\beta}_0)^2 =...
Estou trabalhando no R através de um excelente tutorial de PCA de Lindsay I Smith e estou ficando paralisado na última etapa. O script R abaixo nos leva ao estágio (na p.19) em que os dados originais estão sendo reconstruídos a partir do (principal neste caso) Componente Principal, que deve gerar...
Se uma variável fator (por exemplo, gênero com níveis M e F) for usada na fórmula glm, as variáveis fictícias serão criadas e poderão ser encontradas no resumo do modelo glm juntamente com seus coeficientes associados (por exemplo, gêneroM) Se, em vez de confiar em R para dividir o fator dessa...
Fiz essa pergunta com (n,k)=(400,220)(n,k)=(400,220)(n, k) = (400, 220) em uma entrevista. Existe uma resposta "correta"? Suponha que os lançamentos sejam iid e a probabilidade de cabeças é p=0.5p=0.5p=0.5 . A distribuição do número de cabeças em 400 lançamentos deve ser próxima de Normal (200, 10...
Estive analisando alguns dados usando a modelagem linear de efeitos mistos em R. Estou planejando fazer um pôster com os resultados e fiquei imaginando se alguém com modelos de efeitos mistos poderia sugerir quais gráficos usar para ilustrar os resultados do modelo. Eu estava pensando em gráficos...
Eu tenho um conjunto de dados supondo que os vizinhos mais próximos são os melhores preditores. Apenas um exemplo perfeito de gradiente bidirecional visualizado- Suponha que tenhamos um caso em que faltam poucos valores, podemos prever facilmente com base em vizinhos e tendências. Matriz de...
Eu calculei a autocorrelação em dados de séries temporais sobre os padrões de movimento de um peixe com base em suas posições: X ( x.ts) e Y ( y.ts). Usando R, executei as seguintes funções e produzi os seguintes gráficos: acf(x.ts,100) acf(y.ts,100) Minha pergunta é: como interpreto...
Em geral, estou me perguntando se é sempre melhor não usar polinômios ortogonais ao ajustar uma regressão com variáveis de ordem superior. Em particular, eu estou pensando com o uso de R: Se poly()com raw = FALSEproduz os mesmos valores ajustados como poly()com raw = TRUE, e polycom raw =...
fundo Estou tentando entender o primeiro exemplo de um curso sobre montagem de modelos (então isso pode parecer ridiculamente simples). Fiz os cálculos manualmente e eles correspondem ao exemplo, mas quando os repito em R, os coeficientes do modelo estão desativados. Eu pensei que a diferença pode...
Eu tenho uma distribuição empírica . Eu calculo da seguinte formaG(x)G(x)G(x) x <- seq(0, 1000, 0.1) g <- ecdf(var1) G <- g(x) Denomino , ou seja, h é o pdf enquanto G é o cdf.h(x)=dG/dxh(x)=dG/dxh(x) = dG/dxhhhGGG Agora, quero resolver uma equação para o limite superior de...
Eu tenho dados sobre quantos usuários publicam quantas perguntas. Por exemplo, [UserCount, QuestionCount] [2, 100] [9, 10] [3, 80] ... ... Isso significa que 2 usuários postaram 100 perguntas, 9 usuários postaram 10 perguntas e assim por diante. Então, como posso determinar se a UserCount,...
Eu obtive um modelo de regressão logística (via train) para uma resposta binária e a matriz de confusão logística via confusionMatrixin caret. Isso me fornece a matriz de confusão do modelo logístico, embora não tenha certeza de qual limite está sendo usado para obtê-lo. Como obtenho a matriz de...
Eu tenho lido sobre o cálculo de valores de R2R2R^2 em modelos mistos e depois de ler as perguntas frequentes do R-sig, outras postagens neste fórum (eu vincularia algumas, mas não tenho reputação suficiente) e várias outras referências que eu entendo que usando valores no contexto de modelos...
Estou lendo um comentário em um artigo, e o autor afirma que, às vezes, mesmo que os estimadores (encontrados por ML ou quase-probabilidade máxima) possam não ser consistentes, o poder de um teste de razão de verossimilhança ou razão de quase-verossimilhança ainda pode convergir para 1 como o...