Re: [R] Saving Graphics

2005-09-30 Thread vincent
Mike Jones a écrit : Is there a way to save graphs from the graphics window using commands in the R console? you may try : bmp(myimage.bmp); plot(...); dev.off(); # png(), jpeg() also available hih __ R-help@stat.math.ethz.ch mailing list

Re: [R] Bug in lmer?

2005-09-30 Thread Martin Maechler
Mark == Mark Lyman [EMAIL PROTECTED] on Thu, 29 Sep 2005 14:44:38 -0600 writes: Mark I am relatively new to R so I am not confident enough in what I am doing Mark to be certain this is a bug. Mark I am running R 2.1.1 on a Windows XP Mark machine and the lme4 package

Re: [R] Bug in lmer?

2005-09-30 Thread Martin Henry H. Stevens
To All: Mark and Martin's code ran fine on R : Version 2.1.1 (2005-06-20) on a Macintosh Dual 2 GHz G5 with 1.5 GB DDR SDRAM. Hank On Sep 30, 2005, at 3:21 AM, Martin Maechler wrote: source(ftp://stat.ethz.ch/U/maechler/R/mltloc-ex.R;, echo = TRUE) Dr. Martin Henry H. Stevens, Assistant

[R] cox proportional-hazards regress for interval censor data

2005-09-30 Thread Lixia ZHU
Hi. I used coxph(surv(start,end,event)~~event,data) to deal with interval censor data. Does anyone know similar samples using coxph(surv(start,end,event)~~event,data)? If you knows, can you tell me? I'll really appreciate it. Thank you very much R learner.

[R] List Email Statistic

2005-09-30 Thread José Raul Capablanca
Dear All, How can I can to know a mail list , to speak about exclusively statistic. Thanks. __ Espacio para todos tus mensajes, antivirus y antispam ¡gratis! __ R-help@stat.math.ethz.ch mailing list

Re: [R] Bug in lmer?

2005-09-30 Thread Martin Henry H. Stevens
I forgot to mention the OS - Mac OS 10.4.2 (Tiger with the latest update). {:-) Hank On Sep 30, 2005, at 3:40 AM, Martin Henry H. Stevens wrote: To All: Mark and Martin's code ran fine on R : Version 2.1.1 (2005-06-20) on a Macintosh Dual 2 GHz G5 with 1.5 GB DDR SDRAM. Hank On Sep 30,

Re: [R] List Email Statistic

2005-09-30 Thread Gabor Grothendieck
There is some discussion and data sources of the volume of email on the list in: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/27532.html On 9/30/05, José Raul Capablanca [EMAIL PROTECTED] wrote: Dear All, How can I can to know a mail list , to speak about exclusively statistic. Thanks.

Re: [R] List Email Statistic

2005-09-30 Thread Jim Lemon
José Raul Capablanca wrote: Dear All, How can I can to know a mail list , to speak about exclusively statistic. Thanks. Hola Jose, If you mean a list devoted exclusively to statistics (Si desea una lista solo para la statistica) http://www.jiscmail.ac.uk/lists/allstat.html (y una lista

[R] Searching for specific functions

2005-09-30 Thread Fabrizio Ferri Benedetti
Hello everybody (I'm new here) I am looking for 2 specific functions. The first one is a variance formula that does not divide by n-1 (I don't want the unbiased formula, that's it). In case it does not exist, how do I modify the default function code? Yes, a workaround would be to multiply

Re: [R] List Email Statistic

2005-09-30 Thread Ted Harding
On 30-Sep-05 Jim Lemon wrote: José Raul Capablanca wrote: Dear All, How can I can to know a mail list , to speak about exclusively statistic. Thanks. Hola Jose, If you mean a list devoted exclusively to statistics (Si desea una lista solo para la statistica)

[R] Compare predict and experimental values

2005-09-30 Thread Toni Viúdez
Hi everybody: I just generate different maps of interpolation with different methods but with the same package -(gstat, i used the krige command). And now for evaluate the best method i want to compare the predict values versus experimental, but i don't know how do it. Could anybody tell me how

Re: [R] Easy cut paste from Excel to R

2005-09-30 Thread stefan . albrecht
Maybe this helps, too. readExcel - function(row.names = 1, check.names = FALSE, ...) read.delim(clipboard, row.names = row.names, check.names = check.names, ...) writeExcel - function(x, check.rows = FALSE, check.names = FALSE, ...){ df - data.frame(x, check.rows = check.rows, check.names

Re: [R] Compare predict and experimental values

2005-09-30 Thread Roger Bivand
On Fri, 30 Sep 2005, Toni Viúdez wrote: Hi everybody: I just generate different maps of interpolation with different methods but with the same package -(gstat, i used the krige command). And now for evaluate the best method i want to compare the predict values versus experimental, but i

[R] Creating an array [0 to 101] of a list

2005-09-30 Thread Rainer M. Krug
Hi I looked, but I didn't find it: I need an array [0 to 101] where each element is a list (the result of Kest in spatstat). I.e. I want to do: A[2] - Kest(pp1) A[3] - Kest(pp2) ... A[101] - Kest(pp100) How can I create A ? Rainer -- Rainer M. Krug, Dipl. Phys. (Germany), MSc Conservation

Re: [R] Creating an array [0 to 101] of a list

2005-09-30 Thread TEMPL Matthias
Probably, you could do it like A - list() for( i in 1:101){ A[[i]] - Kest( get( paste(pp,i,sep=) ) ) } Best, Matthias Hi I looked, but I didn't find it: I need an array [0 to 101] where each element is a list (the result of Kest in spatstat). I.e. I want to do: A[2] - Kest(pp1)

Re: [R] Creating an array [0 to 101] of a list

2005-09-30 Thread Rainer M. Krug
Thanks a lot Matthias, Rainer TEMPL Matthias wrote: Probably, you could do it like A - list() for( i in 1:101){ A[[i]] - Kest( get( paste(pp,i,sep=) ) ) } Best, Matthias Hi I looked, but I didn't find it: I need an array [0 to 101] where each element is a list (the result of

Re: [R] Creating an array [0 to 101] of a list

2005-09-30 Thread Christoph Lehmann
a - array(vector(list, 3), dim=c(3)) a[[1]] - list(x = 1, y = 0, z = -1) a[[2]] - list(x = 0, y = 1, z = -1) a[[3]] - list(x = 0, y = -1, z = 0) HTH christoph Rainer M. Krug wrote: Hi I looked, but I didn't find it: I need an array [0 to 101] where each element is a list (the result of

Re: [R] R-help Digest, Vol 31, Issue 30

2005-09-30 Thread John Maindonald
With lme4, use of mcmcsamp can be insightful. (Douglas Bates drew my attention to this function in a private exchange of emails.) The distributions of random effects are simulated on a log scale, where the distributions are much closer to symmetry than on the scale of the random effects

Re: [R] Bug in lmer?

2005-09-30 Thread Renaud Lancelot
I got a crash too with my Win XP config running: platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status beta major2 minor2.0 year 2005 month09 day 24 svn rev 35666 language R Package: Matrix Version: 0.98-7 Date:

[R] correlation question

2005-09-30 Thread vincent
Dear R-helpers, First, double apologies because the question is not directly related to R, and also probably quite simple. I have observed (with many data) that given 2 series X and Y cor(X,Y) and cor(log(X), log(Y)) differs only very slightly. I suspect this is because log(x) is monotonous, ...

Re: [R] List Email Statistic

2005-09-30 Thread Marc Schwartz
On Thu, 2005-09-29 at 23:47 -0500, José Raul Capablanca wrote: Dear All, How can I can to know a mail list , to speak about exclusively statistic. Thanks. If you have access to Usenet either via an NNTP server or via Google Groups, there are three principal groups for general statistics

[R] Clear Console

2005-09-30 Thread Leonardo L Miceli
Hi, What is the function to clear the console? tks [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: [R] Easy cut paste from Excel to R?

2005-09-30 Thread roger bos
Thanks Gabor and Petr. It seems that all you have to do is use 'clipboard' and not 'file(clipboard)' as we were doing. The row.names won't copy if you use file. I don't know if that is a bug or a feature, but if you were debugging someones code who already had 'file()' there, it would be hard to

Re: [R] correlation question

2005-09-30 Thread TEMPL Matthias
First: The question is really not related to R. What you can find very easly in books or in the internet: The correlation coefficient is invariant under a linear transformation (and the proof) #E.g. library(rrcov) data(brain) cor(brain) cor(scale(brain)) #Is the same, BUT cor(log(brain))

Re: [R] Clear Console

2005-09-30 Thread Doran, Harold
Ctrl L -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Leonardo L Miceli Sent: Friday, September 30, 2005 8:36 AM To: r-help@stat.math.ethz.ch Subject: [R] Clear Console Hi, What is the function to clear the console? tks [[alternative HTML

Re: [R] Estimate predictor contribution in GAM models

2005-09-30 Thread Simon Wood
On Tue, 20 Sep 2005, Yves Magliulo wrote: hi, i'm using gam() function from package mgcv. if G is my gam object, then SG=summary(G) Formula: y ~ +s(x0, k = 5) + s(x1) + s(x2, k = 3) Parametric coefficients: Estimate std. err.t ratioPr(|t|) (Intercept) 3.462e+07

Re: [R] Smooth terms significance in GAM models

2005-09-30 Thread Simon Wood
i'm using gam() function from package mgcv with default option (edf estimated by GCV). G=gam(y ~ s(x0, k = 5) + s(x1) + s(x2, k = 3)) SG=summary(G) Formula: y ~ +s(x0, k = 5) + s(x1) + s(x2, k = 3) Parametric coefficients: Estimate std. err.t ratioPr(|t|)

Re: [R] Descriptive statistics for tables

2005-09-30 Thread DAVID CAMACHO
Thanks for your response, What I need is extremely simple. And I suppose there are so many way to do it. But as I have so many files to do it, I am looking for the simplest way (if possible) What I meant was that I have tables with the same numbers of rows and columns, (square form, should I

Re: [R] Bug in lmer?

2005-09-30 Thread Mark Lyman
I am using version 0.98-7 of the Matrix package. I used the RGui Install Packages... menu option to get the lme4 package from CRAN and this version of the Matrix was automatically downloaded as well. Martin Maechler wrote: Mark == Mark Lyman [EMAIL PROTECTED] on Thu, 29 Sep 2005 14:44:38

Re: [R] Descriptive statistics for tables

2005-09-30 Thread Jean Eid
I do not totally understand your question as well. You seem to want a descriptive statistic about a unitary number. What is the sd of a number? or any other descriptive statictic. Maybe you mean for the columns or rows or it could be that these are t-stats or z-stats that you need to get p_values

Re: [R] correlation question

2005-09-30 Thread vincent
TEMPL Matthias a écrit : First: The question is really not related to R. What you can find very easly in books or in the internet: The correlation coefficient is invariant under a linear transformation (and the proof) yes. #E.g. library(rrcov) data(brain) cor(brain) cor(scale(brain))

[R] p-value for non-linear variable in overdispersed glm()

2005-09-30 Thread tsnall
Dear all, I am fitting an nonlinear glm() using optim() by first minimising glm(resp~ var1 + var2, family=binomial, data=data)$deviance where var1= exp(-a1*dist1), and var2= exp(-a2*dist2), where a1 and a2 are parameters and dist1 and dist2 are independent variables. Next, I calculate the

Re: [R] Searching for specific functions

2005-09-30 Thread Adaikalavan Ramasamy
1) I am not sure why you want the biased version but here are two ways to do this. Also if you have large enough sample size, it would not matter much. my.var1 - function(x) mean( (x - mean(x))^2 ) my.var2 - function(x) mean( x^2 - mean(x)^2 ) 2) Look at quantiles function. Your

[R] ast package?

2005-09-30 Thread Rabaa, Maia
According to the extremely helpful reference card for time series analysis provided at: http://cran.r-project.org/doc/contrib/Ricci-refcard-ts.pdf, the ast package is necessary to perform some of the functions. Where can this package be installed from, as I cannot find it from my pulldown list,

[R] Compiling R on OpenSolaris

2005-09-30 Thread Vincent Yau
Hi: I am trying to R (v2.1.1) to compile on OpenSolaris (build 22) using gcc version 3.4.3. But I am getting this error message: gmake[5]: Entering directory `/usr/local/R-2.1.1/src/library/tools/src' ../../../../library/tools/libs/tools.so is unchanged gmake[5]: Leaving directory

Re: [R] ast package?

2005-09-30 Thread Peter Dalgaard
Rabaa, Maia [EMAIL PROTECTED] writes: According to the extremely helpful reference card for time series analysis provided at: http://cran.r-project.org/doc/contrib/Ricci-refcard-ts.pdf, the ast package is necessary to perform some of the functions. Where can this package be installed from,

Re: [R] Descriptive statistics for tables

2005-09-30 Thread Dave Roberts
If I understand the request, he wants to take a large number of matrices of identical size and stack them into a three dimentional array, and then calculate statistics on the the third dimension. If the multiple arrays have object names they can be combined into a 3-d array a -

[R] ACCESS R and dates

2005-09-30 Thread sloan jones
I have used the RODBC package to read in data I have stored in an Access file. When I am using data from files other than ACCESS I have no problem using the survival package to work with dates; however, with the ACCESS data the dates are reading-in in the following format: 2004-02-11 Pacific

[R] Dots in function names

2005-09-30 Thread Mike Prager
Recent R function names seem to be using CaseOfTheLetters to mark words rather than dots as was done previously. Is the use of dots in function names deprecated, or is that simply a style choice? Will function names with dots cause problems in future revisions? Mike Prager

Re: [R] Dots in function names

2005-09-30 Thread Thomas Lumley
On Fri, 30 Sep 2005, Mike Prager wrote: Recent R function names seem to be using CaseOfTheLetters to mark words rather than dots as was done previously. Is the use of dots in function names deprecated, or is that simply a style choice? Will function names with dots cause problems in future

[R] gbm package on Sun sparc

2005-09-30 Thread Wong, Jackson
Hi, Does anyone know how to compile the gbm-1.5-1 package on a Sun sparc platform running Solaris 9? I have been using gcc-3.4.2 and getting an NAN undeclared error. Any help would be great. thanks, Jackson ___ This

[R] .C help

2005-09-30 Thread Elizabeth Lawson
Hi, I am hoping some one can help me. I am learning to use C and would like to learn how to call c code in R. I have look at Writing R Extensions and I tried to copy the example on page 38 void convolve(double *a, int *na, double *b, int *nb, double *ab) { int i, j, nab = *na +

Re: [R] Dots in function names

2005-09-30 Thread Peter Dalgaard
Mike Prager [EMAIL PROTECTED] writes: Recent R function names seem to be using CaseOfTheLetters to mark words rather than dots as was done previously. Is the use of dots in function names deprecated, or is that simply a style choice? Will function names with dots cause problems in future

[R] .C help

2005-09-30 Thread Elizabeth Lawson
Hi, I am hoping some one can help me. I am learning to use C and would like to learn how to call c code in R. I have look at Writing R Extensions and I tried to copy the example on page 38 void convolve(double *a, int *na, double *b, int *nb, double *ab) { int i, j, nab = *na +

[R] C program

2005-09-30 Thread Erin Hodgess
Dear R People: I have R Version 2.1.1. for Windows in binary form. I would like to look at the C code for massdist. It is part of the density function. How would I access this, please? Thank you very much! Sincerely, Erin Hodgess Associate Professor Department of Computer and Mathematical

Re: [R] C program

2005-09-30 Thread Roger Bivand
On Fri, 30 Sep 2005, Erin Hodgess wrote: Dear R People: I have R Version 2.1.1. for Windows in binary form. I would like to look at the C code for massdist. It is part of the density function. How would I access this, please? https://svn.r-project.org/R/trunk/src/appl/massdist.c is

Re: [R] C program

2005-09-30 Thread Thomas Lumley
On Fri, 30 Sep 2005, Erin Hodgess wrote: Dear R People: I have R Version 2.1.1. for Windows in binary form. I would like to look at the C code for massdist. It is part of the density function. The .C() call shows that massdist is in the base package. This means that it is likely to be

Re: [R] Dots in function names

2005-09-30 Thread Mike Prager
on 9/30/2005 2:55 PM Peter Dalgaard said the following: Mike Prager [EMAIL PROTECTED] writes: Recent R function names seem to be using CaseOfTheLetters to mark words rather than dots as was done previously. Is the use of dots in function names deprecated, or is that simply a style choice?

Re: [R] .C help

2005-09-30 Thread Roger Bivand
On Fri, 30 Sep 2005, Elizabeth Lawson wrote: Hi, I am hoping some one can help me. I am learning to use C and would like to learn how to call c code in R. I have look at Writing R Extensions and I tried to copy the example on page 38 void convolve(double *a, int *na,

[R] how to impose the lines (Apply?)

2005-09-30 Thread Lisa Wang
Hello, I have data like the following (38 patients) each have a lot of record (some have more than others). I would like to plot time* pressure for each patient in the same plot. I try to write a function and then split the data by patients id and then sapply the function to the splited list, but

Re: [R] .C help

2005-09-30 Thread Tuszynski, Jaroslaw W.
See http://www.stats.uwo.ca/faculty/murdoch/software/compilingDLLs/ and R Installation and Administration manual for details that will have to be followed very precisely. My own 2 cents below should be taken as fuzzy recollections, compared to detailed directions in 2 sources above. I do not know

Re: [R] how to impose the lines (Apply?)

2005-09-30 Thread Uwe Ligges
Lisa Wang wrote: Hello, I have data like the following (38 patients) each have a lot of record (some have more than others). I would like to plot time* pressure for each patient in the same plot. I try to write a function and then split the data by patients id and then sapply the function

[R] 'mean' function

2005-09-30 Thread Lisa Wang
Hello, Could you please let me know how to see the mean function in R . The following is what I see when type in mean and enter mean function (x, ...) UseMethod(mean) environment: namespace:base I would like to see how the function is writen. Thanks a lot Lisa Wang Toronto, Ca This e-mail

Re: [R] 'mean' function

2005-09-30 Thread Duncan Murdoch
On 9/30/2005 5:57 PM, Lisa Wang wrote: Hello, Could you please let me know how to see the mean function in R . The following is what I see when type in mean and enter mean function (x, ...) UseMethod(mean) environment: namespace:base I would like to see how the function is writen.

Re: [R] 'mean' function

2005-09-30 Thread Søren Højsgaard
Try base::mean.default Søren Fra: [EMAIL PROTECTED] på vegne af Lisa Wang Sendt: fr 30-09-2005 23:57 Til: R-Help Emne: [R] 'mean' function Hello, Could you please let me know how to see the mean function in R . The following is what I see when type in mean

[R] How to get to the varable in a list

2005-09-30 Thread Lisa Wang
Hello, I have a list lis as the following: $1 x1 x2 4 3 1 $2 x1 x2 3 3 2 5 3 2 $3 x1 x2 2 3 3 6 3 3 How do I get the x1 varible? for example ss1 This e-mail may contain confidential and/or privileged information for the sole use of the intended recipient. Any review or

Re: [R] nonparametric 2way repeated-measures anova

2005-09-30 Thread Chris Buddenhagen
Dear R Christoph Lehmann asked about this recently- as seen at http://finzi.psych.upenn.edu/R/Rhelp02a/archive/53191.html but the reply was not detailed. Did Christoph get a solution to this? Can anyone point me toward how to do repeated measures comparison using the the friedman test or

Re: [R] Bug in lmer?

2005-09-30 Thread Horacio Montenegro
Just to add more info: I've got the same error with both lme4 0.95-10 Matrix 0.95-13 and lme4 0.98-1 Matrix 0.98-7, always running on Win98SE, R 2.1.1. So it seems to occur with any Windows version. cheers, Horacio I'm putting the data and an R script up for FTP, so

Re: [R] 'mean' function

2005-09-30 Thread Martin Lam
Type ?mean and hit enter in the console to see the description of the function. # randomly draw 5 numbers between 1 and 10 without replacement randomlist = sample(1:10, 5) # calculate the mean of the list meanoflist = mean(randomlist) # show the result meanoflist HTH, Martin --- Lisa Wang

[R] update MASS

2005-09-30 Thread Ramón Casero Cañas
I'd like to update MASS from version 7.2-11 to version 7-2.19. I am running R 2.0.1 on ubuntu, that installs MASS with the package r-cran-vr. I have tried doing update.packages(), and I get a list of packages that I could update, but none is called MASS. I have tried the following too, right

[R] Inverse autocorrelation function?

2005-09-30 Thread David Hartley
In time series analysis it is helpful to plot the autocorrelation function (ACF), partial autocorrelation function (PACF), and the inverse autocorrelation function (IACF). The stats library provides the ability to compute and plot the ACF and PACF, but I cannot find an [R] procedure to compute and

[R] R and Data Storage

2005-09-30 Thread rab45+
Where I work a lot of people end up using Excel spreadsheets for storing data. This has limitations and maybe some less than obvious problems. I'd like to recommend a uniform way for storing and archiving data collected in the department. Most of the data could be stored in simple csv type files