Displaying 2 results from an estimated 2 matches for "pet3".
Did you mean:
pet
2001 Nov 14
2
lme: how to extract the variance components?
Dear all,
Here is the question:
For example, using the "petrol" data offered with R.
pet3.lme<-lme(Y~SG+VP+V10+EP,random=~1|No,data=petrol)
pet3.lme$sigma gives the residual StdDev.
But I can't figure out how to extract the "(intercept) StdDev",
although it is in the print out if I do "summary(pet3.lme)".
In S-plus3.4 , $var.ran is th...
2011 Dec 08
1
partial duplicates of dataframe rows, indexing and removal
...from my dataframe, those rows in which the first
7 columns are duplicated even if subsequent columns make those rows unique.
df<-data.frame(id=rep(c('amy','bob','joe') , each=5),
pet1=sample(LETTERS[1:3],15, replace=T),
pet2=sample(LETTERS[1:3],15, replace=T),
pet3=sample(LETTERS[1:5],15, replace=T))
>df
id pet1 pet2 pet3
1 amy C B A
2 amy B A A
3 amy A A D
4 amy B C A
5 amy C B B
6 bob B A A
7 bob C A C
8 bob C C A
9 bob B C E
10 bob C B C
11 joe...