search for: rtler

...ructure is present, the formular 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/

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

...y, ylim=range(y, lower, upper)) lines(X, fit$fit) # strange values in fit fit # here is the difference!! predict(mod, data.frame(x=X), se=TRUE) predict(mod, x=X, se=TRUE) <--- end of snip ---> I assume this has some reason but I do not understand this reason. Merci, best regards leo g??rtler

...don't understand this, because the dataset works well with MICE and aregImpute() 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...

The vp3 / theora decoder in FFmpeg project received a major speedup: http://multimedia.cx/eggs/profiling-and-optimizing-theora/#comment-150751 http://git.ffmpeg.org/?p=ffmpeg;a=commit;h=6e8061980c9a18a0eb9802f5dc5becc349066680 http://git.ffmpeg.org/?p=ffmpeg;a=commit;h=8476744ba50d26a3e188def192cd9086b0f781d7 http://git.ffmpeg.org/?p=ffmpeg;a=commit;h=a51faef9b61bd80499a28322802ef4b22be44d43

...are slightly different. Thus my question: 1) Do I have overseen a switch for the structure of the random effects? Is something wrong with the call/ formular? 2) What is the cause of the convergence error which seems to depend on the built of R/nlme? Thank you very much. Best wishes, leo g??rtler -- email: leog at anicca-vijja.de www: http://www.anicca-vijja.de/

...ing the list at r-help-owner@r-project.org When replying, please edit your Subject line so it is more specific than "Re: Contents of R-help digest..." Today's Topics: 1. Re: Implementing a linear restriction in lm() (Ravi Varadhan) 2. p(H0|data) for lm/lmer-objects R (Leo G?rtler) 3. Re: beginner data.frame question (Oliver Bandel) 4. Re: 4 questions regarding hypothesis testing, survey package, ts on samples, plotting (Thomas Lumley) 5. Re: How can I avoid nested 'for' loops or quicken the process? (Oliver Bandel) 6. Re: 4 questions regardin...

..."sqrt(pi)" because of (-> manpage pgamma) differences in formulating the incomplete gamma function between Mathematica and R. However, I tried this in various ways and - at least for me - it does not help. Thus, I assume it is another issue. Thank you very much, best wishes leo g?rtler now the R script to reproduce (can be pasted directly into R): #################### # R-Portierung aus Mathematica (Urban Studer, 90er) # Ursprung: G.L. Bretthorst "On the difference of means" # zuerst: 12-06-05 # zuletzt: 21-06-06 # -------------------------------------------------- #...

...F denDF F-value p-value (Intercept) 1 80 4123.156 <.0001 age 1 80 114.838 <.0001 Sex 1 25 9.292 0.0054 I am confused about that and I did not find any hint in norm, mice/pan/mix or Hmisc. Any help and hints are appreciated, best regards Leo G??rtler / Germany

...ject=Re:%20%5bR%5d%20plot%20of%20mice.mids%20objects> > Date: Fri, 04 Nov 2005 12:57:05 +0100 Hi, it is impossible to plot a subset of the variables at the moment, but we could add it. p.s. The latest version of mice is v1.14, which can be downloaded from CRAN. Roel de Jong. Leo G?rtler wrote: > Hello... > > how can I plot mice.mids objects as described by Buuren (2000) > http://web.inter.nl.net/users/S.van.Buuren/mi/docs/Manual.pdf > page 17, if there are many variables (~80) with NAs included? > mice runs well, but the plot is not possible because it seem...

Dear altogether, is it possible to integrate "weights" arguments within lmer to incorporate statements to handle heteroscedasticity as it is possible with lme? I searched the R-archive but found nothing, insofer I assume it is not possible, but as lmer is under heavy develpoment, maybe something changed or is solved differently. Thus my question: While encountering heavy

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

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 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 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) within smallgroups