[R] MySQL mismatch

2004-07-21 Thread Joerg Schaber
Hi, library(RMySQL) works, but connecting to the database results in Warning message: RS-DBI driver warning: (MySQL mismatch between compiled version 3.23.49 and runtime version 3.23.51) Any idea? Thanks, joerg __ [EMAIL PROTECTED] mailing list

[R] trend turning points

2004-04-14 Thread Joerg Schaber
Hi, does anybody know of a nice test to detect trend turning points in time series? Possibly with reference? Thanks, joerg __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

[R] installing RMySQL

2004-04-07 Thread Joerg Schaber
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

[R] combined random effects

2004-03-24 Thread Joerg Schaber
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(),

[R] wolfinger microarray normalization

2004-03-24 Thread Joerg Schaber
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

[R] F Dist

2004-02-19 Thread Joerg Schaber
Hi, why is it that in tables for the F distribution it is for example F0.95[6,28]=2.45 or F0.95[10,20]=2.35 but pf(6,28,2.45) [1] 0.8854934 pf(10,20,2.35) [1] 0.9300167 Thanks for clarifying, jeorg __ [EMAIL PROTECTED] mailing list

[R] Comparing two regression slopes

2004-02-19 Thread Joerg Schaber
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

[R] levelplot parameters

2003-12-09 Thread Joerg Schaber
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 __ [EMAIL PROTECTED] mailing list

[R] package compilation

2003-11-18 Thread Joerg Schaber
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

[R] \samp tag

2003-11-18 Thread Joerg Schaber
Hi including the \samp tag in Rd files gives me errors creating the tex-documentation. How can I escape literal sequences in Rd files? \preformatted also does not seem to work. joerg __ [EMAIL PROTECTED] mailing list

[R] R CMD check problems

2003-11-14 Thread Joerg Schaber
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

[R] conf int mixed effects

2003-11-13 Thread Joerg 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

Re: [R] conf int mixed effects

2003-11-13 Thread Joerg Schaber
, the computed CI will cover the true population constant. There's no true population constant for random effects, but there is for a variance component. HTH, Andy -Original Message- From: Joerg Schaber [mailto:[EMAIL PROTECTED] Sent: Thursday, November 13, 2003 10:50 AM To: Douglas Bates

[R] test for connectedness

2003-10-29 Thread Joerg Schaber
Hi, in a 2-way linear model with missing data the normal equations can sometimes be separated in 2 or more independent sets of equations brought about by the nature of the missing data. As far as I know this phenomenon is called 'connectedness' of the data. Does anyone know of a test or

[R] Oracle fetch problems

2003-10-27 Thread Joerg Schaber
relating to my former messages concerning the strange fetch problems I have, I found out that 'fetch' fetches only every second row of the 'native' SQL-results and fills up the rest of the rows with zeros. Here a draft for the native SQL-query: 12:43:07 SQL SELECT OBS_DAY, OBS_YEAR, STAT_ID

[R] Oracle fetch problems

2003-10-24 Thread Joerg Schaber
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 =

[R] Oracle fetch problems

2003-10-24 Thread Joerg Schaber
Concerning my former message: I just found out that irrespective of the fetch size, I get always half of the fetch size good results/rows and the other half is corrupted. Strange, isn't? __ [EMAIL PROTECTED] mailing list

[R] 64-bit Oracle interface

2003-09-24 Thread Joerg Schaber
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

Re: [R] robust regression

2003-08-10 Thread Joerg Schaber
User) Joerg Schaber wrote: Hi, trying to do a robust 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

[R] robust regression

2003-07-30 Thread Joerg Schaber
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

[R] linear model coefficients

2003-07-28 Thread Joerg Schaber
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