Re: [R] Issue with 95% CI using MASS{confint}

2015-05-28 Thread Jennifer Sweatman
Hi all, Thank you for your responses. Unfortunately, setting level=0.951 or 0.94 did not work for me. I followed your code, Duncan, and some of my tau values are less than 2. I don't know how to interpret this information, however. Can you point me in the direction to a reference?

Re: [R] Issue with 95% CI using MASS{confint}

2015-05-28 Thread Duncan Murdoch
On 28/05/2015 1:11 PM, Jennifer Sweatman wrote: Hi all, Thank you for your responses. Unfortunately, setting level=0.951 or 0.94 did not work for me. I followed your code, Duncan, and some of my tau values are less than 2. I don't know how to interpret this information, however. Can

Re: [R] Issue with 95% CI using MASS{confint}

2015-05-28 Thread sami . toppinen
I encountered the same problem today. It seems to occur only when you use exactly 95 % as the confidence level. For example: confint(g1) Waiting for profiling to be done... 2.5% 97.5% a0 1.257512 1.330881 KP NA 163.862932 confint(g1, level = 0.95) Waiting for profiling

Re: [R] Issue with 95% CI using MASS{confint}

2015-05-28 Thread Duncan Murdoch
On 28/05/2015 9:08 AM, sami.toppi...@kolumbus.fi wrote: I encountered the same problem today. It seems to occur only when you use exactly 95 % as the confidence level. For example: confint(g1) Waiting for profiling to be done... 2.5% 97.5% a0 1.257512 1.330881 KP NA

[R] Issue with 95% CI using MASS{confint}

2015-05-27 Thread Jennifer Sweatman
I am trying to run a nonlinear model looking at seasonal abundance of plants with the following limitations: fit.nls.s-nls(S_a~beta+m*time+alpha*sin(w*time+phi), data= xy, start=list(beta=2.1, m=0, alpha=2, phi = 1, w = 1), upper=list(beta=inf, m=inf,

Re: [R] Issue with 95% CI using MASS{confint}

2015-05-27 Thread Duncan Murdoch
On 27/05/2015 4:08 PM, Jennifer Sweatman wrote: I am trying to run a nonlinear model looking at seasonal abundance of plants with the following limitations: fit.nls.s-nls(S_a~beta+m*time+alpha*sin(w*time+phi), data= xy, start=list(beta=2.1, m=0, alpha=2,