similar to: The behavior of match function

It is not a pure R question,but I hope some one can give me advices. I want to use analysis my data with the multilevel model.The data has 2 levels---- the second level has 52 units and each second level unit has 19-23 units.I think the sample size is quite small,but just now I can't make the sample size much bigger.So I want to ask if I use the multilevel model to analysis the data set,will

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

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 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())

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

lapply and rownames will do what you want > class(t) [1] "list" > lapply(t,rownames) [[1]] [1] "ca" "cd" [[2]] [1] "cb" [[3]] [1] "cc" ======= 2005-12-11 01:17:47 伳侜佋佢伬伌佇伵佒佇佇伌伒伬仯伜======= >Hi all, > I have a matrix and its cluster...and i like to extract the row names of each cluster....but using the for loop or some other looping

I have a spss format data file which is about 134M,and I want to read it into R and analysis the data. I plan to split the file into 4 parts and then use read.spss to read it,but during the process,the RGui crashes. and I try to convert it into csv file and use read.csv to import it,but after 1 hours,the process is still not completed. So I want to know if my computer can do this job without

> 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"

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

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

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) +

>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

======= 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

Hi,When I read the source of str,i find these code ----- ## Show further classes // Assume that they do NOT have an own Method -- ## not quite perfect ! (.Class = 'remaining classes', starting with current) cl <- oldClass(object); cl <- cl[cl != "data.frame"] #- not THIS class ----- so I use ?oldClass to try to learn more about oldClass.But after I have reading

It's often heard that the S language is a functional language.But What's the exact meaning of this termology in the context of S language? 2005-10-18 ------ Deparment of Sociology Fudan University Blog:http://sociology.yculblog.com

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

> version _ platform i386-pc-linux-gnu arch i386 os linux-gnu system i386, linux-gnu status major 2 minor 1.0 year 2005 month 04 day 18 language R OS:debian linux # R CMD INSTALL /home/ronggui/Matrix_0.96-3.tar.gz ......... s-3 -L/usr/lib/R/lib -lR /usr/bin/ld: cannot find -lblas-3 collect2: ld returned 1 exit status make: *** [Matrix.so] Error 1

>library(foreign) in the foreign package,you can use read.csv command to read the csv file.you can use excel to open the file and save as cvs file.you should read the manuals first,which tells in details how to import your data into R. before you use the command frome the arules package,you should first type the command >library(arules) ======= 2005-07-26 22:18:48 伳侜佋佢伬伌佇伵佒佇佇伌伒伬仯伜=======

i study the code of function ave,but i can understand one line of the syntax. > ave function (x, ..., FUN = mean) { n <- length(list(...)) if (n) { g <- interaction(...) split(x, g) <- lapply(split(x, g), FUN) } else x[] <- FUN(x) x } my question is : what does "split(x, g) <- lapply(split(x, g), FUN)" mean? thank you! --

entercount an error with demo(scoping). > demo(scoping) demo(scoping) ---- ~~~~~~~ ___snip_____ > ross$balance() Your balance is 120 > try(ross$withdraw(500)) Error in ross$withdraw(500) : You don't have that much money! > version _ platform i486-pc-linux-gnu arch i486 os linux-gnu system i486, linux-gnu status beta major 2 minor