estive a ver o help(abline) mas o unico exemplo que encontro é para o modelo y~x onde posso encontrar um exemplo da função abline() para lm(pthi~dose+grupo) ?Ana
Repetindo o que o Walmes disse: O abline() só funciona para modelos y~x, para mais que isso tem que programar. Como a variável é categórica, não é muito difícil, da para fazer usando dois abline() e especificando o coeficiente angular (que é o mesmo para os dois) e o intercepto (que vai variar). Para deixar eles diferentes, tu pode variar os parâmetros col e/ou lty. Repetindo o que o Walmes disse [2]: Não envie mensagens com titulo "dúvida". 99% dos emails aqui são dúvidas, tente mandar e-mails com assuntos descrevendo o problema. Nova sugestão: Se o problema é o mesmo, responda o último e-mail enviado, assim fica tudo junto para quem for consultar posteriormente. 2012/11/30 <alanarocha@sapo.pt>
estive a ver o help(abline) mas o unico exemplo que encontro é para o modelo y~x
onde posso encontrar um exemplo da função abline() para lm(pthi~dose+grupo) ?Ana ______________________________**_________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-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<http://www.leg.ufpr.br/r-br-guia>) e forneça código mínimo reproduzível.
É mais fácil usar o comando predict() e sobrepor os ajustados sobre o gráfico dos pontos. da <- expand.grid(x=1:10, tr=gl(3,1)) da$y <- as.numeric(da$tr)+0.5*da$x+rnorm(da$x,0,0.1) plot(y~x, da, col=da$tr) m0 <- lm(y~tr+x, da) pred <- expand.grid(tr=levels(da$tr), x=seq(1,10,length=30)) pred$y <- predict(m0, newdata=pred) plot(y~x, da, col=da$tr) by(pred, INDICES=pred$tr, function(i) lines(i$x, i$y, col=i$tr)) Um exemplo com gráfico melhor em https://ridiculas.wordpress.com/2012/04/23/ajuste-simultaneo-de-regressaoano... À 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 ==========================================================================
Nao há sentido nisto uma vez que este modelo ajustado nao define nao uma linha Use predict() para obter linhes com valores ajustados On Fri, 30 Nov 2012, alanarocha@sapo.pt wrote:
estive a ver o help(abline) mas o unico exemplo que encontro é para o modelo y~x
onde posso encontrar um exemplo da função abline() para lm(pthi~dose+grupo) ?Ana _______________________________________________ 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.
Sua duvida é como traçar a curva do modelo ajustado no gráfico? Se for você pode tentar o comando curve()caso contrário desconsidere o exemplo. Ex: vd<-c(0.18, 1.40, 6.85, 12.40, 16.15, 20.45, 25.13, 27.43, 29.23, 32.00, 34.75, 36.85)vi<-c(30, 62, 92, 120, 148, 176, 208, 237, 266, 298, 328, 361)a <- max(vd) + 1reg <- nls(vd ~ A/(1 + B * exp(-C * vi)), start = list(A = a, B = exp(coef(lm((log((a - vd)/vd)) ~ vi))[1]), C = -(coef(lm((log((a - vd)/vd)) ~ vi))[2])))summary(reg)par(family = "serif", cex = 0.8)plot(vd ~ vi, ylab = "vd", xlab = "vi" , las = 1)curve(coef(reg)[1]/(1 + coef(reg)[2] * exp(-coef(reg)[3] *x)), add = T)
Date: Fri, 30 Nov 2012 16:02:23 +0000 From: alanarocha@sapo.pt To: R-br@listas.c3sl.ufpr.br Subject: [R-br] dúvida
estive a ver o help(abline) mas o unico exemplo que encontro é para o modelo y~x
onde posso encontrar um exemplo da função abline() para lm(pthi~dose+grupo) ?Ana _______________________________________________ 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.
participantes (5)
-
alanarocha＠sapo.pt -
Paulo Justiniano -
Rodrigo Coster -
tiago souza marçal -
Walmes Zeviani