[R] An entire data frame which is a time-series?

2004-08-17 Thread Ajay Shah
I have : raw - read.table(monthly.text, skip=3, sep=|, col.names=c(junk, junk2, wpi, g.wpi, wpi.primary, g.wpi.primary, wpi.fuel, g.wpi.fuel, wpi.manuf, g.wpi.manuf, cpi.iw, g.cpi.iw, cpi.unme, g.cpi.unme,

[R] Cross-variograms

2004-08-17 Thread Edzer J. Pebesma
Jacques, provided that X and Y are colocated (i.e., have exactly the same observation locations), you get the cross variogram right; the definition of this cross variogram is however: gamma(h)= E[(X(s)-X(s+h))*(Y(s)-Y(s+h))] also, where you select: cv - v$gamma[1:14] you may be better off using

[R] strptime() bug? And additional problem in package tseries

2004-08-17 Thread javier garcia - CEBAS
Hi all, I've got some problems with irts objects, one of which could be a bug: 1) Read a table with several columns from Postgres and the first column is Timestamp with timezone (this is OK). An extract is: raincida$ts: [2039] 25/03/2000 22:00:00 UTC 25/03/2000 23:00:00 UTC [2041] 26/03/2000

NO bug in Re: [R] strptime() bug? And additional problem in package tseries

2004-08-17 Thread Prof Brian Ripley
There is no bug in R here. There was a change to DST in Spain at 2am on 2000-03-26, and they are *printed* as times in your locale, as documented. Please read the posting guide and FAQ about what is a bug. Also, please try not to confuse an object and its printed representation. On Tue, 17 Aug

[R] table and getting rownames

2004-08-17 Thread merser
hi there say that i have this table x-table(adoc, oarb) x oarb 0 1 adoc ab1 0 am5 1 ba 14 1 cc 271 3 ch 87 2 dz 362 6 fl

Re: [R] table and getting rownames

2004-08-17 Thread Prof Brian Ripley
On Tue, 17 Aug 2004 [EMAIL PROTECTED] wrote: say that i have this table x-table(adoc, oarb) x oarb 0 1 adoc ab1 0 am5 1 ba 14 1 cc 271 3 ch 87 2

[R] Re: Thanks Frank, setting graph parameters, and why social scientists don't use R

2004-08-17 Thread david_foreman
First, many thanks to Frank Harrell for once again helping me out. This actually relates to the next point, which is my contribution to the 'why don't social scientists use R' discussion. I am a hybrid social scientist(child psychiatrist) who trained on SPSS. Many of my difficulties in

Re: [R] Re: Thanks Frank, setting graph parameters, and why social scientists don't use R

2004-08-17 Thread Roger D. Peng
I'm just curious, but how do social scientists, or anyone else for that matter, learn SPSS, besides taking a class? -roger [EMAIL PROTECTED] wrote: First, many thanks to Frank Harrell for once again helping me out. This actually relates to the next point, which is my contribution to the 'why

Re: [R] using nls to fit a four parameter logistic model

2004-08-17 Thread Rogers, James A [PGRD Groton]
Shalini, I think your hill equation is meant to just be an alternative parameterization of the four parameter logistic (BTW, the hill *coefficient* is a function of the slope parameter of the FPL, but I don't believe hill equation is standard terminology). Note conc is the input in this

[R] Bug in colnames of data.frames?

2004-08-17 Thread Arne Henningsen
Hi, I am using R 1.9.1 on on i686 PC with SuSE Linux 9.0. I have a data.frame, e.g.: myData - data.frame( var1 = c( 1:4 ), var2 = c (5:8 ) ) If I add a new column by myData$var3 - myData[ , var1 ] + myData[ , var2 ] everything is fine, but if I omit the commas: myData$var4 - myData[ var1

Re: [R] Bug in colnames of data.frames?

2004-08-17 Thread Uwe Ligges
Arne Henningsen wrote: Hi, I am using R 1.9.1 on on i686 PC with SuSE Linux 9.0. I have a data.frame, e.g.: myData - data.frame( var1 = c( 1:4 ), var2 = c (5:8 ) ) If I add a new column by myData$var3 - myData[ , var1 ] + myData[ , var2 ] everything is fine, but if I omit the commas:

[R] Fwd: strptime() problem?

2004-08-17 Thread javier garcia - CEBAS
Hi all; I've already send a similar e-mail to the list and Prof. Brian Ripley answered me but my doubts remain unresolved. Thanks for the clarification, but perhaps I wasn't clear enough in posting my questions. I've got a postgres database which I read into R. The first column is Timestamp

Re: [R] Bug in colnames of data.frames?

2004-08-17 Thread Peter Dalgaard
Arne Henningsen [EMAIL PROTECTED] writes: Hi, I am using R 1.9.1 on on i686 PC with SuSE Linux 9.0. I have a data.frame, e.g.: myData - data.frame( var1 = c( 1:4 ), var2 = c (5:8 ) ) If I add a new column by myData$var3 - myData[ , var1 ] + myData[ , var2 ] everything is

Re: [R] Bug in colnames of data.frames? -- NOT

2004-08-17 Thread Prof Brian Ripley
This is not a bug, and BTW data frames have names not colnames. As I have said already today, don't confuse the printed repesentation of an object with the object itself. On Tue, 17 Aug 2004, Arne Henningsen wrote: I am using R 1.9.1 on on i686 PC with SuSE Linux 9.0. I have a data.frame,

[R] levels of factor

2004-08-17 Thread Luis Rideau Cruz
R-help, I have a data frame wich I subset like : a - subset(df,df$column2 %in% c(factor1,factor2) df$column2==1) But when I type levels(a$column2) I still get the same levels as in df (my original data frame) Why is that? Is it right? Luis Luis Ridao Cruz Fiskirannsóknarstovan Nóatún 1

Re: [R] using nls to fit a four parameter logistic model

2004-08-17 Thread Spencer Graves
In your second model, log(conc) is negative for conc = 0.001. This observation will generate NA for (log(conc)/xmid)^scal unless scal is an integer or xmid is also negative. In the latter case, (log(conc)/xmid)^scal will be NA for all but that last value unless scal is an integer.

Re: [R] Bug in colnames of data.frames?

2004-08-17 Thread Marc Schwartz
On Tue, 2004-08-17 at 09:01, Arne Henningsen wrote: Hi, I am using R 1.9.1 on on i686 PC with SuSE Linux 9.0. I have a data.frame, e.g.: myData - data.frame( var1 = c( 1:4 ), var2 = c (5:8 ) ) If I add a new column by myData$var3 - myData[ , var1 ] + myData[ , var2 ]

Re: [R] Bug in colnames of data.frames?

2004-08-17 Thread Marc Schwartz
On Tue, 2004-08-17 at 09:34, Marc Schwartz wrote: Take a look at the details, value and coercion sections of ?.data.frame This must be my week for typos. That should be: ?[.data.frame (in ESS) or ?[.data.frame (otherwise) Marc __ [EMAIL

Re: [R] levels of factor

2004-08-17 Thread Marc Schwartz
On Tue, 2004-08-17 at 09:30, Luis Rideau Cruz wrote: R-help, I have a data frame wich I subset like : a - subset(df,df$column2 %in% c(factor1,factor2) df$column2==1) But when I type levels(a$column2) I still get the same levels as in df (my original data frame) Why is that? The

Re: [R] Fwd: strptime() problem?

2004-08-17 Thread Gabor Grothendieck
javier garcia - CEBAS rn001 at cebas.csic.es writes: : : Hi all; : I've already send a similar e-mail to the list and Prof. Brian Ripley : answered me but my doubts remain unresolved. Thanks for the clarification, : but perhaps I wasn't clear enough in posting my questions. : : I've got a

RE: [R] Fwd: strptime() problem?

2004-08-17 Thread Whit Armstrong
Javier, I recently had a problem with dates. This example might shed some light on your problem. x - ISOdate(rep(2000,2),rep(3,2),rep(26,2),hour=0) x [1] 2000-03-26 GMT 2000-03-26 GMT unclass(x) [1] 954028800 954028800 attr(,tzone) [1] GMT When one creates a date with ISOdate, the

Re: [R] Re: Thanks Frank, setting graph parameters,and why social scientists don't use R

2004-08-17 Thread Berton Gunter
A few comments: First, your remarks are interesting and, I would say, mainly well founded. However, I think they are in many respects irrelevant, although they do point to the much bigger underlying issue, which Roger Peng also hinted at in his reply. I think they are sensible because R IS

Re: [R] An entire data frame which is a time-series?

2004-08-17 Thread Gabor Grothendieck
Ajay Shah ajayshah at mayin.org writes: : : I have : : : raw - read.table(monthly.text, skip=3, sep=|, : col.names=c(junk, junk2, : wpi, g.wpi, wpi.primary, g.wpi.primary, : wpi.fuel, g.wpi.fuel, wpi.manuf, g.wpi.manuf, :

[R] aov summary to matrix

2004-08-17 Thread Moises Hassan
Is there an easy way of converting an aov.summary into a matrix in which the rows are the factor names and the columns are Df, Sum Sq, Mean Sq, F value and Pr. For example, convert Df Sum Sq Mean Sq F value Pr(F) block5 343.29 68.66 4.4467 0.015939 * N1

Re: [R] aov summary to matrix

2004-08-17 Thread Prof Brian Ripley
On Tue, 17 Aug 2004, Moises Hassan wrote: Is there an easy way of converting an aov.summary into a matrix in which the rows are the factor names and the columns are Df, Sum Sq, Mean Sq, F value and Pr. You are confusing the printed representation with the object (which seems today's

Re: [R] survdiff

2004-08-17 Thread Peter Dalgaard
Krista Haanstra [EMAIL PROTECTED] writes: As I am quitte an ignorant user of R, excuse me for any wrongfull usage of all the terms. My question relates to the statistics behind the survdiff function in the package survival. My textbook knowledge of the logrank test tells me that if I want to

Re: [R] aov summary to matrix

2004-08-17 Thread Sundar Dorai-Raj
Moises Hassan wrote: Is there an easy way of converting an aov.summary into a matrix in which the rows are the factor names and the columns are Df, Sum Sq, Mean Sq, F value and Pr. For example, convert Df Sum Sq Mean Sq F value Pr(F) block5 343.29 68.66 4.4467

Re: [R] Bug in colnames of data.frames?

2004-08-17 Thread David Forrest
On Tue, 17 Aug 2004, Arne Henningsen wrote: Thank you for all your answers! I agree with you that it is not a bug. My mistake was that I thought that a data frame is similar to a matrix, but as ?data.frame says they ... share many of the properties of matrices and of lists. ... I think the

Re: [R] Re: Thanks Frank, setting graph parameters, and why social scientists don't use R

2004-08-17 Thread John
On Tuesday 17 August 2004 06:14, Roger D. Peng wrote: I'm just curious, but how do social scientists, or anyone else for that matter, learn SPSS, besides taking a class? They sit down with a book, a computer, and data they desperately need to analyze and start working. SPSS documentation and

Re: [R] Re: Thanks Frank, setting graph parameters, and why social scientists don't use R

2004-08-17 Thread John
On Tuesday 17 August 2004 09:20, Berton Gunter wrote: A few comments: It has been decades since I used SPSS. At that time, to really work with it you edited a text file program that identified the data file and variable columns you wanted to work with. You assembled the flow of work commands

[R] logistic -normal model

2004-08-17 Thread syl dan
I am working with a logistic-normal model (i.e, GLMM with random intercept model) by Bayesian method. BUt I met some difficulities for programming by R. Is there anyone have experience of this model or the R code I can refer as example? Thanks for your help. Syl

[R] creating a plot

2004-08-17 Thread Konrad Banachewicz
Hi, I have a time series plot to produce, yet I want the x-axis to be labelled with dates (stored on another array) and not with observation numbers. Can anyone suggest me how? Thanks. Konrad

Re: [R] creating a plot

2004-08-17 Thread Tomas Aragon
--- Konrad Banachewicz [EMAIL PROTECTED] wrote: Hi, I have a time series plot to produce, yet I want the x-axis to be labelled with dates (stored on another array) and not with observation numbers. Can anyone suggest me how? Thanks.