similar to: Incomplete Output from lmer{lme4}

Does anyone know if there is a way to run multivariate linear regression in R? I tried using the lm function (e.g., lm(dv1, dv2~iv1+iv2+iv3), but got error messages. Is my syntax wrong, or do I need a particular package? Thanks, Jeff-- ________________________________________________________ Jeffrey J. Lusk, Ph.D. Postdoctoral Research Associate Department of Forestry &

After reading the help files (?Startup) and using RSiteSearch("defaultPackages"), I have been trying to add several packages to my default startup list using the following code: local({ old<-getOption("defaultPackages") options(defaultPackages=c(old,"lattice","RODBC") }) Subsequently, when I query the default list using:

Hi, I was wondering if I might be able to ask some advice about doing residual plots for the lmer function in the lme4 package. Our group's aim is to find if the expression staining of a particular gene in a sample (or "core") is related to the pathology of the core. To do this, we used the lmer function to perform a logistic mixed model below. I apologise in advance

Hi all, I need all your help on this. I have the next part of code:............e1=x1-mean(x1)e2=x2-mean(x2)n1=length(x1)n2=length(x2N=(n1 + n2)nu2=sum( c( ( x1 -mean(x1) )^2 , ( x2-mean(x2) )^2 ) )/Nss=c(e1,e2) b3= N*sum(ss^4)/ (sum( ss^2)^2)............what do lines 6-8 (mathematical notation)?Also the, what means part?............for(j in 1:B){ ss11=sample(x1, n1, replace=TRUE) ss12=sample(x2,

In Julian Faraway's text on pgs 117-119, he gives a very nice, pretty simple description of how a glm can be thought of as linear model with non constant variance. I just didn't understand one of his statements on the top of 118. To quote : "We can use a similar idea to fit a GLM. Roughly speaking, we want to regress g(y) on X with weights inversely proportional to var(g(y).

Dear all, I hope to know how to compile the R help file in LaTeX format under Windows 2000. The TeX/LaTeX system is TeXLive 2002, and the version of R is 1.8.x. Another question, I hope to get the upright Greek letters, just as these in Prof. Faraway's book "Practical Regression and Anova in R", but I don't know which package(s) and/or font(s) should be used in LaTeX. Any

I am running R 2.6.0 on a PC running MS Windows 2000 Professional. After upgrading to this new version of R, I have been having problems opening scripts: opening certain scripts result in R crashing. After some experimentation, I discovered that it is a problem with the number of subdirectories in the path. For example, if the script is located in c:\Documents and

I am trying to replicate the repeated measures example from Dr.Faraway's book (Extending the linear model with R) as follows: data(vision) vision$npower <- rep(1:4,14) mmod <-lmer(acuity~power+(1|subject)+(1|subject:eye),vision) When I look at the fixed effects p-value, it is missing. Am I missing something here? Thanks. Fixed effects: Estimate Std. Error t value

Dear r helpers, I would like to look at the interaction between two two-level factors, one between and one within participants, after accounting for any variance due to practice (31 trials in each of two blocks) in the task. It seems to require treating practice as a covariate. All the examples I noticed for handling covariates (i.e. ANCOVA, including the ones in Faraway's "Practical

Dear R-experts, I am running R version 2.7.1 on Windows Vista. I have a small dataset which consists of ?chick ID?, ?year (0, 1)?, ?hatching order [HO, defined as first, second and third-hatched chick]?, and the binary outcome of interest ?death (0, 1)?. So a subset of my dataset looks like this on a txt file: y ID Yr HO 1 1 1 First 0 2 1 First 0 3 1 Second 0 4 1 First 1 5 1 First 0 6 1 Third

Hi All, I am running R 2.2.0 on Mac OS 10.4.2, dual G5 processors with 8 Gig RAM. I am running a simulation with lsoda that requires ~378 s to complete one set of time intervals. I need to optimize the parameters, and so need to considerably speed up the simulation. I have tried to figure out how to change the appropriate memory allocation and have search R help and Introductory

Hi, I've got Gluster error: Transport endpoint not connected. It came up twice after trying to rsync 2 TB filesystem over; it reached about 1.8 TB and got the error. Logs on the server side (on reverse time order): [2012-12-15 00:53:24.747934] I [server-helpers.c:629:server_connection_destroy] 0-RedhawkShared-server: destroyed connection of

My apologies in advance for my thickness but I can't seem to solve the following, seemingly simple, data manipulation problem: I have a data frame that contains multiple factors and multiple continuous response variables, but duplicates of some factor combinations. The duplicates contain bad data, so I would like to eliminate the duplicates. I would like to retain the entire rows

Hi Folks, How do I get red bold font on my y axis and black standard font on my x axis? plot(runif(10), ylab="Red, Bold?", xlab="Black, standard?") Any pointers or examples would be great. Thanks! Hank Dr. Hank Stevens, Assistant Professor 338 Pearson Hall Botany Department Miami University Oxford, OH 45056 Office: (513) 529-4206 Lab: (513) 529-4262 FAX: (513) 529-4243

All, In the function below I have 24 individuals and 6 calculations per individual. The 6 calculations are collected each time in a 1:24 loop when calculating "delta". I'd like to collect all 144 = 24*6 calculations in one vector ("delta.patient.comb"). The function works as is via indexing, but is there an easier way to collect the measurements via appendinng the 6

I was trying to help someone who used a spanish keyboard on a PC running Windows. he discovered that he had no tilde key. Does anyone know of any simple work-arounds? Thanks, Hank Dr. Hank Stevens, Associate Professor 338 Pearson Hall Botany Department Miami University Oxford, OH 45056 Office: (513) 529-4206 Lab: (513) 529-4262 FAX: (513) 529-4243 http://www.cas.muohio.edu/~stevenmh/

Hi - I am getting different results when I run the numerical integrator function lsoda (odesolve package) on a Mac and a PC. I am trying to simulating a system of 10 ODE's with two exogenous pulsed inputs to the system, and have had reasonably good success with many model parameter sets. Under some parameter sets, however, the simulations fail on the Mac (see error message below). The

R v. 1.7.1, Windows 2000. A particular journal wants me to provide scatter plots with no box, but with axes that meet in the lower left corner. It seems as though there must be an easy way of doing this, but my reading the help on plot.default, axis, and box have not provided any clues. I would be most appreciative of any feedback. Thank you, Hank Stevens Dr. Martin Henry H. Stevens,

How do I cite a package (not R itself - I know how to do that)? Any thoughts or links? Many thanks in advance? Hank Stevens Dr. Martin Henry H. Stevens, Assistant Professor 338 Pearson Hall Botany Department Miami University Oxford, OH 45056 Office: (513) 529-4206 Lab: (513) 529-4262 FAX: (513) 529-4243 http://www.cas.muohio.edu/botany/bot/henry.html http://www.muohio.edu/ecology/

Using R v. 1.7.0 on Windows 2000 I would like to plot the fitted values of a model as a function of a continuous covariate, augmented with data (e.g., augPred) grouping by combinations of fixed effects. I have not been able to use augPred effectively, and am wondering if it does not handle unbalanced data (3 out of 192 missing). I include below the model and an xyplot that almost does the