similar to: batch mode problem

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

Hola, amigos Se me acabo la sapienza. A ver: tengo un sript de Python script.py que neesita sus argumentos arg1 arg2 arg3, si desde la consola del sistema me pongo en la carpeta en la que lo tengo y hago: $ script.py arg1 arg arg3 Funciona perfecto. Pero desde R no lo consigo (entiendo que por no saber) Lo hago así (copio y pego): dir <- getwd() setwd('data/row/directoriocondatos/')

Hi sambalist, There are some questions I want to ask you. I hope that you could give some advice to me. I want to export all the groups membership to a file such as a csv for someone to audit. I know there are some command I can list group members like: list all the group: samba-tool group list and list members in some group: samba-tool group listmembers groupname so how can I merge all the

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

I have two hosts as Xen Server named host A and host B in the same network. A VM is created on Host A. Host A and B have a shared storage using NFS. While doing migration from host A to host B, how can I see the memory pages being migrated. If I run a small continuous program on the VM residing on host A, will the migration complete first or the sript running will complete first? How can I see

Ni idea de eso; te lo agradezco. Seguramente lo voy a terminar usando bastante. Pregunto: ahí debería hacer todas las operaciones propias del sistema? No me encuentra el script, y se me hace raro. Capaz que pasa al sistema desde la carpeta en la que está el script de R, no desde el wd de R. Ahora me toca salir y después intentar probar y que funcione. Gracis de vuelta Hau idatzi du Marcelino De

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

Hello, Situation: We are in a school class every student logs on with the same account. Until now we had Shares that were accessable for everybody and it was therefore possible to look into and edit/delete other's files. Plan: Create a share that can be clicked on which then asks for User/Pass and directly maps to the User's home Directory upon auth. User auth is done through LDAP

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

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

Hello Windows Users! I have finally got my new computer with it''s win32 build working and have just done the first pre-release for 0.3.13 that has win32 support. *** All win32 users should try to grab the pre-release and tell me if it works. *** You can install it with: $ gem install mongrel --source=http://mongrel.rubyforge.org/releases/ >From EITHER win32 or any other platform