[R] Asking for help
Boris Steipe
boris.steipe at utoronto.ca
Wed Apr 19 22:16:10 CEST 2017
What about the advice I sent you this morning?
B.
> On Apr 19, 2017, at 1:49 PM, Boris Steipe <boris.steipe at utoronto.ca> wrote:
>
> No attachments please. Use the dput() function. And keep the correspondence on the list, there are many colleagues who are far more knowledgeable than I am.
>
> Cheers,
> Boris
>
>
>
>
>> On Apr 19, 2017, at 1:47 PM, Santiago Bueno <swbueno at gmail.com> wrote:
>>
>> Thanks a lot Mr. Steipe. I have attached to this email an Excel file with 48 observations. The problem is that when I try to get the prediction plot for the following model:
>> start <- coef (lm(log(Btot)~I(log(dbh**2*haut)),data=dat[dat$Btot>0,]))
>> start[1] <- exp(start[1])
>> names(start) <- c("a","b")
>> model18 <- nls(Btot~a*(dbh**2*haut)**b,data=dat,start=start,weights=1/dat$dbh**4)
>> summary(model18)
>>
>> using this code:
>>
>> with(dat,plot(dbh**2*haut,Btot,xlab="dbh^2*H(cm^2.m)",ylab="Biomass (t)", pch=19))
>> D <- seq(0,180,length=20)
>> H <- seq(0,61,length=20)
>>
>> abline(model18, lwd=3)
>>
>> I get the data points but not the prediction line (see below)
>>
>> <image.png>
>> If you can help me in getting this line, I will be eternaly grateful..
>>
>> Best regards,
>>
>> Santiago Bueno
>>
>> On Wed, Apr 19, 2017 at 12:20 PM, Boris Steipe <boris.steipe at utoronto.ca> wrote:
>> Can you provide a small reproducible example and explain what exactly is going wrong?
>>
>> Just a handful of data points will do.
>>
>>
>> B.
>>
>>
>>
>>> On Apr 18, 2017, at 2:16 PM, Santiago Bueno <swbueno at gmail.com> wrote:
>>>
>>> Thanks Boris, the following is an extract of my data. I have developed
>>> biomass models using codes like:
>>>
>>> start <- coef (lm(log(Btot)~I(log(dbh**2*haut)),data=dat[dat$Btot>0,]))
>>>
>>> start[1] <- exp(start[1])
>>>
>>> names(start) <- c("a","b")
>>>
>>> M1 <- nls(Btot~a*(dbh**2*haut)**b,data=dat,start=start,weights=1/dat$dbh**4)
>>>
>>>
>>> start <- coef(lm(log(Btot)~I(log(dbh))+I(log(haut)),data=dat[dat$Btot>0,]))
>>>
>>> start[1] <- exp(start[1])
>>>
>>> names(start) <- c("a","b1","b2")
>>>
>>> M2 <- nls(Btot~a*dbh**b1*haut**b2,data=dat,start=start,weights=1/dat$dbh**4)
>>>
>>>
>>> Tree No dbh haut Btot
>>> 1 35.00 18.90 0.535
>>> 2 25.00 16.60 0.248
>>> 3 23.00 19.50 0.228
>>> 4 13.50 15.60 0.080
>>> 5 20.00 18.80 0.172
>>> 6 23.00 17.40 0.190
>>> 7 29.00 19.90 0.559
>>> 8 17.60 18.20 0.117
>>> 9 31.00 25.30 0.645
>>> 10 26.00 23.50 0.394
>>> 11 13.00 13.00 0.038
>>> 12 32.00 20.70 0.443
>>> It is my interest to get prediction plots for the models. I have tried to
>>> use the following codes with no success: Let m be one of the fitted models
>>> with dbh as the only entry. To construct a plot of the predictions made by
>>> this model I have tried:
>>> with(dat,plot(dbh,Btot,xlab="Dbh(cm)",ylab="Biomass (t)"))
>>> D <- seq(par("usr")[1],par("usr")[2],length=200)
>>> lines(D,predict(m,newdata=data.frame(dbh=D)),col="red")
>>> For a model m that has dbh and height as entries, I have tried to get its
>>> predictions as follows:
>>> D <- seq(0,180,length=20)
>>> H <- seq(0,61,length=20)
>>> B <- matrix(predict(m,newdata=expand.grid(dbh=D,height=H)),length(D))
>>>
>>> Can someone provide help please!!!
>>>
>>>
>>> Best regards,
>>>
>>> Santiago Bueno
>>>
>>> [[alternative HTML version deleted]]
>>>
>>> ______________________________________________
>>> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
>>> https://stat.ethz.ch/mailman/listinfo/r-help
>>> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
>>> and provide commented, minimal, self-contained, reproducible code.
>>
>>
>> <Biomasa_Total.xls>
>
> On Apr 19, 2017, at 4:13 PM, Ista Zahn <istazahn at gmail.com> wrote:
>
> For help in formulating your question in a way that makes it easy to
> help you see http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example
>
> When asking for help on the R-help mailing list in particular, make
> sure to send plain-text email.
>
> Best,
> Ista
>
> On Wed, Apr 19, 2017 at 3:56 PM, Santiago Bueno <swbueno at gmail.com> wrote:
>> I have tried to ask for help ragarding R software but it seems that I do
>> not know how to send my questions. Can someone "PLEASE" tutor me on how to
>> structure my messages?
>>
>> Best regards,
>>
>> Santiago
>>
>> [[alternative HTML version deleted]]
>>
>> ______________________________________________
>> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
>
> ______________________________________________
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
More information about the R-help
mailing list