search for: lockwood

...would like to produce a (n x k) model matrix with treatment contrasts for this factor, with rows of NAs placeholding the missing observations. If I use model.matrix() I get back a (nobs x k) matrix. Is there an easy way to get the (n x k) without carrying along a row ID and merging? Thanks. J.R. Lockwood 412-683-2300 x4941 lockwood at rand.org http://www.rand.org/methodology/stat/members/lockwood/

...686, linux-gnu status major 1 minor 5.1 year 2002 month 06 day 17 language R Thanks in advance for your help -- I've learned a ton of statistics and computing on this list. J.R. Lockwood 412-683-2300 x4941 lockwood at rand.org http://www.rand.org/methodology/stat/members/lockwood/ -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or &qu...

Can a single random number be generated in R? I have an exercise that wants to simulate coin tosses, and I cannot seem to find a good example of the use of random number generation in R. Any help? Joshua Gramlich Chicago, IL

Dear Lockwood, As you can see, I'm a beginner... But thank you very much! Sincerely, Tord Quoting "J.R. Lockwood" <lockwood at rand.org>: > how about > > buskartant$buskartant <- sapply( group.list, function(x) > sum(!is.na(unique(x))) ) > > OR > > buskarta...

Greetings, This seems too rudimentary to ask but for the life of me I cannot locate a readily easy method to compute the univariate mode. I know "mode" is not correct and "table" provides a reasonable count but I figured there would be an easy way to extract the value from the table after I do something like: max(table(mydadat$myvar)) unfortunately it only returns the max

...That is, if the data are organized by student, the "Z" matrix in the usual linear mixed model notation has every other row equal to a row of zeros. I am wondering whether there is some way to fit this model using lme(). Thanks in advance for your help and patience. best regards, J.R. Lockwood 412-683-2300 x4941 lockwood at rand.org http://www.rand.org/methodology/stat/members/lockwood/

without having to check the vector element by element? Thanks a lot! Jason ===== Jason G. Liao, Ph.D. Division of Biometrics University of Medicine and Dentistry of New Jersey 335 George Street, Suite 2200 New Brunswick, NJ 08903-2688 phone (732) 235-8611, fax (732) 235-9777 http://www.geocities.com/jg_liao

Hi all; This is my first message to the list, and I've got three "basic" questions: How could I insert comments in a file with commands to be used as source in R? Is it possible to quickly display a window with all the colors available in colors()? How? I'm displaying points, but they overlap, wether points() uses triangles, bullets or whatever. Is it possible to change

Dear listers, I can't find the variance or se of the coefficients in a multilevel model using lme. I want to calculate a Chi square test statistics for the variability of the coefficients across levels. I have a simple 2-level problem, where I want to check weather a certain covariate varies across level 2 units. Pinheiro Bates suggest just looking at the intervals or doing a rather

...)" but is not elegant and of only limited use because of the lack of association with the variable name. My question is whether there is a better way to do this, either with "vars" as defined above or through some more clever definition. Thanks to all for your consideration. J.R. Lockwood 412-683-2300 x4941 lockwood at rand.org http://www.rand.org/methodology/stat/members/lockwood/ -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or...

...,d$y)^2) However, the same is not true of m2. print(summary(m2)$r.squared - cor(fitted(m2),d$y)^2) > R.version _ platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 1 minor 9.0 year 2004 month 04 day 12 language R J.R. Lockwood 412-683-2300 x4941 lockwood at rand.org http://www.rand.org/methodology/stat/members/lockwood/ -------------------- This email message is for the sole use of the intended recipient(s) and may contain privileged information. Any unauthorized review, use, disclosure or distribution is prohibited. I...

It seems that glances is not in EPEL for CentOS 8. Is that coming at some point, or has it been deprecated in favor of some other tool?

I want to subset the dataframe based on certain values in a row. for each row in my dataframe if ANY one value of a particular set of columns satisfies cond append a logical value true at the end of the row else append a false at the end of the row in the end I want to be able to subset the whole data based on the appended true or false value. I could literally code like this, but I think

If you found my subject heading to be confusing then I'm sure you'll enjoy the example I've included below. I find the apply type functions to be wonderful for avoiding loops but when I use them with existing functions, I end up using loops anyway to extract the vectors I want. I would appreciate it if someone could show me how to avoid these loops. Thanks. EXAMPLE:

In general this is what "apply" does. In your example you could use rowMeans(ar) as well, I think. Reid Huntsinger -----Original Message----- From: Olivier Martin [mailto:olivier.martin at inrialpes.fr] Sent: Tuesday, July 02, 2002 10:55 AM To: r-help Subject: [R] mean and array Hi all, I have an array, one say ar, with 3 dimensions. dim(ar) >200 3 4 I would like to

Hi Helpers: I have a factor object that has 314k entries of 39 land cover types. (This object can be coerced to characters neatly should that be easier to work with.) > length(foo) [1] 314482 > foo[1:10] [1] Montane Chaparral Barren Red Fir Red Fir [5] Red Fir Red Fir Red Fir Red Fir [9] Red Fir Red Fir 39 Levels:

Dear all, Sorry to bother you with one more newbie question. I have a dataobject with several hundreds of columns. I want to remove columns with a mean of the column values below a certain value: > a<- c(1,2,3,4,5,6) > b<-c(2,4,6,8,10,12) > c<- c(3,6,9,12,15,18) > test<- as.matrix(cbind(a, b, c)) > mean(a) [1] 3.5 > mean(b) [1] 7 > mean(c) [1] 10.5 Say the

Dear all, I observed that the R process is growing in time (I was working with the same session for two days, and the size reached 300 Mb). Explicit calls to 'gc()' did show a trigger around 90 Mb but did not change the RSS (while I remembered that my linux was nicely behaving before and the process size was shrinking when R was freeing memory...). I ended up saving the session, ending it

Hello R-users, I need to create a vector inserting an 1 after each value of another vector. For example: vec1<-c(2,3,4) I need to create a vector with the values 2,1,3,1,4 Does anyone know how create this vector without loops (vec1 could have 1000 elements) Thank you, Juan -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

I have been building up a script to quickly and easily make CentOS/RHEL and Fedora VMs from kickstart files for a long time, and thought I'd see if anyone else was interested. It's especially useful IMHO if you are working on building kickstarts, because you can fairly rapidly iterate and test. I've got it built as an RPM, so if others think this is useful, I might submit it to