[R-br] Usando o termo Error() com model.matrix para simular dados de Medidas Repetidas
Ola, normalmente para tentar entender melhor algumas analises, eu tento gerar exemplo e ver como ficam os resultados e gráficos. Normalmente eu uso o seguinte método: ### #gerando os tratamentos example-data.frame(time=factor(rep(c(0,1,2),each=60)), trat=rep(rep(c(controle,dose2,dose3),each=20),each=60), id=factor(rep(1:60,3))) str(example) #Usando model.matrix para fazer a matriz que multiplica os parâmetros, se quiser simular uma two way matrix.1-model.matrix(~time*trat,data=example) head(matrix.1) #Gerando valores dos parâmetros param.1-c(1:9) #Gerando os dados de resposta, que é a multiplicação da matriz pelos valores dos parâmetros response.1-rnorm(nrow(matrix.1),matrix.1%*%param.1,1) #Assim a analise vai estimar os parâmetros usados, e eu posso ver vários parâmetros diferentes e o que acontece. summary(lm(response.1~example$time*example$trat)) #Mas eu gostaria de simular algo com medidas repetidas, tipo usando o termo Error() que da para usar com a função AOV. #Mas eu fiquei meio perdido como proceder, para fazer o mesmo esquema acima mas com medidas repetidas, tendo o error term. matrix.2-model.matrix(~time*trat+Error(id),data=example) #Olhando como se faz a mesma coisa com a função gls do pacote nlme, suspeito que teria que fazer uma matriz de correlação, mas ai eu não sei como faria aquela multiplicação la em cima. #Alguém saberia como proceder? # Os comando para fazer a medidas repetidas com nlme que me refiro são desse tutorial: http://www.ats.ucla.edu/stat/r/seminars/Repeated_Measures/repeated_measures.htm Olhando isso e a função aov eu vejo que tenho que separar o Error, mas não estou sabendo como. library(nlme) longg - groupedData(pulse ~ as.numeric(exertype) * as.numeric(time) | id, data = exer) fit.cs - gls(pulse ~ exertype * time, data = longg, corr = corCompSymm(, form= ~ 1 | id) ) summary(fit.cs) Agradeço a atenção Abraços Augusto Ribas -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.wordpress.com/ Lattes: http://lattes.cnpq.br/7355685961127056 ___ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
Re: [R-br] Gerar mapa de distribuição por UF
Então, Elias, Eu montei a matriz da seguinte forma: - Uma coluna com os municípios; - Em seguidam uma coluna com o Nr de casos/Município - Depois, umas 6 colunas com variáveis ambientais. E fui repetindo desta forma para cada ano (empilhando) de 2002 a 2011... Minha ideia é montar um modelo que tenha Nr de casos/mun = f( variáveis ambientais) Depois, quero jogar isto num shapefile usando o TerraView para mostra a evolução dos casos naqueles anos, por município. É possível??? Abs. De: Elias Teixeira Krainski eliaskrain...@yahoo.com.br Para: r-br@listas.c3sl.ufpr.br Enviadas: Domingo, 3 de Março de 2013 17:35 Assunto: Re: [R-br] Gerar mapa de distribuição por UF Como eu entro com os dados no R para obter uma modelagem de dados espaço-temporais Tenho um arquivo onde nas colunas estão as variáveis que quero trabalhar e nas linhas os municípios. Só que eu repito estes municípios para os anos de 2002 a 2011. É assim mesmo Que tipo de modelo vc tem em mente? Eu tenho usado INLA, não sei como trabalham outros pacotes. No INLA trabalhamos com indices. No seu caso, é só 'empilhar' os dados. Assim vc pode ter um vetor de dados (em vez de matriz), um índice para área e outro para tempo. Na classe dos modelos dinâmicos, veja exemplos em modelos dinâmicos espaço temporais em http://www.sciencedirect.com/science/article/pii/S0167947311003999 e maiores detalhes da implementação em http://www.r-inla.org/examples/case-studies/ruiz-cardenas-kraonski-and-rue-2010 O INLA esta em constante desenvolvimento e atualmente é possível simplificar um nível nos códigos nesses links. Por exemplo, usando a extratégia de definir matriz de precisão pelo produto de Kronecker, por exemplo, entre um efeito espacial e um temporal. Na classe de modelos de interação espaço temporais do Held vc pode usar a estratégia de Kronecker. Para interação do tipo 4, você pode fazer simplesmente f(i, model='besag', graph=grafo, group=t, control.group=list(model='rw1')) em que 'i' indica as áreas e 't' o tempo. ___ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.___ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
Re: [R-br] Gerar mapa de distribuição por UF
Oi Gilbert, Essa dica que eu dei foi para a estimação de modelos espaço temporais com o INLA. Não sei se será útil para visualização no TerraView. Pois, trabalhei muito pouco com o TerraView e já faz alguns anos... Para visualização, vc pode usar o pacote spacetime. Neste seu caso, veja a função STFDF(). Abs, Elias. ___ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
Re: [R-br] Gerar mapa de distribuição por UF
Eu não conhece direito o INLA. Baixei alguma coisa da internet sobre ele, mas relamente reconheço que foi muito pouco. O que eu tenho visto bastante é o uso do WinBugs e tb do BayesX para modelagem destas dados. Mas estou pensando em apostar minhas fichas no INLA, pois parece ser mais fácil do que outros dois, em virtude da própria modelagem p/ uso destes pacotes. Vc tem mais alguma coisa sobre o INLA??? Abs. De: Elias Teixeira Krainski eliaskrain...@yahoo.com.br Para: r-br@listas.c3sl.ufpr.br Enviadas: Segunda-feira, 4 de Março de 2013 10:23 Assunto: Re: [R-br] Gerar mapa de distribuição por UF Oi Gilbert, Essa dica que eu dei foi para a estimação de modelos espaço temporais com o INLA. Não sei se será útil para visualização no TerraView. Pois, trabalhei muito pouco com o TerraView e já faz alguns anos... Para visualização, vc pode usar o pacote spacetime. Neste seu caso, veja a função STFDF(). Abs, Elias. ___ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.___ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
[R-br] teste de Bartlett
Boa tarde a todos, gostaria de saber se alguem sabe se para teste de Bartlett quando faz para fatorial usa os dois fatores juntos ou faz para cada um dos fatores separados? Obrigado a todos deste já. -- *José Henrique Soler Guilhen* Graduando em Agronomia Universidade Federal do Espírito Santo - CCAUFES ___ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
Re: [R-br] teste de Bartlett
Você pode considerar as combinações entre os níveis como um único fator ou fazer para cada fator em separado desde que faça com os resíduos do modelo ajustado. À disposição. Walmes. == Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: wal...@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 == ___ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
Re: [R-br] teste de Bartlett
Lembrando que todos os testes implementados em softwares que existem são one-way. Em 2007 se não me engano propuseram uma extensão do teste de Levene para análises multi-way. (S,f,P) Allaman ___ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
Re: [R-br] teste de Bartlett
Obrigado pela ajuda Walmes e todos os outros. Em 4 de março de 2013 14:42, Walmes Zeviani walmeszevi...@gmail.comescreveu: Você pode considerar as combinações entre os níveis como um único fator ou fazer para cada fator em separado desde que faça com os resíduos do modelo ajustado. À disposição. Walmes. == Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: wal...@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 == ___ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. -- *José Henrique Soler Guilhen* Graduando em Agronomia Universidade Federal do Espírito Santo - CCAUFES ___ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
