Re: Re: [R] what does the sum of square of Gaussian RVs with differentvariance obey?

2003-09-23 Thread Jean Sun
Thanks a lot. It does work. The fitted data match the simulated ones well. Even no need the shifted or scaled version of Chi-squared pdf. Also, I have tested the case of non-independent RVs,generated by linear combining of independent Gaussian RVs,the result is satisfactory too. Regards, J.Sun

RE: [R] data.frame with duplicated id's

2003-09-23 Thread Andrew Hayen
Try ?reshape A -Original Message- From: Christian Schulz [mailto:[EMAIL PROTECTED] Sent: Wednesday, 24 September 2003 3:42 PM To: [EMAIL PROTECTED] Subject: [R] data.frame with duplicated id's Hi, is there a exstisting function (..i found nothing until now.) what makes it possible tr

[R] data.frame with duplicated id's

2003-09-23 Thread Christian Schulz
Hi, is there a exstisting function (..i found nothing until now.) what makes it possible transfrom a dataset: ID AGE V.MAI V.JUNE 11 20 100 120 12 30 200 90 into IDAGEV 1120 100 1120 120 1230200 123090 ,or have i to programm th

Re: [R] R Production Performance

2003-09-23 Thread Joe Conway
Paul Meagher wrote: Below is the test I ran awhile back on invoking R as a system call. It might be faster if you had a c-extension to R but before I went that route I would want to know 1) roughly how fast Python and Perl are in returning results with their c-bindings/embedded stuff/dcom stuff, 2

Re: [R] Omitting blank lines with read.table

2003-09-23 Thread kjetil brinchmann halvorsen
On 24 Sep 2003 at 11:32, Patrick Connolly wrote: read.table has an argument blank.lines.skip, which is true by default, at least in read.table. You can try to give that to read.delim? Kjetil Halvorsen > Say we have a tab delimited file called bug.txt > > Part Rep CageHb pupae > 1

[R] weird behaviour when calling c++

2003-09-23 Thread Francisco J Molina
I am using 1.7.1 in a PC ( redhat 9, linux ) I created a subroutine in C++, mySubrutine, to be used in R. To debug this subroutine I have a main routine in C++ that calls mySubrutine. The only thing main () does is to provide mySubrutine with its arguments (this is the easiest way for me to deb

Re: [R] confusion about what to expect?

2003-09-23 Thread A.J. Rossini
Marc Schwartz <[EMAIL PROTECTED]> writes: > On Tue, 2003-09-23 at 18:08, A.J. Rossini wrote: <-- about confusion --> Anyway, the mnemonic to use seems to be to remember that complex data structures get simplified whenever possible. Thanks to Doug G. and Patrick C. for that! I'll plead an exc

Re: [R] R Production Performance

2003-09-23 Thread Paul Meagher
Hi Zitan, Below is the test I ran awhile back on invoking R as a system call. It might be faster if you had a c-extension to R but before I went that route I would want to know 1) roughly how fast Python and Perl are in returning results with their c-bindings/embedded stuff/dcom stuff, 2) whether

Re: [R] install problem with R Windows

2003-09-23 Thread Duncan Murdoch
On Tue, 23 Sep 2003 16:16:38 -0500, you wrote: >Dear R People: > >I'm trying to install R 1.7.1 for Windows from Source. > >The error that I get is: >previous declarion of 'ssize_t' >MAKE[2]: ***[internet.o]Error 1 >MAKE[1]: ***[all]Error 1 >MAKE: *** [rmodules] Error 2 > >Any ideas on how to proc

Re: [R] least squares regression line

2003-09-23 Thread Jason Turner
On Wed, 2003-09-24 at 11:04, Carmen Fridell wrote: > I can't seem to find the command to find the least squares regression line > for my bivariate data set. Can you please help? ~Carmen Any time you're lost in R, you can type "help.start()". This will start a web browser, which loads the startin

Re: [R] confusion about what to expect?

2003-09-23 Thread Marc Schwartz
On Tue, 2003-09-23 at 18:08, A.J. Rossini wrote: > In playing around with data.frames (and wanting a simple, cheap way to > use the variable and case names in plots; but I've solved that with > some hacks, yech), I noticed the following behavior with subsetting. > > > testdata <- data.frame(matr

Re: [R] confusion about what to expect?

2003-09-23 Thread Tony Plate
Have you investigated the drop= argument to "["? (as in the expression testdata[,2,drop=F], which will return a dataframe). "[.data.frame" has somewhat different behavior from "[" on matrices with respect to the drop argument: If the result would be a dataframe with a single column, the default

[R] Updating least squares

2003-09-23 Thread David Duffy
Paul Meagher <[EMAIL PROTECTED]> wrote: > > I am looking at developing a user modelling type app (new data points coming > in and wanting to dynamically update regression co-efficients for each user) > which could be viewed as a type of control problem. > Alan Miller's AS274 (in C, f77 or f90) doe

[R] Omitting blank lines with read.table

2003-09-23 Thread Patrick Connolly
Say we have a tab delimited file called bug.txt PartRep CageHb pupae 1 1 S 32 1 M 34 L 42 2 S 36 M 28 L 36

[R] confusion about what to expect?

2003-09-23 Thread A.J. Rossini
In playing around with data.frames (and wanting a simple, cheap way to use the variable and case names in plots; but I've solved that with some hacks, yech), I noticed the following behavior with subsetting. testdata <- data.frame(matrix(1:20,nrow=4,ncol=5)) names(testdata) ## expect labels, ge

[R] least squares regression line

2003-09-23 Thread Carmen Fridell
I can't seem to find the command to find the least squares regression line for my bivariate data set. Can you please help? ~Carmen __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

Re: [R] Typical R installation problem

2003-09-23 Thread Marc Schwartz
On Tue, 2003-09-23 at 15:28, [EMAIL PROTECTED] wrote: > Dear Peter, > > I don't know if this is the proper way to ask for help installing R, but if > not, I presume you can pass this on to the appropriate place. > > I'm trying to install the latest binary on Redhat 9, and keep gett

Re: [R] R Production Performance

2003-09-23 Thread Zitan Broth
Hi James, Thanks for your response :-) - Original Message - > It is like anything else that you want to run as part of web services: what > do you want it to do? Yes, it is fast in doing computations, but what will > you have it do? It is probably as fast as anything else that you will

[R] install R for Windows problem

2003-09-23 Thread Erin Hodgess
Here are all of the messages: In file included from internet.c:858: sock.h.27:conflicting types for 'ssize_t' c:/mingw/include/sys/types.h:119: previous declaration of 'ssize_t' MAKE[2]: ***[internet.o] Error 1 MAKE[1]: ***[all] Error 1 MAKE: ***[rmodules] Error 2 ___

Re: [R] loops in Sweave

2003-09-23 Thread Jason Turner
On Wed, 2003-09-24 at 03:06, Millo Giovanni wrote: > Dear all, > > I was wondering whether there is a way to make loops in Sweave, i.e. for example to: > 1) calculate a parameter, say, a=length(b) > 2) according to that, add #a# chapters to the document, each including some > repetitive analysis,

[R] install problem with R Windows

2003-09-23 Thread Erin Hodgess
Dear R People: I'm trying to install R 1.7.1 for Windows from Source. The error that I get is: previous declarion of 'ssize_t' MAKE[2]: ***[internet.o]Error 1 MAKE[1]: ***[all]Error 1 MAKE: *** [rmodules] Error 2 Any ideas on how to proceed, please? thanks in advance! Sincerely, Erin Hodgess m

[R] Typical R installation problem

2003-09-23 Thread mahustonor
Dear Peter, I don't know if this is the proper way to ask for help installing R, but if not, I presume you can pass this on to the appropriate place. I'm trying to install the latest binary on Redhat 9, and keep getting the same error message no matter what I try. [EMAIL PRO

[R] bug in stack?

2003-09-23 Thread apjaworski
I am posting it here because I am not sure if the behavior described below is actually a bug. If I do something like this: > x1 <- 1:3 > x2 <- 5:9 > x3 <- 21:27 > ll <- list(x1, x2, x3) > stack(x1,x2,x3) I get the following error: Error in rep.int(names(x), lapply(x, length)) : invalid

Re: [R] Plotting multiple lines

2003-09-23 Thread Jonathan Baron
On 09/23/03 15:49, K Skanes wrote: >Hi, > >I have a data set with 7 years worth of data, and 2 different values (a real >value and a model value) of interest in each year for different lengths. I >would like a plot with the year on the y axis and an animal length along the >x axis. For each ye

[R] Plotting multiple lines

2003-09-23 Thread K Skanes
Hi, I have a data set with 7 years worth of data, and 2 different values (a real value and a model value) of interest in each year for different lengths. I would like a plot with the year on the y axis and an animal length along the x axis. For each year I would like to see two lines, one rep

Re: [R] How to extract data from Excel

2003-09-23 Thread Ben Bolker
Your best bet is saving as a comma-separated value file (.csv) and using read.csv to get the data into R. Ben On Tue, 23 Sep 2003 [EMAIL PROTECTED] wrote: > Hi, > > I would like to know how to extract the data from Excel Spreadsheet. > > Thank you very much. > > Melissa > > >

Re: AW: [R] Rank and extract data from a series

2003-09-23 Thread Tony Plate
Using Thomas Unternährer's handy example, one could also do: > X <- c(1, 4.5, 2.3, 1, 7.3) > mean(order(X, decreasing=TRUE)[1:2]) [1] 3.5 > I think this will give the same results as Thomas Unternährer's suggested code in almost all cases, but it is perhaps more concise and direct (provided that

Re: [R] discretization method

2003-09-23 Thread Edgar Acuna
there are plenty of discretization methods, which one are you looking for? (In Spanish hay bastantes metodos de discretizacion cual de ellos estas busacando? Regards (saludos) Edgar Acuna On Tue, 23 Sep 2003, Jaime Lopez Carvajal wrote: > Hi R users > > I need to apply discretization to my conti

Re: [R] How to extract data from Excel

2003-09-23 Thread Uwe Ligges
[EMAIL PROTECTED] wrote: Hi, I would like to know how to extract the data from Excel Spreadsheet. I would like to know whether you have read the manual "R Data Import/Export" before having posted the question. It tells about you more than one way. Uwe Ligges Thank you very much. Melissa ___

[R] How to extract data from Excel

2003-09-23 Thread Melissa_Kuang
Hi, I would like to know how to extract the data from Excel Spreadsheet. Thank you very much. Melissa JLT Risk Solutions Ltd 6 Crutched Friars, London EC3N 2PH. Co Reg No 1536540 Tel: (44) (0)20 7528 4000 Fax: (44) (0)20 7528 4500

Re: [R] filled.contour without box

2003-09-23 Thread Martin Maechler
> "UweL" == Uwe Ligges <[EMAIL PROTECTED]> > on Tue, 23 Sep 2003 11:30:02 +0200 writes: UweL> Jan Kleinn wrote: >> Dear all, >> >> I would like to make a filled contour plot without the >> box R is generating by default around the plotting area, >> i.e. I'm loo

Re: [R] (Fwd) Re: goodfit macro

2003-09-23 Thread Achim Zeileis
On Tuesday 23 September 2003 16:38, Michael Shott wrote: > Dear R-Help: > > As you can see, Prof. Friendly refers me to your site for an > executable version of vcd. I don't mean to be obtuse, but 15 > minutes spent exploring your site failed to locate a downloadable > version of the vcd package

RE: [R] what does the sum of square of Gaussian RVs with differen t variance obey?

2003-09-23 Thread RBaskin
This is a relatively recent article that is somewhat accessible. Jensen, D. R., and Solomon, Herbert (1994), "Approximations to joint distributions of definite quadratic forms", Journal of the American Statistical Association, 89 , 480-486 It has references to previous work. I also have an old pa

AW: [R] weighted standard deviation

2003-09-23 Thread Schnitzler, Johannes
Thank you all for the reply, the Hmisc library is exactly what i was looking for. Johannes __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

Re: [R] (Fwd) Re: goodfit macro

2003-09-23 Thread Uwe Ligges
Michael Shott wrote: Dear R-Help: As you can see, Prof. Friendly refers me to your site for an executable version of vcd. I don't mean to be obtuse, but 15 minutes spent exploring your site failed to locate a downloadable version of the vcd package to which he referred. I know plainly what t

[R] loops in Sweave

2003-09-23 Thread Millo Giovanni
Dear all, I was wondering whether there is a way to make loops in Sweave, i.e. for example to: 1) calculate a parameter, say, a=length(b) 2) according to that, add #a# chapters to the document, each including some repetitive analysis, each time done on a particular subset of the data indexed by t

Re: [R] discretization method

2003-09-23 Thread Ben Bolker
Consider: ?cut ?round With more detail we might be able to help more ... On Tue, 23 Sep 2003, Jaime Lopez Carvajal wrote: > Hi R users > > I need to apply discretization to my continuous data. > Is there a method in R to do this? > > Thanks in advance, > > Jaime > -- 620B Bartram Ha

Re: [R] discretization method

2003-09-23 Thread Thomas W Blackwell
On Tue, 23 Sep 2003, Jaime Lopez Carvajal wrote: > I need to apply discretization to my continuous data. > Is there a method in R to do this? See help("cut"). - tom blackwell - u michigan medical school - ann arbor - __ [EMAIL PROTECTED] maili

[R] (Fwd) Re: goodfit macro

2003-09-23 Thread Michael Shott
Dear R-Help: As you can see, Prof. Friendly refers me to your site for an executable version of vcd. I don't mean to be obtuse, but 15 minutes spent exploring your site failed to locate a downloadable version of the vcd package to which he referred. I know plainly what this application can do

[R] discretization method

2003-09-23 Thread Jaime Lopez Carvajal
Hi R users I need to apply discretization to my continuous data. Is there a method in R to do this? Thanks in advance, Jaime __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

Re: [R] problems installing Design and Hmisc libs

2003-09-23 Thread Thomas Lumley
On Tue, 23 Sep 2003, Federico Calboli wrote: > Dear All, > > when I try to: > > install.packages("Design"); install.packages("Hmisc") > > I get the following error messages: > > * Installing *source* package 'Design' ... > ** libs > g77 -mieee-fp -O2 -fomit-frame-pointer -pipe -march=i586 > -mcpu

Re: [R] what does the sum of square of Gaussian RVs with different variance obey?

2003-09-23 Thread Thomas Lumley
On Tue, 23 Sep 2003, Jean Sun wrote: > >From basic statistics principle,we know,given several i.i.d Gaussian > >RVs with zero or nonzero mean,the sum of square of them is a central or > >noncentral Chi-distributed RV.However if these Gaussian RVs have > >different variances,what does the sum of sq

Re: [R] filled.contour without box

2003-09-23 Thread Roger D. Peng
If you just want to get rid of the axes, you can do filled.contour(x, plot.axes = { }) -roger Uwe Ligges wrote: Jan Kleinn wrote: Dear all, I would like to make a filled contour plot without the box R is generating by default around the plotting area, i.e. I'm looking for an option in fille

[R] Marginal Means with the lme()

2003-09-23 Thread [EMAIL PROTECTED]
I need help computing model-predicted estimated marginal means of the dependent variable in the cells defined by the fixed factors, including interactions. By estimated marginal means, I just mean cell means collapsed over all other factors and adjusted for the mean values of other covariates i

AW: [R] Rank and extract data from a series

2003-09-23 Thread "Unternährer Thomas, uth"
Hi, >I would like to rank a time-series of data, extract the top ten data items from this >series, determine the >corresponding row numbers for each value in the sample, and take a mean of these *row >numbers* (not the data). >I would like to do this in R, rather than pre-process the data o

RE: [R] Rank and extract data from a series

2003-09-23 Thread Liaw, Andy
Here's one way. Suppose your "time series" is in a vector called "x". top10 <- sort(x, decreasing=TRUE)[1:10] mean.index <- mean(which(x %in% top10)) HTH, Andy > -Original Message- > From: James Brown [mailto:[EMAIL PROTECTED] > Sent: Tuesday, September 23, 2003 7:51 AM > To: [EMAIL PR

RE: [R] R-project [.com?] [.net?]

2003-09-23 Thread Liaw, Andy
> From: Murray Jorgensen [mailto:[EMAIL PROTECTED] > > I got a shock a few days ago when I accidentally visited > www.r-project.com . I thought that the r-project site had been hacked This one seems to be about some sort of city revival projects in Japan. (The introduction starts with "Recycle

[R] Rank and extract data from a series

2003-09-23 Thread James Brown
I would like to rank a time-series of data, extract the top ten data items from this series, determine the corresponding row numbers for each value in the sample, and take a mean of these *row numbers* (not the data). I would like to do this in R, rather than pre-process the data on the UNIX comm

Re: [R] problems installing Design and Hmisc libs

2003-09-23 Thread Federico Calboli
> You need the Fortran compiler (g77) as well (as indicated by the error > message). Is it installed? Is it in your path? > > Uwe Ligges > It was not installed. Now everything works. Cheers, Federico -- = Federico C. F. Calboli Department of Biology Univ

RE: [R] Date on x-axis of xyplot

2003-09-23 Thread Heywood, Giles
One thing you could do is to use the 'its' (irregular time-series) package on CRAN. e.g. using a trivial dataset require(its) its.format("%Y-%m-%d") #defines text format of dates in dimnames df <- data.frame(1:3,(1:3)^2) dimnames(df) <- list(c("2003-01-03","2003-01-06","2003-01-07"),letters[1:2]

Re: [R] problems installing Design and Hmisc libs

2003-09-23 Thread Uwe Ligges
Federico Calboli wrote: Dear All, when I try to: install.packages("Design"); install.packages("Hmisc") I get the following error messages: * Installing *source* package 'Design' ... ** libs g77 -mieee-fp -O2 -fomit-frame-pointer -pipe -march=i586 -mcpu=pentiumpro -O2 -fomit-frame-pointer -p

[R] problems installing Design and Hmisc libs

2003-09-23 Thread Federico Calboli
Dear All, when I try to: install.packages("Design"); install.packages("Hmisc") I get the following error messages: * Installing *source* package 'Design' ... ** libs g77 -mieee-fp -O2 -fomit-frame-pointer -pipe -march=i586 -mcpu=pentiumpro -O2 -fomit-frame-pointer -pipe -march=i586 -mcpu=pen

Re: [R] Plotting of the lm

2003-09-23 Thread Gavin Simpson
Try this: plot(x, y) abline(lm(y ~ x)) #plots the fitted line or this if you need the model results elsewhere: mod.lm <- lm(y ~ x) # store the model plot(x, y)# plot the data abline(mod.lm)# plot the fitted line see ?plot.lm and ?abline HTH Gav [EMAIL PROTECTED] wrote: H

Re: [R] filled.contour without box

2003-09-23 Thread Uwe Ligges
Jan Kleinn wrote: Dear all, I would like to make a filled contour plot without the box R is generating by default around the plotting area, i.e. I'm looking for an option in filled.contour similar to 'axes=F' in 'contour' or in 'plot'. I couldn't find any option to get rid of the box, any help

Re: [R] Plotting of the lm

2003-09-23 Thread Uwe Ligges
[EMAIL PROTECTED] wrote: Hi, I would like to enquire if by typing plot (lm(y~x)) would this show me the plot of the fitted line? I tried this function previously but I was only able to get the last 4 plots starting with Residuals vs fitted. No, it shows plots for analyses of the residuals. Try

[R] filled.contour without box

2003-09-23 Thread Jan Kleinn
Dear all, I would like to make a filled contour plot without the box R is generating by default around the plotting area, i.e. I'm looking for an option in filled.contour similar to 'axes=F' in 'contour' or in 'plot'. I couldn't find any option to get rid of the box, any help is welcome. Thank

[R] Plotting of the lm

2003-09-23 Thread Melissa_Kuang
Hi, I would like to enquire if by typing plot (lm(y~x)) would this show me the plot of the fitted line? I tried this function previously but I was only able to get the last 4 plots starting with Residuals vs fitted. Thank You. Melissa **

Re: [R] R Production Performance

2003-09-23 Thread Laurent Faisnel
Zitan Broth wrote: Greetings All, Been playing with R and it is very easy to get going with the UI or infile batch commands :-) What I am wondering is how scalable and fast R is for running as part of a web service. I believe R is written in C which is a great start, but what are peoples general

Re: [R] Very small estimated random effect variance (lme)

2003-09-23 Thread Peter Dalgaard BSA
"Remko Duursma" <[EMAIL PROTECTED]> writes: > Dear R-helpers, > > i get some strange results using a linear mixed-effects model (lme), of the type: > > lme1 <- lme(y ~ x, random=~x|group, ...) > > For some datasets, i obtain very small standard deviations of the random effects. I > compared th