[R] RdbiPgSQL and POSIXct

2003-12-10 Thread Miha STAUT
Hi, Does the Bioconductor library RdbiPgSQL support the POSIXlt and POSIXct (R time format)?. I tried to import an R data frame which had also some columns of class POSIXct into PostgreSQL 7.2.1 with the command dbWriteTable without success. str(jame) `data.frame': 8123 obs. of 27 variables:

Re: [R] pvalues

2003-12-10 Thread Uwe Ligges
Maya Sanders wrote: dear all- If I have a vector of numbers (not necessarily normally distributed) how can I get the p-value of a number in this distribution. I am interested in the inverse of 'quantile' . thank you- Maya I'm not sure whether I got the point, but my guess is that you are

Re: [R] problem with pls(x, y, ..., ncomp = 16): Error in inherit s( x, data.frame) : subscript out of bounds

2003-12-10 Thread Ron Wehrens
On Wednesday 10 December 2003 05:14, Liaw, Andy wrote: My first guess was wrong (should have RTFM...). Try: m - pls(x, y, validation = CV, niter = 68, ncomp = 1:16) I'm not familiar enough with SIMPLS to know whether it makes sense to pass ncomp of length 1 to pls(), but the code as is

[r] interfacing r and python in ms windows

2003-12-10 Thread Hogendoorn, Rene
I have a windows rpy for python 2.3 if you like. Regards, René Hogendoorn __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

[R] How to calculate standard error for a vector?

2003-12-10 Thread Laurent Houdusse
Hi all! I 'm beginner and i develop a bio-application with VB and i need some statistic functions! could i calculate StdError, CoeffOfVariance, SumSquared with R langage? if yes, what are functions to use? I need also to use ANOVA and t-test... Thanks for your help! Laurent Houdusse

Re: [R] How to calculate standard error for a vector?

2003-12-10 Thread Spencer Graves
Have you read An Introduction to R and The R Language Definition that comes with the R installation? Try help.start(). Also, there are several good introductions downloadable from www.r-project.org. I suggest you also consider lm = linear model, e.g., fit - lm(y~x, dataFrame) #

[R] Error in model.frame

2003-12-10 Thread Göran Broström
I recently ran into the following: x - c(1,2,4,3) lm(x ~ 1:length(x)) Error in model.frame(formula, rownames, variables, varnames, extras, extranames, : variable lengths differ but lm(x ~ c(1:length(x))) Call: lm(formula = x ~ c(1:length(x))) Coefficients: (Intercept)

RE: [R] pvalues

2003-12-10 Thread Adaikalavan RAMASAMY
If you want to generate a p-value, what is the null hypothesis ? -- Adaikalavan Ramasamy -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Uwe Ligges Sent: Wednesday, December 10, 2003 3:43 PM To: Maya Sanders Cc: [EMAIL PROTECTED] Subject: Re: [R]

Re: [R] Error in model.frame

2003-12-10 Thread Martyn Plummer
On Wed, 2003-12-10 at 10:34, Göran Broström wrote: I recently ran into the following: x - c(1,2,4,3) lm(x ~ 1:length(x)) Error in model.frame(formula, rownames, variables, varnames, extras, extranames, : variable lengths differ but lm(x ~ c(1:length(x))) Call:

RE: [R] How to calculate standard error for a vector?

2003-12-10 Thread Adaikalavan RAMASAMY
sd - sqrt(var(x)) # standard deviation cv - sd / mean(x) # coefficient of variation ss - mean( x^2 ) # sum of squares - definitions vary ! If these are the statistics you are going to be using, it might be more efficient to code it directly. -- Adaikalavan Ramasamy -Original Message-

Re: [R] Error in model.frame

2003-12-10 Thread Uwe Ligges
Göran Broström wrote: I recently ran into the following: x - c(1,2,4,3) lm(x ~ 1:length(x)) Error in model.frame(formula, rownames, variables, varnames, extras, extranames, : variable lengths differ but lm(x ~ c(1:length(x))) Call: lm(formula = x ~ c(1:length(x))) Coefficients:

RE : [R] How to calculate standard error for a vector?

2003-12-10 Thread Laurent Houdusse
Thanks for this formulae For stdError, i found this: stdError - sd(x)/sqrt(length(x)) Laurent Houdusse Analyste Programmeur -Message d'origine- De : Adaikalavan RAMASAMY [mailto:[EMAIL PROTECTED] Envoyé : mercredi 10 décembre 2003 10:52 À : Laurent Houdusse; [EMAIL PROTECTED] Objet

Re: [R] Error in model.frame

2003-12-10 Thread Göran Broström
On 0, Uwe Ligges [EMAIL PROTECTED] wrote: Göran Broström wrote: I recently ran into the following: x - c(1,2,4,3) lm(x ~ 1:length(x)) Error in model.frame(formula, rownames, variables, varnames, extras, extranames, : variable lengths differ but lm(x ~ c(1:length(x)))

[R] NLME and limits on parameter space

2003-12-10 Thread Karsten D Bjerre
Dear R-users, How can I impose limits on parameter estimates in NLME as indicated in the second #'ed line below? Can I use the naPattern, some thing like this: nlme(response = yftp(x,ii,cc,b,dd),... , naPattern= is.na(response) ) ? yftp-function(x,ii,cc,b,d){ tpx-tpx(ii,cc,b)# x

[R] How to remove extra spaces introduced (?) by write.table

2003-12-10 Thread Valtteri Wirta
Dear Group, My question relates to the write.table function. I have created a data.frame containing both character and numeric columns. When I use the write.table function to write this table into a text file spaces are added into the numeric columns so that the number of digits in each column

[R] Re: p-value from chisq.test working strangely on 1.8.1

2003-12-10 Thread Federico Calboli
Mostly for the record, it happens also on Mandrake 9.1 with R 1.8.1 installed from rpm. chisq.test(matrix(c(0, 1, 1, 12556), 2, 2), simulate.p.value=TRUE) Pearson's Chi-squared test with simulated p-value (based on 2000 replicates) data: matrix(c(0, 1, 1, 12556), 2, 2)

Re: [R] How to remove extra spaces introduced (?) by write.table

2003-12-10 Thread Prof Brian Ripley
Just convert to character before using write.table: as.character(c(1, 10)) [1] 1 10 whereas format(c(1,10)) [1] 1 10 is what you are seeing. On Wed, 10 Dec 2003, Valtteri Wirta wrote: Dear Group, My question relates to the write.table function. I have created a data.frame

Re: [R] setMethod(min, myclass, ...)

2003-12-10 Thread Peter Ruckdeschel
Hello, first of all thank you for your reply to our help request. On Wed, 3 Dec 2003, John Chambers wrote: Thomas Stabla wrote: I have defined a new class setClass(myclass, representation(min = numeric, max = numeric)) and want to write accessor functions, so that for foo =

[R] Scatterplot axes

2003-12-10 Thread Ellen Kalmbach
Please, could someone help me figure out what seems to be a very simple problem (and is still taking me hours...). I want to draw a simple scatterplot but with 'equal' axes, i.e. I want both axes to go from -3 to 3. Values for x lie between -2 and 0.5, values for y between -2.2 and 3. I have tried

RE: [R] Scatterplot axes

2003-12-10 Thread Pfaff, Bernhard
Please, could someone help me figure out what seems to be a very simple problem (and is still taking me hours...). I want to draw a simple scatterplot but with 'equal' axes, i.e. I want both axes to go from -3 to 3. Values for x lie between -2 and 0.5, values for y Hello Ellen, how

Re: [R] Scatterplot axes

2003-12-10 Thread Giovanni Petris
Have you tried specifying xlim and ylim? See ?par. -- __ [ ] [ Giovanni Petris [EMAIL PROTECTED] ] [ Department of Mathematical Sciences ] [ University of Arkansas -

Re: [R] Scatterplot axes

2003-12-10 Thread Thomas W Blackwell
Ellen - plot(my.x.vector, my.y.vector, xlim=c(-3,3), ylim=c(-3,3)) It's the named arguments xlim and ylim that you were looking for. I frequently set them as xlim = 3 * c(-1,1), ylim = 3 * c(-1,1) so that I can change the range by editing just one number rather than two. For an added fillip,

[R] custom kernel

2003-12-10 Thread Miha STAUT
Hi, I would like to estimate a 2D density on a matrix but not with a predefined kernel shape and density. I am sure there is a way to do this only I did not find it. Any help is kindly appreciated. Thanks, Miha Staut __ [EMAIL PROTECTED] mailing

Re: [R] PROC MIXED vs. lme()

2003-12-10 Thread Christophe Pallier
Hi I'm trying to learn how to do a repeated measures ANOVA in R using lme(). In SAS, I use the code: PROC MIXED DATA=[data set below]; CLASS sub group trial; MODEL dv = group trial group*trial; REPEATED trial / SUBJECT=sub TYPE=CS; In R, I'm trying the code: results.cs - lme(DV ~

[R] factorial experiments and repeated measures

2003-12-10 Thread eesteves
Dear All, I'm new to the list and relatively new to R (but not so to S-Plus) and I´ve been asked to help study some experimental data. In a factorial experiment (5 factors at 2 levels) the response-variable was measured at various times (0, 2, 5, 7, ..., 30 days after start). There's just one

Re: [R] OT: BibTex year-only citation in text?

2003-12-10 Thread Andrew Perrin
Generally, use natbib style and then \citep- (Anderson 1992) \citet- Anderson (1992) \citeyearpar - (1992) \citealt - Anderson 1992 -- Andrew J Perrin - http://www.unc.edu/~aperrin Assistant Professor of Sociology, U of

RE: [R] Windows Memory Issues

2003-12-10 Thread JFRI (Jesper Frickmann)
I recommend you get the latest version 1.8.1 beta. I also had some memory problems and that fixed it. Kind regards, Jesper Frickmann Statistician, Quality Control Novozymes North America Inc. Tel. +1 919 494 3266 Fax +1 919 494 3460 __ [EMAIL

Re: [R] Windows Memory Issues

2003-12-10 Thread Duncan Murdoch
On Wed, 10 Dec 2003 10:43:12 -0500, JFRI (Jesper Frickmann) [EMAIL PROTECTED] wrote : I recommend you get the latest version 1.8.1 beta. I also had some memory problems and that fixed it. 1.8.1 has been released, so there's no more beta. The release is available on CRAN (see

Re: [R] Font problem

2003-12-10 Thread John Dougherty
Thank you for the help. The problem was, as Peter and Marc suggested, that the 100 dpi fonts were missing. Installing them made the difference. John __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

Re: [R] OT: BibTex year-only citation in text?

2003-12-10 Thread Torsten Hothorn
Jason, For many bibliography styles, the command \citeyear{key} will work. If this doesn't work for the style you are using, you can investigate style-specific methods or consider other styles. I find that natbib is good for author-year formats. Yup. \usepackage[round]{natbib} is

[R] configure stuck in checking leap seconds (R-1.8.1 on AIX)

2003-12-10 Thread Liaw, Andy
Dear R-help, We've been trying, so far unsuccessfully, to compile R as 64-bit under AIX 5.1. Following the recent post by Dr. Christoph Pospiech, I started with R-1.8.1 and manually edited the configure script according to the .diff file. Part of the diff has: *** *** 24446,24452

[R] an eval/parse trivia

2003-12-10 Thread Vincent Detours
Dear all, Any cue on how to evaluate x? x [1] structure(c(GO:0004707, GO:0005524, GO:0004674, GO:0006468, GO:074, GO:0008372, GO:0016740), .Names = c(NAS, NAS, IEA, IDA, NAS, ND, IEA)) 6204 Levels: GO:074 GO:158 GO:163 GO:166 ... structure(c(GO:0019538, GO:0016706), .Names =

Re: [R] an eval/parse trivia

2003-12-10 Thread Thomas W Blackwell
Vincent - From the values shown, this looks like a Bioconductor question, rather than base R. You might try the maintainers of whatever package the function comes from. Is 2287 the index in levels for one of the character strings shown ? - tom blackwell - u michigan medical school - ann

[R] (no subject)

2003-12-10 Thread Dan Olsen
What commands can I use to find a confidence interval for a coefficient for a quadratic regression test? Dan Olsen Brigham Young University Continuing Education Computer Operations 191 HCEB (801)378-3232 [EMAIL PROTECTED]

Re: [R] an eval/parse trivia

2003-12-10 Thread rpugh
Is the data structure a factor variable containing character strings? That would account for the lack of quotes and the Levels bit. If so, try eval(parse(text=as.character(...))). x - factor('structure(1:3, .Names=LETTERS[1:3])') eval(parse(text=as.character(x))) A B C 1 2 3 Weird data

[R] extensive grid docs?

2003-12-10 Thread Jeff D. Hamann
I'm looking for extensive docs on using grid (for the somewhat newbie). I'm attempting to create a set of graphics that look similar to the attached image (I hope this doesn't get bounced) and have only come across the R newsletters and it appears that grid was new as of 1.8.0? I think the best

[R] RODBC with RAqua

2003-12-10 Thread Peter McMahan
Hi, I use RAqua on mac os x 10.3. I am trying to connect to an odbc dsn, but RODBC can't seem to find it. I have it installed in my ODBC Administrator in /Applications/Utilities as a system-wide data source. does RODBC need more than just the dsn name (like a path or something)? I have

[R] How to start RMySQL

2003-12-10 Thread Jeff Lee
I am a newbie in R and want to use MySQL database. Here is what I have and done. R 1.8.0 mysql 4.0.16 Both running in Windows XP. I download the RMySQL.zip from the internet and used the Packages installer in the RGui to install RMySQL. I also installed DBI packages. I ran the following

[R] dyn.load for c code

2003-12-10 Thread K Okamoto
I am learning how to load C code into R-1.8.0 on Windows 98. To this end I wrote a small c program, downloaded the tools, perl, and mingw on the building R for windows page, and proceeded to create libR.a libRblas.a as explained in the readme.packages. I started with a simple c program called

[R] How can I call R from Java or C?

2003-12-10 Thread Aylin Yilmaz
Hi all, I am new to R. I have been looking for a way to call R functions from Java or C. What I have found so far is: I can execute R commands in batch mode from a file and redirect the output to another file. However, I am not sure this is the way to do it. About C, is it true that I can only

Re: [R] How to start RMySQL

2003-12-10 Thread David James
Jeff Lee wrote: I am a newbie in R and want to use MySQL database. Here is what I have and done. R 1.8.0 mysql 4.0.16 Both running in Windows XP. I download the RMySQL.zip from the internet and used the Packages installer in the RGui to install RMySQL. I also installed DBI

Re: [R] RODBC with RAqua

2003-12-10 Thread Prof Brian Ripley
On Wed, 10 Dec 2003, Peter McMahan wrote: Hi, I use RAqua on mac os x 10.3. I am trying to connect to an odbc dsn, but RODBC can't seem to find it. I have it installed in my ODBC Administrator in /Applications/Utilities as a system-wide data source. does RODBC need more than just the dsn