[R-br] construção de grafico complexo...

[Ayuni Sena]

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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110419/c9cb9a97/attachment-0001.html>