Re: [R] Intercepts is coming as Zero in the Mixed Models

2011-11-07 Thread Ben Bolker
arunkumar akpbond007 at gmail.com writes:
 
 I'm getting the intercepts of  the Random effects as 0. Please help me to
 understand why this is coming Zero
 
 This is my R code
 
 Data- read.csv(C:/FE and RE.csv)
 Formula=Y~X2+X3+X4 + (1|State) + (0+X5|State)
 fit=lmer(formula=Formula,data=Data)
 ranef(fit).
 

  This question is more suited for the r-sig-mixed-models mailing
list ...

   You are getting an estimate of zero variance because lmer
is computing that as the best estimate.  The reason is that it's
really completely impractical to try to estimate the variance among
levels of a factor with only two levels.  There is more discussion
of this issue at

http://glmm.wikidot.com/faq#fixed_vs_random

  Thanks for including your data so we could see the problem.

__
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.


[R] Intercepts is coming as Zero in the Mixed Models

2011-11-06 Thread arunkumar1111
Hi

I'm getting the intercepts of  the Random effects as 0. Please help me to
understand why this is coming Zero

This is my R code

Data- read.csv(C:/FE and RE.csv)
Formula=Y~X2+X3+X4 + (1|State) + (0+X5|State)
fit=lmer(formula=Formula,data=Data)
ranef(fit).

My sample Data

State Year Y X2 X3 X4 X5 X6
S2 1960 27.8 397.5 42.2 50.7 78.3 65.8
S1 1960 29.9 413.3 38.1 52 79.2 66.9
S2 1961 29.8 439.2 40.3 54 79.2 67.8
S1 1961 30.8 459.7 39.5 55.3 79.2 69.6
S2 1962 31.2 492.9 37.3 54.7 77.4 68.7
S1 1962 33.3 528.6 38.1 63.7 80.2 73.6
S2 1963 35.6 560.3 39.3 69.8 80.4 76.3
S1 1963 36.4 624.6 37.8 65.9 83.9 77.2
S2 1964 36.7 666.4 38.4 64.5 85.5 78.1
S1 1964 38.4 717.8 40.1 70 93.7 84.7
S2 1965 40.4 768.2 38.6 73.2 106.1 93.3
S1 1965 40.3 843.3 39.8 67.8 104.8 89.7
S2 1966 41.8 911.6 39.7 79.1 114 100.7
S1 1966 40.4 931.1 52.1 95.4 124.1 113.5
S2 1967 40.7 1021.5 48.9 94.2 127.6 115.3
S1 1967 40.1 1165.9 58.3 123.5 142.9 136.7
S2 1968 42.7 1349.6 57.9 129.9 143.6 139.2
S1 1968 44.1 1449.4 56.5 117.6 139.2 132
S2 1969 46.7 1575.5 63.7 130.9 165.5 132.1
S1 1969 50.6 1759.1 61.6 129.8 203.3 154.4
S2 1970 50.1 1994.2 58.9 128 219.6 174.9
S1 1970 51.7 2258.1 66.4 141 221.6 180.8
S2 1971 52.9 2478.7 70.4 168.2 232.6 189.4 

--
View this message in context: 
http://r.789695.n4.nabble.com/Intercepts-is-coming-as-Zero-in-the-Mixed-Models-tp3997498p3997498.html
Sent from the R help mailing list archive at Nabble.com.

__
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.