similar to: sorting Trellis panels

Hello all, I am trying to use xyplot to give a six panel plot and to highlight only points (in any panel) that meet a certain criterion. With the plot command I would do something like: plot.default(filein$Site,filein$circ.conc) points(filein$Site,filein$circ.conc,type="p", pch=ifelse(filein$p.value<5e-02,19,21)) I had thought I could just stick in the pch line from above into

Hi, I would like to compare two models in R with the same dependant variable but different predictors (two different types of frequency and reaction times as RT). An editor told me to have a look at Lorch and Myers 1990. Lorch and Myers use the following technique: 1) they perform regressions on individual subjects' data 2) they extract the beta weights 3) they run t-test on these beta

How does one go about doing a repeated measure regression? The documentation I have on it (Lorch & Myers 1990) says to use linear / (subj x linear) to get your F. However, if I put subject into glm or lm I can't get back a straight error term because it assumes (rightly) that subject is a nominal predictor of some sort. In looking at LME it seems like it just does the right thing

Hello, I am trying to run a mixed effects nested ANOVA but none of my codes are giving me any meaningful results and I am not sure what I am doing wrong. I am a new user on R and would appreciate some help. The experimental design is that I have some frogs that have been exposed to three acoustic Treatments and I am measuring neural activity (egr), in 12 brain regions. Some frogs also called

It was suggested to me that I forward this message to you. ---------- Forwarded Message ---------- Subject: privsep in ssh Date: Fri, 19 Sep 2003 12:22 From: Russell Coker <russell at coker.com.au> To: SE Linux <selinux at tycho.nsa.gov> Cc: Colin Watson <cjwatson at debian.org> #ifdef DISABLE_FD_PASSING if (1) { #else if (authctxt->pw->pw_uid == 0 ||

i'm trying to build a little summary table for the contents of a data frame. t<-sapply(macro, data.class) c<-sapply(macro, length) m<-sapply(macro, mean, na.rm=T, digits=2) cbind(type=t, n=c , mean=m) I want to drop the variables that are factors so I can include -max- and -min- in my table. -macro- contacts the data--how do I drop the variables according to their data.class

Hello, I have conducted an SEM in which the resultant standardized path coefficients are much higher than would be expected from the raw correlation matrix. To explore further, I stripped the model down to a simple bivariate relationship between two variables (NDVI, and species richness), where it's my understanding that the SEM's standardized path coefficient should equal the correlation

Hi, I'm having trouble with an IRQ conflict I can't figure out how to resolve. I tried KDE (which is great), but this required the NVIDIA driver. However, this driver conflicts with my modem. If the driver is loaded, then the system locks up when accessing the modem. If I remove the NVIDIA driver from /boot/loader.conf, modem works fine. I tried moving the modem to different slots, but

If libxl_create_logfile fails for some reason then libxl__create_qemu_logfile previously just carried on and dereferenced the uninitialised logfile. Check for the error from libxl_create_logfile, which has already logged for us. This was reported as Debian bug #784880. Reported-by: Russell Coker <russell at coker.com.au> Signed-off-by: Ian Campbell <ian.campbell at citrix.com> Cc:

Is it expected that running btrfsck more than once will keep reporting errors? Below is the end of a btrfsck output when run the second time. backpointer mismatch on [111942471680 32768] owner ref check failed [111942471680 32768] ref mismatch on [111942504448 40960] extent item 1, found 0 Incorrect local backref count on 111942504448 root 5 owner 160739 offset 3440640 found 0 wanted 1 back

Dear group: I wrote a code to iterate a non-linear fit with a set of data. The entire process didn't implemented to the end because an error message, "singular gradient". I knew that some sub-sets (columns) do not fit my formula well and may result in parameters going to infinity. It is pretty hard to remove those sub sets before running the code since that will take a lot of

Greetings- I'm venturing into a new (for me) area of using R: some presentation graphics based on R data. I'd like to compare two distributions of ordinal variables in relatively small (N=203 and 207) subsets of an R dataset. I can get a reasonably good picture of the distributions' differences with: x11() plot(table(hcd.df$auth.sum[hcd.df$datecat==1]),

Hello, I have the following model: lm.7 <- lm(Y ~ F + C1 + C2 , data = EM4) F is a 4-level factor, the rest are covariates centered at their mean (Y is a two-column matrix). I have tried to find functions to give the model-adjusted means (adjusted at the covariates'means) and their standard deviations for each. (That is, what I believe is called in SAS "least square or LS-means,

Hello, I don't know if anyone gets familiar with the llvm Chapel front end, it's version 3.3 and comes within Chapel source package. I'm trying to parse the .bc code generated from a simple helloworld chapel program. Just run a basic Hello pass on the code (default one in the llvm source code) that prints out all the function names in the module. However, it prints a bunch of

Greetings- I am in the process of constructing a nonlinear model using nlme. The model is attempting to fit a nested data structure from some public-opinion data (the data are from individuals nested within organizations). The question I have is fairly simple (I hope). The data were collected in two stages: a 15,000-subject randomly-sampled telephone interview (the SCREEN), with 2,517 subjects

Greetings- A relatively small request, but here goes. I have some crosstabs, generated with table(), in R: > table(pg.df$grouptype,pg.df$age) 0 1 2 3 4 B 1 3 6 8 3 C 1 3 4 10 7 P 0 2 9 9 8 S 2 9 4 5 16 U 0 15 15 3 0 X 0 0 0 0 0 I would like to export these in latex format using xtable(). But it appears there's no appropriate method for

Greetings- I'm working through Pinheiro and Bates' _Mixed Effects Models in S and S-Plus_ using R (1.3.1 for linux). On page 13 (okay, so I haven't got that far :)) is: plot.design( ergoStool) which returns on my system: > plot.design(ergoStool) Error: couldn't find function "plot.design" any ideas? Thanks.

Greetings. I'm planning to make a point in a meeting tomorrow and want to make sure I've got my ducks in a row. I've got four distinct vectors: years, fec, gss, and nes. Years is simply seq(1972,1996,by=4), and the other three are measures of voter participation in US elections in those years. fec is the Federal Election Commission's numbers, gss and nes are the General Social

Attempting to install the most recent R on a Windows/XP machine, I get what appears to be a normal (completely default) installation. Then when I attempt to load rgui I get the fatal error message: Invalid HOMEDRIVE. Any explanation of the error or suggestions of fixes to files or registry would be appreciated. -- Jim Stimson Email: <jstimson at email.unc.edu> Web:

I'm working on a package for R to create graphs for modeling interactions between latent variables in SEM. At the moment, I have one function which performs the necessary steps to prepare the data for graphing, and then separate graphing functions depending on what graph type (or types) the user wants. To get the prepped data from the setup function to the graphing functions, I