Pessoal, lá vai mais uma mega duvida.... Eu preciso gerar um gráfico muito parecido com a figura disponivel no link: http://www.datafilehost.com/download-f1ba4ceb.html Existe apenas uma pequena diferença: No meu caso, o eixo Y é repesentado pela variável N o eixo X pela variavel D Foram ajustados os seguintes modelos: 1 - N~SSasymp(D) #Regressao assintotica do livro do Pinheiro e Bates (2000) - Mixed effect models 2 - D ~ b0*N^b1*hdc^b2 3 - V ~ b3*D^b4*hdc^b5*N^b6 A partir destes modelos, quero plotar no gráfico as isolinhas que correspondem a: 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. 2 - a partir do modelo 2, obtemos a seguinte equação: N = (D/(b0*hdc^b2))^(1/b1) 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) 3 = a partir do modelo 3, obtemos a seguinte equação: N = (V/(b3*D^b4*hdc^b5))^(1/b6) 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) os passos nº 1 e nº 2 eu consegui fazer com auxilio dos comandos "seq" e "lines", da seguinte forma: #passo 1: N100<-779.36+(26947.92)*(exp(-exp(-1.29)*D)) N55<-428.648+(14821.356)*(exp(-exp(-1.29)*D)) N40<-311.744+(10779.168)*(exp(-exp(-1.29)*D)) N15<-116.904+(4042.188)*(exp(-exp(-1.29)*D)) plot(D, N, xlab=" Dq (cm)", ylab="Número de plantas por hectare (N/ha)", ylim=c(0,5000)) lines(dq,n100,type="l", lty=1, lwd=1) lines(dq,n55,type="l", lty=2, lwd=1) lines(dq,n40,type="l", lty=1, lwd=1) lines(dq,n15,type="l", lty=2, lwd=1) #passo 2: b0=25.360699 b1=-0.303046 b2=0.51502 b3=0.0000377 b4=2.2286023 b5=0.7798658 b6=0.9907735 D<-seq(from=5.5, to=21, lenght.out=1061) hdc<-seq(from=7, to=31,5, lenght.out=1061) N1<-(D/(b0*15^b2))^(1/b1) N2<-(D/(b0*20^b2))^(1/b1) #etc etc etc plot(DQ, NHA, ylab="Densidade de plantas por hectare (N/ha)", xlab="Dq (cm)") lines(D, N1,type="l", lty=1, lwd=1) lines(D, N2,type="l", lty=1, lwd=1) #etc etc etc Com auxílio desta rotina eu nao consegui fazer o passo 3.... Alguem teria uma sugestão melhor???????????????????????????? Agradeço a atenção de todos! Ayuni