Dados de contagem com análise de experimentos

Bom dia a todos. Estou tentando analisar o banco de dados anexo com a rotina abaixo: jat1<-read.csv2("...\\jg_09.csv") names(jat1) stopifnot(complete.cases(jat1$ngalhos) !=is.na(jat1$ngalhos))#eliminar os NA's dos dados jat<-jat1[complete.cases(jat1$ngalhos),] Gostaria de verificar o efeito da aração e tplan na contagem de galhos(ngalhos) de uma árvore. Tentei fazer a transformação para "normalidade" a transormação não foi possível. Tentei fazer um teste para verificar se os dados tem distribuição de "poisson" com o pacote vcd, não sei se fiz o correto. Alguém poderia dar algumas sugestões? []'s Edson Lira Estatístico Manaus-Amazonas

A primeira sugestao e' seguir o guia de postagem e nao enviar arquivos anexados. A segunda, vc pode fazer um teste de adequacidade de ajuste (qui-quadrado - que e' o q o goodfit nao-mencionado faz - ou kolmogorov-smirnov). Se vc teve algum problema com ele, apresente-nos os detalhes do que vc tentou, juntamente com o erro... b 2011/5/16 Edson Lira <edinhoestat@yahoo.com.br>:
Bom dia a todos. Estou tentando analisar o banco de dados anexo com a rotina abaixo:
jat1<-read.csv2("...\\jg_09.csv") names(jat1) stopifnot(complete.cases(jat1$ngalhos) !=is.na(jat1$ngalhos))#eliminar os NA's dos dados jat<-jat1[complete.cases(jat1$ngalhos),]
Gostaria de verificar o efeito da aração e tplan na contagem de galhos(ngalhos) de uma árvore.
Tentei fazer a transformação para "normalidade" a transormação não foi possível.
Tentei fazer um teste para verificar se os dados tem distribuição de "poisson" com o pacote vcd, não sei se fiz o correto. Alguém poderia dar algumas sugestões? []'s Edson Lira Estatístico Manaus-Amazonas _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-- Successful people ask better questions, and as a result, they get better answers. (Tony Robbins)

Edson, Ajuste um modelo glm(..., family=poisson) aos dados, faça a análise gráfica de resíduos usual, verifique se a presença de equidispersão, se tudo ok até aqui, peça o quadro de análise de deviance (anova(...)), dependendo do resultado, compare as estimativas dos efeitos. À 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 (3)
-
Benilton Carvalho
-
Edson Lira
-
Walmes Zeviani