<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:times new roman, new york, times, serif;font-size:12pt"><DIV>Pessoal, lá vai mais uma mega duvida....</DIV>
<DIV style="FONT-FAMILY: times new roman, new york, times, serif; FONT-SIZE: 12pt">
<DIV style="FONT-FAMILY: times new roman, new york, times, serif; FONT-SIZE: 12pt">
<DIV style="FONT-FAMILY: times new roman, new york, times, serif; FONT-SIZE: 12pt">
<DIV> </DIV>
<DIV>Eu preciso gerar um gráfico muito parecido com a figura disponivel no link: <A href="http://www.datafilehost.com/download-f1ba4ceb.html" rel=nofollow target=_blank>http://www.datafilehost.com/download-f1ba4ceb.html</A></DIV>
<DIV> </DIV>
<DIV>Existe apenas uma pequena diferença:</DIV>
<DIV>No meu caso, o eixo Y é repesentado pela variável N</DIV>
<DIV>o eixo X pela variavel D</DIV>
<DIV>Foram ajustados os seguintes modelos:</DIV>
<DIV>1 - N~SSasymp(D) #Regressao assintotica do livro do Pinheiro e Bates (2000) - Mixed effect models</DIV>
<DIV>2 - D ~ b0*N^b1*hdc^b2</DIV>
<DIV>3 - V ~ b3*D^b4*hdc^b5*N^b6</DIV>
<DIV> </DIV>
<DIV>A partir destes modelos, quero plotar no gráfico as isolinhas que correspondem a:</DIV>
<DIV>1 - 100%, 55%, 40% e 15% da curva N~SSasymp(Dq) - isso vai gerar apenas um deslocamento da curva para baixo, mantendo o mesmo comportamento.</DIV>
<DIV>2 - a partir do modelo 2, obtemos a seguinte equação:</DIV>
<DIV>N = (D/(b0*hdc^b2))^(1/b1)<BR>e assim, quero plotar no grafico as isolinhas correspondente a hdc=5, hdc=10, hdc=15, hdc=20, e assim por diante... (existirá o mesmo paralelismo entre as isolinhas relatado no passo acima)</DIV>
<DIV>3 = a partir do modelo 3, obtemos a seguinte equação:</DIV>
<DIV>N = (V/(b3*D^b4*hdc^b5))^(1/b6)<BR>e assim, quero plotar as isolinhas correspondentes a V=50, V=100, V=150, e assim por diante.... (existirá o mesmo paralelismo entre as isolinhas relatado no passo acima)</DIV>
<DIV> </DIV>
<DIV>os passos nº 1 e nº 2 eu consegui fazer com auxilio dos comandos "seq" e "lines", da seguinte forma:</DIV>
<DIV>#passo 1:</DIV>
<DIV>N100<-779.36+(26947.92)*(exp(-exp(-1.29)*D))<BR>N55<-428.648+(14821.356)*(exp(-exp(-1.29)*D))<BR>N40<-311.744+(10779.168)*(exp(-exp(-1.29)*D))<BR>N15<-116.904+(4042.188)*(exp(-exp(-1.29)*D))<BR>plot(D, N, xlab=" Dq (cm)", ylab="Número de plantas por hectare (N/ha)", ylim=c(0,5000))<BR>lines(dq,n100,type="l", lty=1, lwd=1)<BR>lines(dq,n55,type="l", lty=2, lwd=1)<BR>lines(dq,n40,type="l", lty=1, lwd=1)<BR>lines(dq,n15,type="l", lty=2, lwd=1)</DIV>
<DIV> </DIV>
<DIV>#passo 2:</DIV>
<DIV>b0=25.360699<BR>b1=-0.303046<BR>b2=0.51502<BR>b3=0.0000377<BR>b4=2.2286023<BR>b5=0.7798658<BR>b6=0.9907735<BR>D<-seq(from=5.5, to=21, lenght.out=1061)<BR>hdc<-seq(from=7, to=31,5, lenght.out=1061)<BR>N1<-(D/(b0*15^b2))^(1/b1)<BR>N2<-(D/(b0*20^b2))^(1/b1) #etc etc etc<BR>plot(DQ, NHA, ylab="Densidade de plantas por hectare (N/ha)", xlab="Dq (cm)")<BR>lines(D, N1,type="l", lty=1, lwd=1)<BR>lines(D, N2,type="l", lty=1, lwd=1) #etc etc etc</DIV>
<DIV> </DIV>
<DIV>Com auxílio desta rotina eu nao consegui fazer o passo 3.... Alguem teria uma sugestão melhor????????????????????????????</DIV>
<DIV> </DIV>
<DIV>Agradeço a atenção de todos!</DIV>
<DIV> </DIV>
<DIV>Ayuni</DIV>
<DIV> </DIV></DIV></DIV></DIV></div></body></html>