similar to: Making Home Directory available for Windows Users

Hi, I have a linear mixed-effects model object and want to extract the 95% confidence intervals for the fixed and random effects, respectively. I found the function intervals() for confidence intervals for the fixed effects but no corresponding function for the random effects. Does it exist or do I have to calculate the confidence intervals for the random effects myself? Greetings, joerg

Hi, does anybody know of a nice test to detect trend turning points in time series? Possibly with reference? Thanks, joerg

Hi, concerning my earlier mail, maybe someone has noted from the variable names that I try to analyse mircoarrary experiments. Does anybody know of a R-implementation of the two-step mixed-model normalization procedure proposed by Wolfinger et al. (2001) J. Comput. Biol. 8:625-637? That would be great, best, joerg

Hi, I have two non-normal distributions and use interquartile ranges as a dispersion measure. Now I am looking for a test, which tests whether the interquartile ranges from the two distributions are significantly different. Any idea? Thanks, joerg [[alternative HTML version deleted]]

Hi, trying to do a robudt regression of a two-way linear model, I keep getting the following error: > lqs(obs ~ y + s -1,method="lms", contrasts=list(s=("contr.sum"))) Error: lqs failed: all the samples were singular Robust regression with M-estimators works (also regular least square fits, of course): rlm.formula(formula = obs ~ y + s - 1, method = "M",

Hi, I have the following linear mixed model: y(g,i,j,k,l)=u + L(g) + T(i) + D(j) + S(k) + (TS)(i,k) + error(g,i,j,k,l) where S(k) and the combined effect (TS)(i,k) are random effects whereas the rest are fixed effects. How do I specifiy the random part of the model formula in lme(), especially concerning the combined effect (TS)? Moreover, when I run the model as a fixed effect model I get

Hi, I am trying to create a R-package. When I load my R source files using source() for debugging reasons everthing works fine and all function also do want they are supposed to do. However, when I run 'R CMD check' it stop with * checking R files for syntax errors ... ERROR Syntax error in file Is there a way to further investigate where the error is? greetings, joerg

Hello, I am using maptools for plooting geographical data. The colour of the region indicates some region dependent value (population for example). I pass the colours of the regions to the plot.Map function by defining the foreground colour: jet.colors = colorRampPalette(c("#00007F", "blue", "#007FFF", "cyan", "#7FFF7F", "yellow",

Hello, I am looking for the XML Package for Windows or any package, which would enable me to use the function "xmlTreeParse". In a thread dating from 2001, the same question is raised: >>Griffith Feeney wrote: >> >> Has anyone gotten the XML package (R News 1-1:24) with R running on Windows >> (95 or 98 or 2000)? I've looked at the omegahat, expat and libxml

Hi, I have problems installing RMySQL under R1.8.1 and RHlinux (kernel2.4.7). During compilation I get: ** libs gcc -I/usr/lib/R/include -I/usr/include/mysql -I/usr/local/include -D__NO_MATH_INLINES -mieee-fp -fPIC -O2 -m486 -fno-strength-reduce -g -c RS-DBI.c -o RS-DBI.o gcc -I/usr/lib/R/include -I/usr/include/mysql -I/usr/local/include -D__NO_MATH_INLINES -mieee-fp -fPIC -O2 -m486

Hi, I wonder if there is a possibility to avoid that R sets one level of a factor equal zero in a model fit. More precisely, I want to fit a two-way unbalanced linear model: o ~ 0 + x + y x is a factor with 10 levels, y is a factor with 9 levels. In order to get a unique solution i set the intercept =0 and impose that sum(y)=0 i.e. res <- lm(o ~ 0 + x + y,

Hi, I am running the 64bit version of R on a AIX RS6000 cluster on which is also running the ' SQL*Plus: Release 9.2.0.2.0' in the 'ORACLE-environment Version 64 bit AIX 5.1'. I have problems setting up an Oralce database connection. I can instantiate a driver using drv <- Oracle(). But when I try to connect to my database using dbConnect (dbConnect(drv,

Hi, I recently wrote a new R package, which I could 'R CMD check' and 'R CMD build' nicely on my Intel Linux PC. Loading the new library as a normal user on my PC also works nicely. This packages includes some compiled C-code in a shared library. However, before submitting the package to CRAN I wanted to install it also under an other platform in order to check if everything

Hi, I have a levelplot with one panel. I just can't find out how I can manipulate the size of the axis lables. e.g. scales.cex doesn't work, the usual par-parameters either. Any hint? joerg

Hi, ok, I know I should be using a later version than 1.7.1 (64 bit) but it's not in my power. So here is the problem: In my R script I declare a data.frame that consists of 40 vectors, each having 125 numeric elements. This is no problem as long as I run the sript in interactive mode, but running it in batch mode I get strange error messages. Apparently, it has to do with the size of

I would suggest the method of Sokal and Rholf (1995) S. 498, using the F test. Below I repeat the analysis by Spencer Graves: Spencer: > df1 <- data.frame(x=1:3, y=1:3+rnorm(3)) > df2 <- data.frame(x=1:3, y=1:3+rnorm(3)) > fit1 <- lm(y~x, df1) > s1 <- summary(fit1)$coefficients > fit2 <- lm(y~x, df2) > s2 <- summary(fit2)$coefficients > db <-

Hokay. It appears that enumdrivers / getdriver in rpcclient is broken, and won't be fixed soon. I'd fix it meself, but hey, I'm too busy clearing skyscrapers in a single bound and moving faster than a speeding bullet. Or maybe I looked at the cvs once and got REALLY scared. Unfortunately, I can't wait, so I've come up with a nasty hack to help us with querying samba for the

I thought I had Dovecot configured to use Maildir instead of mbox, but apparently it is still using mbox. Excerpt from dovecot.conf: mail_location = maildir:~/Maildir Excerpt from sendmail.mc: MAILER(dovecot)dnl dovecot.m4: ######################################## ### DOVECOT Mailer Specification ### ######################################## Mdovecot,

Hi, I use the ROracle interface in the R 64bit version. When I fetch a large number of rows from the database I get alternating 250 rows of what I want to see and than 250 rows junk and so on. It does not matter whether I use the fetch command or dbGetquery or whatever. My fetch_default_rec = 500 and I can not change it by setting it explicitly in Oracle(). Any idea what's the problem

Hi, I fitted a linear model: fit <- lm(y ~ a * b + c - 1 , na.action='na.omit') Now I want to extract only the a * b effects with confidence intervals. Of course, I can just add the coefficients by hand, but I think there should an easier way. I tried with predict.lm using the 'terms' argument, but I didn't manage to do it. Any hints are appreciated, best, joerg