Prezados,
Fiz minha "inscrição" a pouco e não sei se fiz correto, mas vamos lá.
Sou iniciante no R, e estou fazendo umas análises de regressões nonlinear , LL, testando doses de um herbicida no comprimento da raiz do arroz (cm).
Como fiz dois experimentos, estou vendo se ambos são iguais para discutir em conjunto. Já vi que não as curvas fitted não são iguais, tanto utilizando os dados originais (cm) ou usando dados transf. em %. Ou seja, o p <0.05.
Tentei fazer uma transformação com box-cox, mas não fui"eficiente" em usar o R. Poderia me dar uma dica como usar esta ferramenta?
Tenho no file uma cultivar, 7 doses, root , divididos em dois estudos O que fiz inicialmente foi testar a logistica LL.3:
cult<drm(root~dose, data=cultivar, fct=LL.3())
mas com os dados originais ou em % não consigo P>0.05. POr isto gostaria de testar box-cox.
obrigado..Andre
Em quarta-feira, 10 de abril de 2013 12:35:38, Andre Andres escreveu:
Prezados,
Fiz minha "inscrição" a pouco e não sei se fiz correto, mas vamos lá.
Sou iniciante no R, e estou fazendo umas análises de regressões nonlinear , LL, testando doses de um herbicida no comprimento da raiz do arroz (cm).
Como fiz dois experimentos, estou vendo se ambos são iguais para discutir em conjunto. Já vi que não as curvas fitted não são iguais, tanto utilizando os dados originais (cm) ou usando dados transf. em %. Ou seja, o p <0.05.
Tentei fazer uma transformação com box-cox, mas não fui"eficiente" em usar o R. Poderia me dar uma dica como usar esta ferramenta?
Tenho no file uma cultivar, 7 doses, root , divididos em dois estudos O que fiz inicialmente foi testar a logistica LL.3:
cult<drm(root~dose, data=cultivar, fct=LL.3())
mas com os dados originais ou em % não consigo P>0.05. POr isto gostaria de testar box-cox.
obrigado..Andre
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
-- Qual o objetivo de você está querendo utilizar um box-cox. Os pressupostos do modelo não foram atendidos? Se seu objetivo é comparar os modelos você deve fazer um teste de identidade de modelos. Não está claro sua pergunta. Recomendo que você coloque o CMR para que possa obter melhores respostas. Leia o manual do R-BR para saber como fazer isto. Eu geralmente utilizo a função boxcox do pacote MASS. boxcox(root ~ dose , data = cultivar, plotit=T) Se vc quiser ver os valores ao invés do gráfico, altere o plotit para F Fernando Antônio de Souza Zootecnista, Dsc. Nutricão Animal, UFMG lattes: http://lattes.cnpq.br/6519538815038307
Bom dia a todos, Tenho duas curvas (dose-resposta) e usei o modelo de Brain-Cousens para fittar as curvas. E foram criados os parâmetros b, d, e (BC.4), o lower limit é zero. Mas quando tentei gerar os ED para 10, 25, 50 e 95%, no pct drc, através de:
ED(clt, c(10, 25, 50, 95), (interval="delta")
nao aceita. Pela modelo LL.3 deu tudo certo. Li em um paper que realmente não seria possível (se li certo) fazer por BC.4, mas tenho um colega que fez e consegui nos seus dados gerar os EDs e depois tentamos com os meus e não conseguimos.. Alguém pode dar um help? obrigado........André
Seu exemplo não é reproduzível, você não forneceu o modelo nem o que significa LL.3, BC, BC.4, ED, a, b e c, etc. À 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 skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
participantes (3)
-
Andre Andres -
Fernando A. Souza -
Walmes Zeviani