search for: schaber

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

...+ s - 1, method = "M", contrasts = list(s = ("contr.sum"))) I tried an exact sampling (psamp="exact"), but I keep getting syntax errors. Any idea how I can make the first one work? Thanks, joerg -- ---------------------------------------------------------- J?rg Schaber Instituto Cavanilles de Biodiversidad y Biologia Evolutiva Universidad de Valencia Tel.: ++34 96 354 3666 A.C. 22085 Fax.: ++34 96 354 3670 46071 Valencia, Espa?a email : jos at uv.es

...but behind that there is the user's homedir or a share that lists /home and asks for a User/Pass for each dir you click on. I know this is party done by setting appropriate rights on the home dirs (700). I hope I made everything clear :) Thanks alot for your ideas! -- Greetings Alexander Schaber http://www.alexanderschaber.de/

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

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

...;http://cran.r-project.org/bin/windows/contrib >Uwe Ligges As for Feeney compiling from source is not a preferable option for me, as I have no experience on that. In the link provided by Uwe Ligges, I could not locate the file I am looking for (xml.** zip). Thanks a lot for your help, Katrin Schaber

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",

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, 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, 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 <-

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