search for: anicca

...lmer is under heavy develpoment, maybe something changed or is solved differently. Thus my question: While encountering heavy heteroscedasticity within data, lmer is not the right application to use, but use instead lme? Thanks in advance for a short statement, best , leo -- email: leog at anicca-vijja.de www: http://www.anicca-vijja.de/

Dear list members, the following hlm was constructed: hlm <- groupedData(laut ~ design | grpzugeh, data = imp.not.I) the grouped data object is located at and can be downloaded: www.anicca-vijja.de/lg/hlm_example.Rdata The following works: library(nlme) summary( fitlme <- lme(hlm) ) with output: ... AIC BIC logLik 425.3768 465.6087 -197.6884 Random effects: Formula: ~design | grpzugeh Structure: General positive-definite StdDev Corr...

...owing draws each plot without waiting for a signal (mouse-klick), so par() does not work for that: library(nlme) for(i in 1:3) { fitlme <- lme(Orthodont) par(ask=TRUE) # does not work with trellis.... print( plot(augPred(fitlme)) ) } thanks, leo -- email: leog at anicca-vijja.de www: http://www.anicca-vijja.de/

Dear altogether, I tried local regression with the following data. These data are a part of a bigger dataset for which loess is no problem. However, the plot shows extreme values and by looking into the fits, it reveals very extreme values (up to 20000 !) although the original data are > summary(cbind(x,y)) x y Min. :1.800 Min. :2.000 1st Qu.:2.550

...mular of CohensD (see above) has to be corrected to give respect to the grouping structure. Is there any equivalent measure like eta^2, partial-eta^2, etc.? Can anybody help me with formulas, R code or some references? Thank you very much, thanks in advance, leo g??rtler -- email: leog at anicca-vijja.de www: http://www.anicca-vijja.de/

I just set up the latest dashboard release candidate in hopes of trying out the new inventory service. After modifying the puppet- dashboard config to use the inventory service and after changing my auth.conf I gave it a whirl but ran into some issues.. I initially tried using fact_terminus = yaml and was able to do some queries via curl. I''m currently using storeconfigs so setting up

...ute() and other statistics _and_ I checked the manpages and it does not seem that the calls are wrong. Thus, either it depends on the dataset (but why?) or it is maybe a bug. I appreciate every help, thanks, leo g??rtler <---snip---> library(norm) rngseed(1234) load(url("http://www.anicca-vijja.de/lg/dframe.Rdata")) # load object "dframe" dim(dframe) apply(dframe,2,function(x) sum(is.na(x))) # check how many NAs in the dataset #dframe <- subset(dframe,select=-c(alter,grpzugeh,is1,is4,is6,klassenstufe,mmit,vorai,vorap,voras,vorkf,vorsg,vorvb)) s1 <- prelim.n...

Dear alltogether, two lme's, the data are available at: http://www.anicca-vijja.de/lg/hlm3_nachw.Rdata explanations of the data: nachw = post hox knowledge tests over 6 measure time points (= equally spaced) zeitn = time points (n = 6) subgr = small learning groups (n = 28) gru = 4 different groups = treatment factor levels: time (=zeitn) (n=6) within subject (n=4) w...

Hello alltogether, how is it possible to assign the colnames of a data.frame to a function called by apply, e.g. for labeling a plot? Example: I want to plot several qqnorm-plots side by side and there should be a maintitle for each qqnorm-plot which is identical to the respective colname. I checked, but the column which is processed by the function called by apply does not contain a colname

...____________________________ > R-help@r-project.org mailing list > > PLEASE do read the posting guide > and provide commented, minimal, self-contained, reproducible code. ------------------------------ Message: 2 Date: Thu, 25 Dec 2008 19:51:36 +0100 From: Leo G?rtler <leog@anicca-vijja.de> Subject: [R] p(H0|data) for lm/lmer-objects R To: r-help@stat.math.ethz.ch Message-ID: <4953D638.4090507@anicca-vijja.de> Content-Type: text/plain; charset=ISO-8859-15 Dear R-List, I am interested in the Bayesian view on parameter estimation for multilevel models and ordinary r...

Dear R-list, I try to transform a mathematica script to R. #######relevant part of the Mathematica script (* p_sv *) dd = NN (DsD - DD^2); lownum = NN (L-DD)^2; upnum = NN (H-DD)^2; low = lownum/(2s^2); up = upnum/(2s^2); psv = NIntegrate[1/(s^NN) Exp[-dd/(2s^2)] (Gamma[1/2,0,up] + Gamma[1/2,0,low]),{s,sL,sH}, MinRecursion->3]; PSV = psv/Sqrt[2NN]; Print["------------- Results

Dear list members, how can multiple imputation realized for anova tables in R? Concretely, how to combine F-values and R^2, R^2_adjusted from multiple imputations in R? Of course, the point estimates can be averaged, but how to get standarderrors for F-values/R^2 etc. in R? For linear models, lm.mids() works well, but according to Rubins rules, standard errors have to be used together with

Dear alltogether, I tried pan() to impute NAs for longitudinal data. The terminology in the following output follows the pan manpage. No data are attached to this script as this may be too huge. y = 15 responses pred = at first just intercept was tried (later on covariates should follow) subj = 168 different subjects with 4 to 6 observations for each subject at time points t1, t2, ..., t6 #

Dear altogether, I want to blank the lower (or upper) part of a correlation matrix as it is done by dist() example: ( d <- cor(matrix(runif(12),nrow=4)) ) If I do the following d[lower.tri(d)] <- "" of course everything is changed to character - that's not what should be. Additionally, it does not work to assign "0" or anything else. The same is true for

Dear alltogether, I want to use CGIwithR in conjunction with R2HTML. A small example called 'test.R': ##### #! /usr/bin/R invisible(capture.output(library(R2HTML))) HTML(summary(as.numeric(scanText(formData$numbers))), file=stdout()) ##### The script gets its input via 'CGIwithR.cgi' and contains the variable "numbers." The 'HTML' output (-> summary() in

Dear R-list members, I have a problem with translating a mathematica script into R. The whole script is at the end of the email (with initial values for easy reproduction) and can be pasted directly into R. The problematic part (which is included below of course) is <--- Original Mathematica ---> (* p_svbar *) UiA = Ni (Dsi - 2Di A + A^2)/2; UiiA = Nii (Dsii - 2Dii A + A^2)/2; psvbar =