Re: [R] Re : Regarding lme functions

[ratna ghosal]

If I give the following command:
 
myfit<-lme(...)
intervals(myfit)
 then its returning the following error
 
"Cannot get confidence intervals on var-cov components:Non-positive definite 
approximate variance-covariance"
 
I would like to extract the confidence interval for the fixed effects......so 
which command to use???
 
Ratna


 
Ratna Ghosal
Research Scholar
Centre for Ecological Sciences
Indian Institute of Sciences
Bangalore-12

--- On Mon, 10/11/08, justin bem <[EMAIL PROTECTED]> wrote:

From: justin bem <[EMAIL PROTECTED]>
Subject: Re : [R] Regarding lme functions
To: [EMAIL PROTECTED]
Cc: "R Maillist" <[EMAIL PROTECTED]>
Date: Monday, 10 November, 2008, 10:35 PM





in nlme package you the intervals function.

myfit<-lme(...)
intervals(myfit)


 Justin BEM
BP 1917 Yaoundé
Tél (237) 99597295
(237) 22040246






De : ratna ghosal <[EMAIL PROTECTED]>
À : r-help@r-project.org
Envoyé le : Lundi, 10 Novembre 2008, 17h51mn 51s
Objet : [R] Regarding lme functions

Hi all,

I recently got accquainted with the R stats and using the 'lme' function
for some of my analyses. Since I am a novice, I am facing a lot of
problems regarding the analyses. 

I have applied the Linear Mixed Effect model for my dataset and would like
to know the commands for extracting the confidence interval for the fixed
effect. 
 
Also would like to request if you can kindly guide me the command
for plotting the observed value vs the value predicted from the mixed
effect model.

Thanking you

Regards

Ratna




Ratna Ghosal
Research Scholar
Centre for Ecological Sciences
Indian Institute of Sciences
Bangalore-12

--- On Mon, 10/11/08, John Poulsen <[EMAIL PROTECTED]> wrote:

From: John Poulsen <[EMAIL PROTECTED]>
Subject: [R] Labeling points with xYplot
To: r-help@r-project.org
Date: Monday, 10 November, 2008, 9:56 PM

Hello R-Users,

I am trying to label points on an xYplot.  I can label the points, but then I
lose the error bars.  I suspect that I should be using panel.xYplot rather than
panel.xyplot, but that didn't work either...

Can anybody help me out?  Example code is below.  First xYplot has error bars,
second tries to label the points and the bars disappear.

Thanks,
John


a=c(1:4)
b=rnorm(4,15,5)
blab=b+1

xYplot(Cbind(b, b-1.5, b+1.5)~a,
      ylim=c(0,max(b)+3))

xYplot(Cbind(b, b-1.5, b+1.5)~a, ylim=c(0,max(b)+3),
      panel=function(x,y,...){
        panel.xyplot(x,y,...)
        ltext(x=list(1,2,3,4),
              y=list(13,14,15,16),
              
labels=list("A","B","C","D"), pos=2)})

______________________________________________
R-help@r-project.org mailing list
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.



      Share files, take polls, and make new friends - all under one roof. G
    [[alternative HTML version deleted]]





      Add more friends to your messenger and enjoy! Go to 
http://messenger.yahoo.com/invite/
        [[alternative HTML version deleted]]


______________________________________________
R-help@r-project.org mailing list
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.