[R-br] transf Box-Cox no R

[Fernando A. Souza]

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 em 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