similar to: multilevel models and sample size

data from multilevel units,first sample the class ,and then the student in calss.following is the 2-level model. and the level-1 model deals with the student,and the level-2 model deals with the class level the students belong to. Level-1 Model Y = B0 + B1*(ZLEAD) + B2*(ZBUL) + B3*(ZSHY) + R Level-2 Model B0 = G00 + U0 B1 = G10 + G11*(ZWARMT) + U1 B2 = G20 + G21*(ZWARMT) + G22*(ZABLET) +

the data file has such structure: 1992 6245 49 . . 20 1 0 0 8.739536 0 . . . . . . . . "alabama" . 0 . 1993 7677 58 . . 15 1 0 0

The green book tells:"The basic technique is classic :keep it simple ."A long ,complicated expression or function is less fravorable than" a relatively small computations that combines calls to a few other functions to perform its tasks." But I don't get the point totally.Can anyone give me an example to make me understand this rules totally? ps: Is it mean that f1 is

Hi, Given a frame with calendar date's: "2005-07-01", "2005-07-02","2005-07-03","2005-07-04","2005-07-05",etc. I want to extract the following from these dates: week number month number year number Any ideas how to accomplish this? Many thanks. Regards, Richard

it takes a long time to load the lme4 package.anyone else encounter this problem? > system.time(library(lme4)) 佋佖伻侂佇佽佉仾伒伳伋佁伡伃伆侟仯伜Matrix 佋佖伻侂佇佽佉仾伒伳伋佁伡伃伆侟仯伜lattice [1] 19.90 0.30 25.56 NA NA > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status Patched major 2 minor

>On Wed, 3 Aug 2005, Bernd Weiss wrote: > >> I am trying to replicate some multilevel models with binary outcomes >> using R's "lmer" and "glmmPQL" and Stata's gllmm, respectively. > >That's not going to happen as they are not using the same criteria. the glmmPQL and lmer both use the PQL method to do it ,so can we get the same result by

> library(foreign) > da<-read.dta(file.choose()) > da startdat starttim enddate endtime days hoursmin secused 1 2005-01-11 2 2005-12-15 20.19 NA NA 9 > attributes(da) $datalabel [1] "Example of use of date and time functions" $time.stamp [1] "15 伿伄伓侢佋伮 2005 20" $names [1] "startdat" "starttim"

======= 2005-12-06 22:16:17 伳侜佋佢伬伌佇伵佒佇佇伌伒伬仯伜======= >Martin Maechler a 侀crit : > >> please, please, these trailing ";" are *so* ugly. >> This is GNU S, not C (or matlab) ! >> >> but I'll be happy already if you could >> drop these ugly empty statements at the end of your lines... > >May I disagree ? >I find missing ";" at

the help page says: 'is.list' returns 'TRUE' iff its argument is a 'list' _or_ a 'pairlist' of 'length' > 0, whereas 'is.pairlist' only returns 'TRUE' in the latter case. does the "latter case" mean a 'pairlist' of 'length' > 0? but > is.pairlist(pairlist()) [1] TRUE > length(pairlist())

I use lmer to fit a mixed effect model.It give some warnings.what does this warnings mean? and what should I do? > (fm2.mlm <- lmer(qd ~ edu + jiankang + peixun +hunyin + cadcj + age + age2 + sex + dangyuan + Comp.1 + Comp.2+trust.cz1 +(trust.cz1|commid), data = individual,na.action = "na.exclude",family="quasibinomial")) Generalized linear mixed model fit using PQL

In this case the environment contains the variables local to the function, and its enclosure is the environment of the enclosing function.(R-lang:p11) I want to know if the "enclosing function" means the closure of the function? for example ,if I call function mean(),and the create an environment,say e1,then the enclosure of e1 is the namespace of base package? Right? Thank you!

I want to get the logLik to calculate McFadden.R2 ,ML.R2 and Cragg.Uhler.R2, but the value from multinom does not have logLik.So my quetion is : is logLik meaningful to multinomial logistic model from multinom?If it does, how can I get it? Thank you! ps: I konw VGAM has function to get the multinomial logistic model with logLik, but I prefer use the function from "official" R

Hi, I wonder how well SPSS and R communicate, because I need SPSS but would like to do some data manipulations in R. However I am very afraid of never ending import-export-complications ? especially with all those labels and extra information my SPSS files contain. My data come from SPSS and have to be exported to SPSS again (because I need to produce special output tables, which would be hard to

Dear R-list, I have a dataset, say (the real dataset is 20 columns,110200 rows). > my.reducedID V1 V2 V3 V4 V5 V6 V7 V8 V9 [1,] 1 0 0 1 14 3 1 0 2 [2,] 2 0 0 1 14 3 1 0 2 [3,] 0 1 0 1 14 2 1 0 2 [4,] 0 0 1 1 14 3 1 0 2 [5,] 0 1 1 0 14 2 1 0 2 [6,] 0 0 0 1 14 3 1 0 2 [7,] 0 0 0 1 0 3 1 0 2 [8,] 0

x<-1 f<-function(){ x<-3 eval(substitute(x+y,list(y=10))) } f() #13 x<-1 f<-function(){ x<-3 eval(substitute(x+y,list(y=10)), envir = sys.frame(sys.parent())) } f() #11 x<-1 f<-function(){ x<-3 eval.parent(substitute(x+y,list(y=10))) } f()#11 the help page says: "If 'envir' is not specified, then 'sys.frame(sys.parent())', the

> x<-1:10 > y<-x+1e-20 > x [1] 1 2 3 4 5 6 7 8 9 10 > y [1] 1 2 3 4 5 6 7 8 9 10 > identical(x,y) [1] FALSE > match(x,y) [1] 1 2 3 4 5 6 7 8 9 10 What's the principle the function use to determine if x match y? Thank you! 2005-10-21 ------ Deparment of Sociology Fudan University My new mail addres is ronggui.huang at gmail.com

Hi, I'm using R for some arima models. In the past I used for arima models Rats and Tsp. Using the R arima function, I get only the statistics sigma^2 and log likelihood; with Rats and Tsp it is possible to obtain more statistics, such as R, R square, Durbin Watson, standard error, etc. Is it possible using R to have the statistics mentioned? My best regards,

Hi guys, I have a large number of columns of a dataframe that I want to apply a common factorization to; the columns are all numeric and the factorization collapses some of these values into common groups (labels). How can I do this systematically? Is there an analog to sapply that can pass columns as vectors to FUN? Thanks in advance, Gregory Gentlemen ---------------------------------

Hello, I'm currently using a tool that provides a Loess fitting, but I obtained results that are slightly different from those provided by R implementation of the Loess. That's why I would like to know if you could give me a source (bibliography or web) that explains in a clear way each step of the algorithm, with the possible options to choose, etc.. in order for me to understand those

My question is when the object argument of NexthMethod be used? In the following example, weather object argument is used will not affects the result. ### foo=function(x) {UseMethod("foo")} foo.cls1=function(x) { x=x+1;class(x)<-"ncls" NextMethod() } foo.ncls=function(x) { cat("ncls\n") } foo.cls2=function(x) { cat("cls2\n");print(x) }