[R] blank upper or lower triangle of cor-matrix

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

Re: [R] CGIwithR and visible output of 'invisible(capture.output(library(...)))'

Eric schrieb: library(R2HTML, verbose = FALSE) Dear Eric, thanks - I tried that, but that does not work for me. And it does not depend whether R2HTML is loaded via script/ batch job (- webserver, cgi) or directly within R. The same happens (of course) by using require(R2HTML, quietly=TRUE)

Re: [R] CGIwithR and visible output of 'invisible(capture.output(library(...)))'

platform, version of R, R2HTML or anything else. But almost certainly the problem is the command given by getOption(unzip), which looks like an unzip clone not respecting the -q flag. Try using the 'real' unzip () instead, from http://www.info-zip.org/. On Sat, 27 Jan 2007, Leo Gürtler

[R] integration problem with gamma function

, 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

[R] mathematica - r (gamma function + integration)

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] +

[R] effect sizes in lme/ multi-level models

in advance, leo gürtler -- email: [EMAIL PROTECTED] www: http://www.anicca-vijja.de/ __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Re: [R] norm package prelim.norm

(Ted Harding) schrieb: Dear Elizabeth, the R-package mice has ?md.pattern and Hmisc also contains some software to describe missing data. I think both do not have limitations like prelim.norm. best wishes, leo On 01-Feb-06 Ted Harding wrote: On 01-Feb-06 Elizabeth Lawson wrote:

[R] lmer and handling heteroscedasticity

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

[R] how to obtain par(ask=TRUE) with trellis-plots

Dear alltogether, how can a delay like possible with par(ask=TRUE) be attained while using trellis-plots within a loop or something like that? the following 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 -

[R] decide between polynomial vs ordered factor model (lme)

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

[R] convergence error (lme) which depends on the version of nlme (?)

of the convergence error which seems to depend on the built of R/nlme? Thank you very much. Best wishes, leo gürtler -- email: [EMAIL PROTECTED] www: http://www.anicca-vijja.de/ __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman

[R] strange behavior of loess() predict()

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 __ R-help@stat.math.ethz.ch mailing list https

Re: [R] strange behavior of loess() predict()

- that predict.loess predicts values that are almost 2/max(y) ~ 5000 times higher than expected? best, leo gürtler On Tue, 2005-12-06 at 18:09 +0100, Leo Gürtler wrote: Dear altogether, snip # here is the difference!! predict(mod, data.frame(x=X), se=TRUE) predict(mod, x=X, se=TRUE) --- end

[R] analzying multiple variables (dv) in a sequence by using fit.mult.impute together with a MICE.object

not work, because it is just the name of the list element. Otherwise by specifying one of the columns with the list element, the multiple datasets would be missed. I thought about using eval() or expr() but without real succes. Every hint is appreciated, best regards leo gürtler

[R] multiple imputation of anova tables

, mice/pan/mix or Hmisc. Any help and hints are appreciated, best regards Leo Gürtler / Germany __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting

[R] error in NORM lib

. 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

[R] plot of mice.mids objects

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 seems that there are too many variables.

[R] problems with pan(): Indizierung ausserhalb der Grenzen = subscript out of bounds

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