similar to: lme intervals

Dear r-helpers, I can run the examples on the mcsamp help page. For example: **************************************** > M1 <- lmer (y1 ~ x + (1|group)) > (M1.sim <- mcsamp (M1)) fit using lmer, 3 chains, each with 1000 iterations (first 500 discarded) n.sims = 1500 iterations saved mean sd 2.5% 25% 50% 75% 97.5% Rhat n.eff beta.(Intercept)

Hi, Can someone please offer a procedure for going from CIs produced by intervals.lme() to fixed-effects response CIs. Here's a simple example: library(mlmRev) library(nlme) hsb.lme <- lme(mAch ~ minrty * sector, random = ~ 1 | cses, Hsb82) (intervals(hsb.lme)) (hsb.new <- data.frame minrty = rep(c('No', 'Yes'), 2), sector = rep(c('Public',

Dear r-helpers, I would very much appreciate help with the following problem: The following command (in a .Rnw file) latex(anova(e7.lmer3, e7.lmer4), file = 'e7lmer34.tex', rowname = c ('nonlinear', 'linear'), longtable = FALSE, dcolumn = T, booktabs = T, table.env = F) produces the following output after running Sweave: % latex.default(anova(e7.lmer1, e7.lmer2),

I am using xYplot to plot data and CIs. How do I add several lines to the figure? _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA 22904-4400 Parcels: Room 102 Gilmer Hall McCormick Road Charlottesville, VA 22903 Office: B011 +1-434-982-4729 Lab: B019

In the course of exploring response prediction, I stumbled upon a small discrepancy between the CIs produced by predict.lm() and all.effects() require(mlmRev) require(effects) hsb.lm <- lm(mAch ~ minrty * sector, Hsb82) hsb.new <- data.frame( minrty = rep(c('No', 'Yes'), 2), sector = rep(c('Public', 'Catholic'), each = 2)) hsb.eff <-

Dear R-helpers, When I try > anova(txtE2.lme, txtE2.lme1) Model df AIC BIC logLik Test L.Ratio p-value txtE2.lme 1 10 8590 8638 -4285 txtE2.lme1 2 7 8591 8624 -4288 1 vs 2 6.79 0.0789 > latex(anova(txtE2.lme, txtE2.lme1)) Error: object "n.group" not found I don't even see n.group as one of the arguments of latex() I checked to see >

I would like to suppress the text 'Scatter Plot Matrix' that appears under the plot. Could someone please suggest how? _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA 22904-4400 Parcels: Room 102 Gilmer Hall McCormick Road Charlottesville, VA 22903 Office: B011

The following plot is a first approximation to what I need: *********************************** mu1 <- 0 mu2 <- 5 s <- 1 x <- seq(-2.5, 7.5, length = 41) y <- seq(-2.5, 2.5, length = 41) f <- function(x,y){ term1 <- 1/(2*pi*sqrt(s*s)) term2 <- -1/2 term3 <- (x - mu1)^2/s term4 <- (y - mu1)^2/s term5 <- (x - mu2)^2/s term1*(.5 * exp(term2*(term3 + term4)) + .5 *

We have a wireless network that uses 802.1x authentication, in which domain joined computers use their machine credentials to connect. Windows machines do this automatically, and until recently Linux computers could join using wicd, wpa-supplicant, and a simple script that would retrieve the machine password with tdbdump. ( specifically tdbdump -k SECRETS/MACHINE_PASSWORD/DOMAIN

> (fm1 <- lmer2(Reaction ~ Days + (Days|Subject), sleepstudy)) Error in as.double(start) : Calloc could not allocate (888475968 of 4) memory ************************* > sessionInfo() R version 2.4.1 (2006-12-18) powerpc-apple-darwin8.8.0 locale: C attached base packages: [1] "grid" "datasets" "stats" "graphics" "grDevices"

excuse for to find a desperate solution... but i'm boried to spent hours... ;) not in asterisk !!! ;) i use the festial-script.pl of Donny Kavanagh... but i want to change the language that festival uses, depending on a variable for the callerid. English/Spanish How i can tell the script the correct voice that festival needs to use?? like --language spanish --language english ...in normal

Dear R People: I'm running from R to a grid computer and getting some unusual results with the system command: > e7 <- system("globus-job-run xxxxx /bin/sh -c 'cd $OSG_APP/engage;chmod 777 oops'",intern=TRUE) chmod: changing permissions of `oops': Operation not permitted > e7 character(0) > I xx'ed out the site name. Anyhow, I would expect that the e7

Dear Dieter, sorry for not being more specific. I would like to use R to get a prediction (with standard error) of the response in a mixed model at selected values of the fixed-effects factors. Hence, in a mixed model, say, for response body size with, say, fixed factors sex and age, I would like to get a prediction of size for each sex and at selected ages such as 5, 10, 15; and I want a SE

This is a follow-up mail of "the effect of blocking on the size of confidence intervals - analysis using aov". In both mails I pursue the idea of using blocking factors in order to reduce the width of confidence intervals. My dataset comprises, a quantitative response variable, namely: "response", and three categorical eplanatory variables, namely: "method",

We are communicating with Asterisk via AMI. Running Asterisk version 13.18.5 on an Ubuntu box. If you look at the event response, the Result field is filled with random characters. I'm not sure what to do because that is a completely random result. It makes no sense. We send the following command to asterisk via AMI Action: AGI ActionID: C44415 Channel: SIP/192.168.40.105-00001338

Dear R-helpers, How do I transform v <- c('insd-otsd', 'sppr-unsp') into c('insd--otsd', 'sppr--unsp') ? _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA 22904-4400 Parcels: Room 102 Gilmer Hall McCormick Road Charlottesville, VA 22903

Dear r-helpers, Does anyone have a straightforward example of putting together three unrelated (expect for a common y-axis) xyplot() figures in what would be in base graphics a par(mfrow = c(1, 3)) arrangement? _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA 22904-4400 Parcels: Room 102

Hi, Does anyone know how to do a coord_trans() in which the y-axis is tranformed into (for example) -1000/y? Thanks, _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA 22904-4400 Parcels: Room 102 Gilmer Hall McCormick Road Charlottesville, VA 22903 Office: B011

I am merging two dataframes using a relational key (incident number and incident year), but not all the records match up. I want to be able to review only the records that cannot be merged for each individual dataframe (essentially trying to select records from one dataframe using a multi-value relational key from the other dataframe). The following code shows what I am trying to do. The final

Dear Friends, Here is an algorithm for finding concave hulls: http://get.dsi.uminho.pt/local/ Has anyone implemented such an algorithm in R? RSiteSearch('concave hull') didn't reveal one (I think). _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology Postal Address: P.O.Box 400400, Charlottesville, VA 22904-4400 Express Parcels