[R] Solving a maximization problem using QUADPROD

2006-11-08 Thread Serguei Kaniovski
Hello, here is an example from the manual. How to turn this minimization problem into maximization problem, i.e. -(0 5 0) %*% b - 1/2 b^T b? # Assume we want to minimize: -(0 5 0) %*% b + 1/2 b^T b # under the constraints: A^T b = b0 # with b0 = (-8,2,0)^T # and (-4 2 0) # A = (-3 1 -2) # ( 0 0

[R] Cost function in cv.glm for count data

2006-11-08 Thread [EMAIL PROTECTED]
Hi everybody, I'm trying to use cross-validation (cv.glm) for count data. Does someone know which is the appropriate cost function for Poisson distribution? Thank you in advance. Valerio. -- Scopri il nuovo MOTOROLA K1, il cellulare con

Re: [R] [Rd] How to link to vignette from .Rd file

2006-11-08 Thread Antonio, Fabio Di Narzo
2006/11/8, Agner Fog [EMAIL PROTECTED]: Thanks for the tip. I had thought of using the demo facility for this but I didn't know how to check for the operating system. See help about '.Platform'. I.e., .Platform$OS.type can be useful. This is a clumsy solution, but better than nothing. I

Re: [R] CPU or memory

2006-11-08 Thread Stefan Grosse
64bit does not make anything faster. It is only of use if you want to use more then 4 GB of RAM of if you need a higher precision of your variables The dual core question: dual core is faster if programs are able to use that. What is sure that R cannot make (until now) use of the two cores if you

Re: [R] gamm(): nested tensor product smooths

2006-11-08 Thread Simon Wood
I'd like to compare tests based on the mixed model representation of additive models, testing among others y=f(x1)+f(x2) vs y=f(x1)+f(x2)+f(x1,x2) (testing for additivity) In mixed model representation, where X represents the unpenalized part of the spline functions and Z the wiggly

[R] get compressed data via a socket connection

2006-11-08 Thread Simon Penel
Dear R developers I am currently working on the seqinR package. The seqinR package allows a remote access to biological databases via a socket connection. We are using the functions socketConnection, writeLines and readLines to open the socket, send request to the server and receive

Re: [R] CPU or memory

2006-11-08 Thread Prof Brian Ripley
On Wed, 8 Nov 2006, Stefan Grosse wrote: 64bit does not make anything faster. It is only of use if you want to use more then 4 GB of RAM of if you need a higher precision of your variables The dual core question: dual core is faster if programs are able to use that. What is sure that R

Re: [R] snow's makeCluster hanging (using Rmpi)

2006-11-08 Thread Ramon Diaz-Uriarte
On Tuesday 07 November 2006 19:28, Randall C Johnson [Contr.] wrote: On 11/7/06 11:28 AM, Ramon Diaz-Uriarte [EMAIL PROTECTED] wrote: On Tuesday 07 November 2006 15:56, Randall C Johnson [Contr.] wrote: Hello everyone, I've been fiddling around with the snow and Rmpi packages on my new

[R] Convert ordinary dates into POSIX

2006-11-08 Thread Ivan Kalafatic
Does anyone know where I can find any tool for Windows that converts dates from ordinary formats into POSIX? I need it to import sime time series from Excel into R and use them with its package. Thanks! __ R-help@stat.math.ethz.ch mailing list

[R] Fw: data frames re-ordering and naming columns

2006-11-08 Thread Jabez Wilson
Thank you for help with the first part. The second part I do not understand. All I get is 3 3. Presumably this is the factor that you refer to. Is there no simple way of extracting the text alone - Forwarded Message From: Peter Alspach [EMAIL PROTECTED] To: Jabez Wilson [EMAIL

[R] axis command and excel time format

2006-11-08 Thread Carmen Meier
Hi to all, I have some problems to get the times-scale to the x-axis the times are coming from an excel sheet f. e [1] 0:01:00 0:02:00 0:03:00 0:04:00 0:05:00 0:06:00 0:07:00 [8] 0:08:00 0:09:00 0:10:00 0:11:00 0:12:00 0:13:00 0:14:00 [15] 0:15:00 0:16:00 0:17:00 0:18:00 0:19:00 0:20:00 0:21:00

Re: [R] Convert ordinary dates into POSIX

2006-11-08 Thread Prof Brian Ripley
?strptime in R. However, if you want to convert dates (and not date-times), see ?as.Date, using the 'format' argument. On Wed, 8 Nov 2006, Ivan Kalafatic wrote: Does anyone know where I can find any tool for Windows that converts dates from ordinary formats into POSIX? I need it to import

Re: [R] query in R

2006-11-08 Thread Mark Wardle
Xiaodong Jin wrote: how to realize the following SQL command in R? select distinct A, B, count(C) from TABLE group by A, B ; quit; Best Regards You don't say which database you are using, but I have found RODBC to be very effective. library(RODBC) Have a look at

Re: [R] axis command and excel time format

2006-11-08 Thread Gabor Grothendieck
Try this: plot(z, xaxt = n) xt - paste(23, seq(5, 50, 5), sep = :) axis(1, times(paste(xt, 0, sep = :)), xt) On 11/8/06, Carmen Meier [EMAIL PROTECTED] wrote: Hi to all, I have some problems to get the times-scale to the x-axis the times are coming from an excel sheet f. e [1] 0:01:00

Re: [R] Convert ordinary dates into POSIX

2006-11-08 Thread Gabor Grothendieck
There is a discussion and code in the R News 4/1 help desk article. On 11/8/06, Ivan Kalafatic [EMAIL PROTECTED] wrote: Does anyone know where I can find any tool for Windows that converts dates from ordinary formats into POSIX? I need it to import sime time series from Excel into R and use

Re: [R] changing image dimensions

2006-11-08 Thread Oleg Sklyar
'EBImage' R-package from www.bioconductor.org. It's based on ImageMagick and will allow you to do this and much more with images of any time. The prerequisite: either Linix/Unix or Mac (from source) -- will not work on Windows. On 07/11/06, Milton Cezar Ribeiro [EMAIL PROTECTED] wrote: Hi

Re: [R] axis command and excel time format

2006-11-08 Thread Carmen Meier
Thank you for your reply Gabor, sure, the manually written axis works fine in any configuration. but I would prefer an automatic input. That means that I would like to use the datafield[1] for the minimum time and the datafield[max] (means the last one) for the maximum time. divided into x

Re: [R] axis command and excel time format

2006-11-08 Thread Gabor Grothendieck
Is the problem how to produce an axis with a given minimum tick, maximum tick and given number of ticks? In that case try this (if not explain further): # input data # z is from original example mn - times(23:00:00) mx - times(23:55:00) n - 12 xt - times(seq(mn, mx, length = n)) plot(z, xaxt =

[R] Building Problem

2006-11-08 Thread JSteiner
Hi, I'ma having some troubles installing R under Unix. As it is written in the how to install R the ./configure passed on but when I type make I recieve this error message : gmake[1]: Entering directory `/usr/local/R_HOME/R-2.4.0/m4' gmake[1]: Nothing to be done for `R'. gmake[1]: Leaving

Re: [R] axis command and excel time format

2006-11-08 Thread Carmen Meier
Gabor Grothendieck schrieb: Is the problem how to produce an axis with a given minimum tick, maximum tick and given number of ticks? In that case try this yes but ... ;-) I started with an plain R gui library(zoo) library(chron) # input data # z is from original example mn - times(23:00:00) mx

Re: [R] Extracting parameters for Gamma Distribution

2006-11-08 Thread Thomas Lumley
On Tue, 7 Nov 2006, yongchuan wrote: I'm doing a cox regression with frailty: model - coxph(Surv(Start,Stop,Terminated)~ X + frailty(id),table) I understand that model$frail returns the group level frailty terms. Does this mean this is the average of the frailty values for the respective

[R] Sweave and font problems

2006-11-08 Thread Mark Wardle
Dear All, Having now successfully started using Sweave, I have just noticed an odd side effect with fonts. Using plain LaTeX, I have had no problems using \usepackage{palatino} or \usepackage{times} and the font correctly changes. However, once I convert the tex document into Snw, and run it

Re: [R] axis command and excel time format

2006-11-08 Thread Gabor Grothendieck
My understanding is that the main point of your post was how to get times on the X axis. hopefully at this point its clear how to do that and you can come up with some algorithm to put whatever points you want on. Here is a slight generalization although you will likely have to generalize it

Re: [R] query in R

2006-11-08 Thread Gregor Gorjanc
Jin, Xiaodong Jin close2ceo at yahoo.com writes: how to realize the following SQL command in R? select distinct A, B, count(C) from TABLE group by A, B ; quit; if you want to perform such query on a data.frame! then table(table$a, table$b) will produce similar output, but

[R] interprete wilcox.test results

2006-11-08 Thread Baoqiang Cao
Dear All, I am using wilcox.test to test two samples, data_a and data_b, earch sample has 3 replicates, suppose data_a and data_b are 20*3 matrix. Then I used the following to test the null hypothesis (they are from same distribution.): wilcox.test(x=data_a, y=data_b, alternative=g) I got

Re: [R] get compressed data via a socket connection

2006-11-08 Thread Prof Brian Ripley
?gzcon may help, depending exactly what you mean by `zlib-compressed data'. If not, its code will provide you with a prototype to work on. There are also private entry points in connections.c to (de)compress blocks of data to/from a raw vector which again may be a useful statring point for

Re: [R] Sweave and font problems

2006-11-08 Thread Roger D. Peng
Try \usepackage[noae]{Sweave} -roger Mark Wardle wrote: Dear All, Having now successfully started using Sweave, I have just noticed an odd side effect with fonts. Using plain LaTeX, I have had no problems using \usepackage{palatino} or \usepackage{times} and the font correctly changes.

Re: [R] axis command and excel time format

2006-11-08 Thread Carmen Meier
Gabor Grothendieck schrieb: My understanding is that the main point of your post was how to get times on the X axis. hopefully at this point its clear how to do that and you can come up with some algorithm to put whatever points you want on. That´s right thank you Here is a slight

Re: [R] interprete wilcox.test results

2006-11-08 Thread Torsten Hothorn
On Wed, 8 Nov 2006, Baoqiang Cao wrote: Dear All, I am using wilcox.test to test two samples, data_a and data_b, earch sample has 3 replicates, suppose data_a and data_b are 20*3 matrix. Then I used the following to test the null hypothesis (they are from same distribution.):

Re: [R] axis command and excel time format

2006-11-08 Thread Gabor Grothendieck
You indicated that one of the solutions did not work for you but I am not clear which one you were referring to. To avoid confusion I have deleted all my and your comments except for the following which works for me if I paste it into a fresh session using Edit | Paste commands only menu in

Re: [R] Comparing models in multiple regression and hierarchical linear regression

2006-11-08 Thread Spencer Graves
The questions you ask about the interactions in the model not making sense relates, I believe, to a multiple comparisons issue that is not adequately addressed by the stepAIC analysis you did. To understand this, note first that you've got something close to 2^(2^5) possible models:

[R] Splitting matrix into several small matrices

2006-11-08 Thread Milton Cezar Ribeiro
Hi R-friends I need to read a big matrix (something like 1024000x1024) into 1000 small matrices. After that, I need write each small matrix (1024x1024) into a asc file, putting de index of matrix (1:1000) on the file name. I´m makint it by a hard and unelegant way:

Re: [R] interprete wilcox.test results

2006-11-08 Thread Baoqiang Cao
Thanks Torsten! I got it! Best, Baoqiang Original message Date: Wed, 8 Nov 2006 17:26:18 +0100 (CET) From: Torsten Hothorn [EMAIL PROTECTED] Subject: Re: [R] interprete wilcox.test results To: Baoqiang Cao [EMAIL PROTECTED] Cc: r-help@stat.math.ethz.ch On Wed, 8 Nov 2006,

[R] can Rgui.exe be launched with a script as an argument?

2006-11-08 Thread Bill Bachman
Ultimately, I would like to launch R (from another application) with a script as an argument that creates a plot. I am running on Windows XP with R 2.4.0. I currently can do this with S-Plus from my application by giving this command: “C:\Program Files\Insightful\splus62\cmd\splus.exe”

Re: [R] Install RMySQL with R 2.4.0

2006-11-08 Thread Joe Byers
All, I have tried to compile the RMySQL from source as specified in the documentation. I am receiving the following errors *** beginning of cygwin window [EMAIL PROTECTED] /cygdrive/d/Program Files/R/bin $ Rcmd build --binary ../library/RMySQL *

[R] Numerical Integration

2006-11-08 Thread Xiaofan Cao
Hi everyone, I'm trying to integrate f(x) over x where f(x) does not have a close form but only numerical values at certurn knots of x. Is there a way that I can use any generical R function (such as integrate) or any package to do so? Thanks! I appreciate your time. Best Regards, Martha Cao

Re: [R] multivariate splines

2006-11-08 Thread Spencer Graves
You say, I have looked at a few spline packages in R, but didn't find what I was looking for. Could somebody please point me in the right direction? It's difficult to comment without more information on what you've tried and why you thought it was not adequate. Please provide

Re: [R] axis command and excel time format

2006-11-08 Thread Carmen Meier
Yes this one works, but (sorry) in the OP there was a plot without data to define the range So I tried to use your working suggestion in that manner: library(zoo) library(chron) time - c(2:25:00,2:26:00,2:27:00,2:28:00,2:29:00,2:30:00,2:31:00,

Re: [R] Numerical Integration

2006-11-08 Thread Doran, Harold
You might try the statmod package which provides nodes and weights for gaussian quadrature. -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Xiaofan Cao Sent: Wednesday, November 08, 2006 12:43 PM To: r-help@stat.math.ethz.ch Subject: [R] Numerical

Re: [R] CPU or memory

2006-11-08 Thread Christos Hatzis
Prof. Ripley, Do you mind providing some pointers on how coarse-grained parallelism could be implemented on a Windows environment? Would it be as simple as running two R-console sessions and then (manually) combining the results of these simulations. Or it would be better to run them as batch

Re: [R] nls

2006-11-08 Thread Jost Burkardt
Xiaodong Jin [EMAIL PROTECTED] writes: y [1] 1 11 42 64 108 173 214 t [1] 1 2 3 4 5 6 7 nls(1/y ~ c*exp(-a*b*t)+1/b, start=list(a=0.001,b=250,c=5), trace=TRUE) 29.93322 :0.001 250.000 5.000 Error in numericDeriv(form[[3]], names(ind), env) : Missing

Re: [R] can Rgui.exe be launched with a script as an argument?

2006-11-08 Thread Lanre Okusanya
I assume this is windoze. Have you considered using the Rterm? You can have in your script the plot to be created as an eps file Lanre Bill Bachman wrote: Ultimately, I would like to launch R (from another application) with a script as an argument that creates a plot. I am running on Windows

Re: [R] Numerical Integration

2006-11-08 Thread Gabor Grothendieck
You could approximate it with splines and then integrate that: x - 0:10/10 f - function(x) 1 + x + x^2 + x^3 + x^4 y - f(x) fs - splinefun(x, y) integrate(fs, 0, 1) integrate(f, 0, 1) On 11/8/06, Xiaofan Cao [EMAIL PROTECTED] wrote: Hi everyone, I'm trying to integrate f(x) over x where f(x)

Re: [R] CPU or memory

2006-11-08 Thread Prof Brian Ripley
On Wed, 8 Nov 2006, Christos Hatzis wrote: Prof. Ripley, Do you mind providing some pointers on how coarse-grained parallelism could be implemented on a Windows environment? Would it be as simple as running two R-console sessions and then (manually) combining the results of these

Re: [R] axis command and excel time format

2006-11-08 Thread Gabor Grothendieck
Your code plots x which has nothing to do with xt. On 11/8/06, Carmen Meier [EMAIL PROTECTED] wrote: Yes this one works, but (sorry) in the OP there was a plot without data to define the range So I tried to use your working suggestion in that manner: library(zoo) library(chron) time -

Re: [R] Install RMySQL with R 2.4.0

2006-11-08 Thread Joe Byers
All, After staring at this error message for an hour or so yesterday and this morning. I decided to try something else. Low and behold trying to build the package in cygwin causes R to try and build under linux/unix not windows. I went to the command prompt and was able to build the

Re: [R] CPU or memory

2006-11-08 Thread Christos Hatzis
Great. I will try it. Thank you. -Christos -Original Message- From: Prof Brian Ripley [mailto:[EMAIL PROTECTED] Sent: Wednesday, November 08, 2006 1:21 PM To: Christos Hatzis Cc: 'Stefan Grosse'; r-help@stat.math.ethz.ch; 'Taka Matzmoto' Subject: RE: [R] CPU or memory On Wed, 8 Nov

Re: [R] get compressed data via a socket connection

2006-11-08 Thread Jeffrey Horner
Prof Brian Ripley wrote: ?gzcon may help, depending exactly what you mean by `zlib-compressed data'. If not, its code will provide you with a prototype to work on. There are also private entry points in connections.c to (de)compress blocks of data to/from a raw vector which again may be a

Re: [R] have I an actual matrix?

2006-11-08 Thread Ricardo Rodríguez - Your EPEC ICT Team
On Nov 8, 2006, at 11:12 AM, Patrick Burns wrote: My guess is that you have a data frame and not a matrix. Try barplot(as.matrix(rs)) See chapter 1 of S Poetry for an explanation of data structures in R. Thanks, Patrick, You are completely right. I've been trying to plot a data frame.

Re: [R] can Rgui.exe be launched with a script as an argument?

2006-11-08 Thread Hans-Peter
2006/11/8, Lanre Okusanya [EMAIL PROTECTED]: I assume this is windoze. Have you considered using the Rterm? You can windows... afaic try with R.exe or, as said and imho better, Rterm.exe. Examples testR.bat: set cmd=C:\Programme\R\R-2.3.1\Bin\R.exe call %cmd% CMD BATCH --slave --vanilla

Re: [R] have I an actual matrix?

2006-11-08 Thread Ricardo Rodríguez - Your EPEC ICT Team
On Nov 8, 2006, at 11:41 AM, Michael Kubovy wrote: On Nov 8, 2006, at 4:44 AM, Ricardo Rodríguez wrote: Hi Ricardo, You probably have two columns in 'rs'. You need to do the barplot on one of them and use the other as the vector of labels. Assuming that the first column is n and the

[R] Failing to install on Linux FC4

2006-11-08 Thread Benilton Carvalho
Hi everyone, I downloaded the source code available at: http://cran.fhcrc.org/src/base/R-2/R-2.4.0.tar.gz to a linux machine (Linux 2.6.11-1.1369_FC4smp #1 SMP). I successfully configured and compiled it, which means I'm able to run it from the R-2.4.0/bin directory. I want to do a system

Re: [R] can Rgui.exe be launched with a script as an argument?

2006-11-08 Thread Bill Bachman
Thank you to Paul Murrell for a solution: # test.R require(graphics) a-c(1,2,3,4) b-c(1,4,9,16) plot(a,b) -Original Message- From: Lanre Okusanya [mailto:[EMAIL PROTECTED] Sent: Wednesday, November 08, 2006 1:04 PM To: Bill Bachman Cc: r-help@stat.math.ethz.ch Subject: Re: [R] can

[R] size restrictions ??

2006-11-08 Thread Bond, Stephen
should I be able to read in a table of 1mln by 100 into R or it's too big? Thank you very much stephen [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

[R] Reading bzip2 and tar files in R

2006-11-08 Thread DrakeGis
Is there any way to read and decompress inside R tar and bzip2 files ? Thanks D - Stay ahead of the information curve. Receive EDA news and jobs on your desktop daily. Subscribe today to the EDA CafeNews newsletter. [

[R] Making a case for using R in Academia

2006-11-08 Thread Charilaos Skiadas
Hello, new to the list, first message. This question perhaps might be more appropriate to R-sig-teaching, and I'd be happy to take it there if this is not the right place for it. I am teaching applied statistics at a small liberal arts college with limited resources, and we are currently

Re: [R] axis command and excel time format

2006-11-08 Thread Carmen Meier
Gabor Grothendieck schrieb: Your code plots x which has nothing to do with xt. The same result if you change xt to x: 02:25 at the origin nothing else - I do not know why #-- your suggestion mn - times(min_time) mx - times(max_time) n - 12 t - times(seq(mn, mx, length =

[R] superpc.plotcv

2006-11-08 Thread Haleh Yasrebi
Hi, Why do I get error msg by calling superpc.plotcv (see below)? I use superpc v. 1.0.4 and R 2.2.1 thanks Haleh -- set.seed(332) x-matrix(rnorm(1000*40),ncol=40) y-10+svd(x[1:60,])$v[,1]+ .1*rnorm(40) censoring.status-

Re: [R] multivariate splines

2006-11-08 Thread Tamas K Papp
On Wed, Nov 08, 2006 at 09:49:50AM -0800, Spencer Graves wrote: You say, I have looked at a few spline packages in R, but didn't find what I was looking for. Could somebody please point me in the right direction? It's difficult to comment without more information on what you've tried

Re: [R] question on multilevel modeling

2006-11-08 Thread David Duffy
Christine A. Calmes [EMAIL PROTECTED] wrote: [...] Basically, I am trying to test a psychological model in which a person's level of passive, negative communication with others in their environment (CORUMTO) on one week predicts their level of depression on the following week (BDIAFTER)

[R] a list with row.names = c(NA, 1000)

2006-11-08 Thread Taka Matzmoto
R users, I posted a message on error message for row names. I looked at a saved file that has a list structure. I found that some of list components have row.names = c(NA, 1000). I compared what I did with R.2.3.0 with what I did with R.2.4.0. Previously, with R.2.3.0, row.names = c(1, 2, 3,

Re: [R] Numerical Integration

2006-11-08 Thread Xiaofan Cao
Hi Ravi and Harold, Thanks for the input. I'm using trapezoidal rule and like to know if there's other alternatives. This f(x) is the kernel density estimator and thus we can get an estimate of f(x) at any given x in theory. Thanks again, Martha On Wed, 8 Nov 2006, Ravi Varadhan wrote:

Re: [R] CPU or memory

2006-11-08 Thread John C Frain
I would like to thank all who replied to my question about the efficiency of various cpu's in R. Following the advice of Bogdan Romocea I have put a sample simulation and the latest version of R on a USB drive and will go to a few suppliers to try it out. I will report back if I find anything of

Re: [R] Reading bzip2 and tar files in R

2006-11-08 Thread Roger D. Peng
bzip2-ped files can be read via the bzfile() connection function. I don't think there's a way to read from tar files though. -roger DrakeGis wrote: Is there any way to read and decompress inside R tar and bzip2 files ? Thanks D - Stay

Re: [R] axis command and excel time format

2006-11-08 Thread Gabor Grothendieck
Please provide a complete self contained example. I can't follow the partial code below; however, its likely you are plotting one thing and creating axes using another so there is no reason it should come out right. On 11/8/06, Carmen Meier [EMAIL PROTECTED] wrote: Gabor Grothendieck schrieb:

Re: [R] Sweave and font problems

2006-11-08 Thread Mark Wardle
Roger D. Peng wrote: Try \usepackage[noae]{Sweave} -roger Works perfectly - many thanks! Best wishes, Mark -- Dr. Mark Wardle Clinical research fellow and Specialist Registrar in Neurology, C2-B2 link, Cardiff University, Heath Park, CARDIFF, CF14 4XN. UK

Re: [R] error message for row names

2006-11-08 Thread Prof Brian Ripley
On Wed, 8 Nov 2006, Taka Matzmoto wrote: Hi R users - I got an error message when reading in a saved file (a list structure) dget(REPLICATION001) Error in attributes(.Data) - c(attributes(.Data), attrib) : row names must be 'character' or 'integer', not 'double' I don't need row

[R] Cox model using mfp library

2006-11-08 Thread Inman, Brant A. M.D.
I am running R 2.3.1 on a Windows XP machine. I have a large dataset of over 13 000 cases of a disease for which I am attempting to build a prognostic model using Cox proportional hazards regression. Some of the continuous covariates are skewed and therefore require transformation for use in

[R] convert factor p000345 to numeric

2006-11-08 Thread Marco Boks
Dear All, I am lost about the following. I have got a large dataframe (largeset) with in the first column identification numbers as factors largeset$ID p000345 p000356 p000569 etc -- in order to use them to merge with another dataframe with numerical values (000345, 000356) I

Re: [R] a list with row.names = c(NA, 1000)

2006-11-08 Thread Peter Dalgaard
Taka Matzmoto [EMAIL PROTECTED] writes: R users, I posted a message on error message for row names. I looked at a saved file that has a list structure. I found that some of list components have row.names = c(NA, 1000). I compared what I did with R.2.3.0 with what I did with R.2.4.0.

Re: [R] Comparing models in multiple regression and hierarchical linear regression

2006-11-08 Thread Andrew Robinson
Hello Jenifer, your question reflects some very interesting statistical problems, surrounding the effects of subset selection upon estimation and inference. There does not seem to be a right answer, as far as I am aware, so I will offer an opinion. I'd welcome further discussion. I don't have

Re: [R] Making a case for using R in Academia

2006-11-08 Thread Mitchell Maltenfort
Check http://cran.r-project.org/doc/FAQ/R-FAQ.html#R-and-S and following links. R has a 95% overlap with S and S+, and those two are popular enough that statistics books target them (e.g., Venables and Ripley). I am teaching applied statistics at a small liberal arts college with limited

Re: [R] convert factor p000345 to numeric

2006-11-08 Thread Marc Schwartz
On Wed, 2006-11-08 at 23:16 +0100, Marco Boks wrote: Dear All, I am lost about the following. I have got a large dataframe (largeset) with in the first column identification numbers as factors largeset$ID p000345 p000356 p000569 etc in order to use them to merge with another

Re: [R] Install RMySQL with R 2.4.0

2006-11-08 Thread Frank McCown
Joe Byers wrote: All, After staring at this error message for an hour or so yesterday and this morning. I decided to try something else. Low and behold trying to build the package in cygwin causes R to try and build under linux/unix not windows. I went to the command prompt and was

Re: [R] axis command and excel time format

2006-11-08 Thread Robert Baer
Carmen, Gabor has already given you the detail you ask for, but might try the following plot to see what is going wrong: plot(times(tt), x, type='l') This does not give you the EXACT control of the axis you asked for, but this simple plot command gives you a fairly nice result. It

Re: [R] Making a case for using R in Academia

2006-11-08 Thread Marc Schwartz
On Wed, 2006-11-08 at 15:55 -0500, Charilaos Skiadas wrote: Hello, new to the list, first message. This question perhaps might be more appropriate to R-sig-teaching, and I'd be happy to take it there if this is not the right place for it. I am teaching applied statistics at a small

[R] Trying to create multiple stripcharts

2006-11-08 Thread Marc Feuerstein
Hi everybody, I'm a new member of the list. I'm just starting to use R. I'm looking for information about stripcharts. I didn't find a lot on the web. Outside http://www.maths.lth.se/help/R/.R/library/graphics/html/stripchart.html there isn't much I found. I'm trying to create three vertical

[R] parameter estimation -- Censored mean to real mean

2006-11-08 Thread Jenny Stadt
Hi List, Is there anyway to estimate the parameters with integral forms? Here is the function. The idea is to calculate the real mean from the censored mean given normal distribution. I have the function as trun.mean- function(mu,sigma,t) # t is the threshold { f - function(x)

Re: [R] Trying to create multiple stripcharts

2006-11-08 Thread Deepayan Sarkar
On 11/8/06, Marc Feuerstein [EMAIL PROTECTED] wrote: Hi everybody, I'm a new member of the list. I'm just starting to use R. I'm looking for information about stripcharts. I didn't find a lot on the web. Outside http://www.maths.lth.se/help/R/.R/library/graphics/html/stripchart.html there

[R] Can I set the proxy server and authentification from within R?

2006-11-08 Thread Werner Wernersen
Hi, I am using Tinn-R and somehow I cannot point it to start a R shortcut with http_proxy settings including authentification. It always shows the .exe. Is there any way to set the proxy from within R or any other ideas how to do it? I am using R 2.4 for Windows. Thanks a million, Werner

Re: [R] Comparing models in multiple regression and hierarchical linear regression

2006-11-08 Thread hadley wickham
On the other hand, at least one _other_ authoritative statistician I know claims that he will never bother to test a term that he can't interpret - he usually draws the line at three-way interactions. So, opinions are divided among the authorities. This leaves you very vulnerable to your

Re: [R] Comparing models in multiple regression and hierarchical linear regression

2006-11-08 Thread Andrew Robinson
Good point, Hadley. But for the purpose of discussion, then, what would you think about 5-way interactions? Or ten-way interactions? Surely one has to draw the line somewhere. I suppose that a generalization of my summary of that position would be: ask yourself if you can interpret a term,

Re: [R] CPU or memory

2006-11-08 Thread Liaw, Andy
My understanding is that it doesn't have much to do with 32- vs. 64-bit, but what the instruction sets of the CPUs. If I'm not mistaken, at the same clock speed, a P4 would run slower than PIII simply because P4 does less per clock-cycle. Also, I believe for the same architecture, single core

Re: [R] Install RMySQL with R 2.4.0

2006-11-08 Thread Joe W. Byers
Frank, I will post a zip file of the RMySQL 5.10 on my university web site www.cba.utulsa.edu/byersj under the Research and Analytics link in the next day or two. Here at the University we had a major technology problem in the college of business causing IT to move up a migration time table

Re: [R] Which genetic optimization package allows customized crossover and mutation operation

2006-11-08 Thread Spencer Graves
Are you familiar with www.bioconductor.org? They have a listserve with people who may know more about the question you asked. Spencer Graves sun wrote: Hi, I am looking for genetic optimization package that allow to define my own chromosome solution and crossover/mutation

Re: [R] Install RMySQL with R 2.4.0

2006-11-08 Thread Joe W. Byers
Frank, While I was replying to you our IT guys sent me a message as said I could connect. I have posted the RMySQL zip file at http://bus.cba.utulsa.edu/byersj/Research.asp Good luck to all Joe Frank McCown wrote: Joe Byers wrote: All, After staring at this error message for an hour or

Re: [R] Install RMySQL with R 2.4.0

2006-11-08 Thread Joe W. Byers
Frank, While I was replying to you our IT guys sent me a message as said I could connect. I have posted the RMySQL zip file at http://bus.cba.utulsa.edu/byersj/Research.asp Good luck to all Joe Frank McCown wrote: Joe Byers wrote: All, After staring at this error message for an hour or

[R] Repeated Measures MANOVA in R

2006-11-08 Thread A. Bolu Ajiboye
Can R do a repeated measures MANOVA and tell what dimensionality the statistical variance occupies? I have been using MATLAB and SPSS to do my statistics. MATLAB can do ANOVAs and MANOVAs. When it performs a MANOVA, it returns a parameter d that estimates the dimensionality in which the means

Re: [R] Can I set the proxy server and authentification from within R?

2006-11-08 Thread Prof Brian Ripley
On Thu, 9 Nov 2006, Werner Wernersen wrote: I am using Tinn-R and somehow I cannot point it to start a R shortcut with http_proxy settings including authentification. It always shows the .exe. Is there any way to set the proxy from within R or any other ideas how to do it? From the help

Re: [R] axis command and excel time format

2006-11-08 Thread Carmen Meier
Gabor Grothendieck schrieb: Please provide a complete self contained example. I can't follow the partial code below; however, its likely you are plotting one thing and creating axes using another so there is no reason it should come out right. You are right .. seems to be that it was too

[R] Meta-regression with lmer() ? If so, how ?

2006-11-08 Thread Emmanuel Charpentier
Dear List, I am (again) looking at meta-regression as a way to refine meta-analytic results. What I want to do is to assess the impact of some fixed factors on the results of a meta-analysis. Some of them may be crossed with the main factor of the meta-analysis (e. g. clinical presentation of a