Blocos incompletos no R
Pessoal, tomei conhecimento dessa lista através do site do leg. Por acaso alguém saberia me informar se no R existem pacotes ou funções já prontas para analisar os diferentes tipos de Blocos incompletos? Tipo I, Tipo II, Tipo III... Andei conversando com alguns colegas e eles apontaram como saída fazer tudo a mão, montar aquelas tabelas auxiliares e por aí vai. Eu queria saber se não tem como utilizar os comandos normais para os outros delineamentos, criando modelos via LM ou aov. Alguém poderia me ajudar? Grata Julia -- View this message in context: http://r-br.2285057.n4.nabble.com/Blocos-incompletos-no-R-tp3402407p3402407.... Sent from the R-br mailing list archive at Nabble.com.
Júlia, O pacote agricolae tem a função BIB.test() para análise intra de blocos incompletos tipo III. Os demais delineamento são possíveis de analisar usando aov() ou lme() e operando matricialmente com o vetor de estimativas dos efeitos, matriz de delineamento e matriz de covariância das estimativas. Inclusive é possível fazer os testes de médias. Quando precisei analisar dados assim eu fiz uma busca e não encontrei funções. Eu tenho umas rotinas não organizadas para isso. Vou providenciar e colocar à disposição nas Ridículas do R (www.leg.ufpr.br/ridiculas) o mais breve possível. Qual o tipo de experimento que você tem? À disposição. Walmes. Em 24 de março de 2011 09:55, julia_est <juliam_est@yahoo.com.br> escreveu:
Pessoal, tomei conhecimento dessa lista através do site do leg. Por acaso alguém saberia me informar se no R existem pacotes ou funções já prontas para analisar os diferentes tipos de Blocos incompletos? Tipo I, Tipo II, Tipo III...
Andei conversando com alguns colegas e eles apontaram como saída fazer tudo a mão, montar aquelas tabelas auxiliares e por aí vai.
Eu queria saber se não tem como utilizar os comandos normais para os outros delineamentos, criando modelos via LM ou aov.
Alguém poderia me ajudar?
Grata
Julia
-- View this message in context: http://r-br.2285057.n4.nabble.com/Blocos-incompletos-no-R-tp3402407p3402407.... Sent from the R-br mailing list archive at Nabble.com. _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-- ========================================================================== 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: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
Julia, Fico realmente interessado pelas soluções que a lista lhe recomendar. Já adianto que se o Walmes comprou o seu problema alguma coisa interessante vamos poder usar! Volte e meia também me vejo tendo que analisar dados de blocos incompletos (mais especificamente latice, resolvíveis). De qualquer maneira, sempre é bom apreender mais sobre esses delineamentos. No caso específico dos latices resolvíveis em outra oportunidade eu e o Walmes demos uma destrinchada e coneguimos chegar em resultados interessantes do modelo intrabloco e interbloco, o 1º usando o aov e o 2º pelo lme. att, Fernando H On 24/03/2011, Walmes Zeviani <walmeszeviani@gmail.com> wrote:
Júlia,
O pacote agricolae tem a função BIB.test() para análise intra de blocos incompletos tipo III. Os demais delineamento são possíveis de analisar usando aov() ou lme() e operando matricialmente com o vetor de estimativas dos efeitos, matriz de delineamento e matriz de covariância das estimativas. Inclusive é possível fazer os testes de médias. Quando precisei analisar dados assim eu fiz uma busca e não encontrei funções. Eu tenho umas rotinas não organizadas para isso. Vou providenciar e colocar à disposição nas Ridículas do R (www.leg.ufpr.br/ridiculas) o mais breve possível. Qual o tipo de experimento que você tem?
À disposição. Walmes.
Em 24 de março de 2011 09:55, julia_est <juliam_est@yahoo.com.br> escreveu:
Pessoal, tomei conhecimento dessa lista através do site do leg. Por acaso alguém saberia me informar se no R existem pacotes ou funções já prontas para analisar os diferentes tipos de Blocos incompletos? Tipo I, Tipo II, Tipo III...
Andei conversando com alguns colegas e eles apontaram como saída fazer tudo a mão, montar aquelas tabelas auxiliares e por aí vai.
Eu queria saber se não tem como utilizar os comandos normais para os outros delineamentos, criando modelos via LM ou aov.
Alguém poderia me ajudar?
Grata
Julia
-- View this message in context: http://r-br.2285057.n4.nabble.com/Blocos-incompletos-no-R-tp3402407p3402407.... Sent from the R-br mailing list archive at Nabble.com. _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
--
========================================================================== 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: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
Obrigada pelas respostas gente. Eu particularmente tenho os 3 tipos ( I, II e III) de experimento e estou interessadíssima em achar uma sáida pelo R. Walmes gostaria de dar uma olhada nos seus códigos, quando tiver um tempinho coloca lá sim. Gostaria de saber tb como voces fizeram com os Lattices. Eu estou com um conjunto de dados de uma tese do professor Iemma, em que ele utilizou blocos incompletos dentro de parcelas subdivididas. Os tratamentos primários eram dispostos em blocos incompletos. Gostaria muito de tentar uma saída disso pelo R, mas não consigo enxergar um modelo para análise. Estou tendo problemas ao montar via lme() ou aov(). Não sei nem se tem jeito. Obrigada mais uma vez. --- Em qui, 24/3/11, Walmes Zeviani <walmeszeviani@gmail.com> escreveu: De: Walmes Zeviani <walmeszeviani@gmail.com> Assunto: Re: [R-br] Blocos incompletos no R Para: r-br@listas.c3sl.ufpr.br Data: Quinta-feira, 24 de Março de 2011, 10:36 Júlia, O pacote agricolae tem a função BIB.test() para análise intra de blocos incompletos tipo III. Os demais delineamento são possíveis de analisar usando aov() ou lme() e operando matricialmente com o vetor de estimativas dos efeitos, matriz de delineamento e matriz de covariância das estimativas. Inclusive é possível fazer os testes de médias. Quando precisei analisar dados assim eu fiz uma busca e não encontrei funções. Eu tenho umas rotinas não organizadas para isso. Vou providenciar e colocar à disposição nas Ridículas do R (www.leg.ufpr.br/ridiculas) o mais breve possível. Qual o tipo de experimento que você tem? À disposição. Walmes. Em 24 de março de 2011 09:55, julia_est <juliam_est@yahoo.com.br> escreveu: Pessoal, tomei conhecimento dessa lista através do site do leg. Por acaso alguém saberia me informar se no R existem pacotes ou funções já prontas para analisar os diferentes tipos de Blocos incompletos? Tipo I, Tipo II, Tipo III... Andei conversando com alguns colegas e eles apontaram como saída fazer tudo a mão, montar aquelas tabelas auxiliares e por aí vai. Eu queria saber se não tem como utilizar os comandos normais para os outros delineamentos, criando modelos via LM ou aov. Alguém poderia me ajudar? Grata Julia -- View this message in context: http://r-br.2285057.n4.nabble.com/Blocos-incompletos-no-R-tp3402407p3402407.... Sent from the R-br mailing list archive at Nabble.com. _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br -- ========================================================================== 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: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== -----Anexo incorporado----- _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Diga primeiramente como foi o experimento! Valeu!!!! auauauauau Fábio Mathias Corrêa Departamento de Estatística Universidade Estadual de Santa Cruz Tel.: 73-3680-5076 Cel.: 73-9991-8155 ________________________________ De: Julia Maria <juliam_est@yahoo.com.br> Para: r-br@listas.c3sl.ufpr.br Enviadas: Quinta-feira, 24 de Março de 2011 14:48:33 Assunto: Re: [R-br] Blocos incompletos no R Obrigada pelas respostas gente. Eu particularmente tenho os 3 tipos ( I, II e III) de experimento e estou interessadíssima em achar uma sáida pelo R. Walmes gostaria de dar uma olhada nos seus códigos, quando tiver um tempinho coloca lá sim. Gostaria de saber tb como voces fizeram com os Lattices. Eu estou com um conjunto de dados de uma tese do professor Iemma, em que ele utilizou blocos incompletos dentro de parcelas subdivididas. Os tratamentos primários eram dispostos em blocos incompletos. Gostaria muito de tentar uma saída disso pelo R, mas não consigo enxergar um modelo para análise. Estou tendo problemas ao montar via lme() ou aov(). Não sei nem se tem jeito. Obrigada mais uma vez. --- Em qui, 24/3/11, Walmes Zeviani <walmeszeviani@gmail.com> escreveu:
De: Walmes Zeviani <walmeszeviani@gmail.com> Assunto: Re: [R-br] Blocos incompletos no R Para: r-br@listas.c3sl.ufpr.br Data: Quinta-feira, 24 de Março de 2011, 10:36
Júlia,
O pacote agricolae tem a função BIB.test() para análise intra de blocos incompletos tipo III. Os demais delineamento são possíveis de analisar usando aov() ou lme() e operando matricialmente com o vetor de estimativas dos efeitos, matriz de delineamento e matriz de covariância das estimativas. Inclusive é possível fazer os testes de médias. Quando precisei analisar dados assim eu fiz uma busca e não encontrei funções. Eu tenho umas rotinas não organizadas para isso. Vou providenciar e colocar à disposição nas Ridículas do R (www.leg.ufpr.br/ridiculas) o mais breve possível. Qual o tipo de experimento que você tem?
À disposição. Walmes.
Em 24 de março de 2011 09:55, julia_est <juliam_est@yahoo.com.br> escreveu:
Pessoal, tomei conhecimento dessa lista através do site do leg.
Por acaso alguém saberia me informar se no R existem pacotes ou funções já prontas para analisar os diferentes tipos de Blocos incompletos? Tipo I, Tipo II, Tipo III...
Andei conversando com alguns colegas e eles apontaram como saída fazer tudo a mão, montar aquelas tabelas auxiliares e por aí vai.
Eu queria saber se não tem como utilizar os comandos normais para os outros delineamentos, criando modelos via LM ou aov.
Alguém poderia me ajudar?
Grata
Julia
-- View this message in context: http://r-br.2285057.n4.nabble.com/Blocos-incompletos-no-R-tp3402407p3402407.... Sent from the R-br mailing list archive at Nabble.com. _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
--
==========================================================================
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: walmes@ufpr.br twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218
==========================================================================
-----Anexo incorporado-----
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Júlia et. al, Está disponível nas Ridiculas do LEG (http://www.leg.ufpr.br/ridiculas) rotina para análise intra/interbloco de experimento com blocos incompletos tipo III. O código é reproduzível e comentado. Então, peço que avaliem meu código, usem os procedimentos em outros dados e/ou comparem com os resultados de outros programas. Espero o retorno de vocês para aprimorar o conteúdo fornecido. Qualquer erro nas análises é responsabilidade minha e não do R. Assim que possível disponibilizo o tipo I e tipo II. À disposição. Walmes.
Pessoal, não sei se a lsta esta recebendo minhas respostas. TEm algum site que eu possa acompanhar? Pois no meu email tem alguns emails falando que minha resposta aguarda aprovação. Não sei se ja fora ou nao. A mensagem que eu mandei com os dados, retornou que estava muito grande, nao sei se foi ou não. Só hoje me liguei pra isso.
http://www.leg.ufpr.br/doku.php/software:rbr (copie e cole... informacoes gerais sobre a lista) http://r-br.2285057.n4.nabble.com/ (arquivo eletronico da lista q vc pode conultar) b 2011/3/30 Julia Maria <juliam_est@yahoo.com.br>
Pessoal, não sei se a lsta esta recebendo minhas respostas. TEm algum site que eu possa acompanhar? Pois no meu email tem alguns emails falando que minha resposta aguarda aprovação. Não sei se ja fora ou nao.
A mensagem que eu mandei com os dados, retornou que estava muito grande, nao sei se foi ou não. Só hoje me liguei pra isso.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
*conSultar por sinal, ao enviar respostas para a lista, confirme que: 1) vc esta' mandando para r-br@listas.c3sl.ufpr.br 2) vc nao esta' enviando anexos 3) vc nao esta' enviando em HTML b
E complementando, que esteja cadastrada na lista e mandando do email com o qual está cadastrado(a) On Wed, 30 Mar 2011, Benilton Carvalho wrote:
*conSultar por sinal, ao enviar respostas para a lista, confirme que:
1) vc esta' mandando para r-br@listas.c3sl.ufpr.br 2) vc nao esta' enviando anexos 3) vc nao esta' enviando em HTML
b
Obrigado pessoal.. Então uma forma de conferir se a mensagem foi postada é pelo nabble né? Eu estava em duvida, pois o Walmes pediu um .xls pra mim que eu enviei respondendo o email. Mas se ta aparecendo no Nabble deve estar tudo certo.
Pessoal, desculpem se estou insistindo muito, mas agora que comecei a mexer com o R. Eu fiz o procedimento que um usuario me indicou, que me enviou por email porque ele nao conseguiu responder na mensagem aqui na lista. O problema é que eu recebo uma mensagem do tipo:
modelo<-aov(y ~ factor(rep) + factor(A)*factor(B) + factor(bloco)*factor(rep) + Error(factor(bloco)/factor(A)), data=dadosbib2) Mensagens de aviso perdidas: In aov(y ~ factor(rep) + factor(A) * factor(B) + factor(bloco) * : Error() model is singular
Aí eu pergunto se na construção de um modelo, a ordem em que são colocados os fatores alterariam a forma da análise. Eu rodei aqui o procedimento de subdividida que tem na pagina do Eder, e deu a mesma mensagem de Error() model is singular. No meu caso, estou aplicando aos dados que ja enviei no post anterior
Júlia, Essa mensagem de erro ocorre devido a algo relacionado ao subespaço das matrizes envolvidas. Porém não causa problemas na análise (pelo menos no caso balanceado). Uma forma de evitar a mensagem de erro é escrever um novo fator, veja o CMR abaixo da <- expand.grid(A=factor(1:4), B=factor(1:3), rept=factor(1:5)) da$y <- rnorm(da$A) m0 <- aov(y~A+B+Error(A:B), data=da) summary(m0) # anova com mensagem de erro da$parcela <- with(da, factor(paste(A, rept))) # escrevendo um novo fator cujos níveis correspondem à A:rept str(da$parcela) m1 <- aov(y~rept+A*B+Error(parcela), data=da) summary(m1) # mesma anova sem a mensagem de erro À 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: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
participantes (7)
-
Benilton Carvalho -
Fabio Mathias Corrêa -
Fernando Henrique Toledo -
Julia Maria -
julia_est
-
Paulo Justiniano -
Walmes Zeviani