[R] [R-pkgs] New package aod: Analysis of Overdispersed Data

2005-06-10 Thread Renaud Lancelot
Information on package 'aod' Description: Package: aod Version: 1.1-2 Date: 2005-06-08 Title: Analysis of Overdispersed Data Author:Matthieu Lesnoff [EMAIL PROTECTED] and Renaud Lancelot [EMAIL PROTECTED] Maintainer:Renaud

Re: [R] Plot/manage spatial boundary data

2005-06-10 Thread Roger Bivand
On Thu, 9 Jun 2005, David Forrest wrote: I have some disconnected boundary data from a finite element ocean model and I'd like to make a plot. Maptools looks promising, but since my data is not in a shapefile or a map, I'm unclear on what the best way to approach the problem. If the line

[R] Contour Plot

2005-06-10 Thread WeiQiang . Li
Hi ALL, How can I differentiate the area that does not have data or the data is below the minimum value using color when I plot Contour chart. For example, my data value range is from 5.0 to 7.5, but I am interested in the data value from 5.5 to 7.5, so I change levels in the

[R] How to read a row dataset one by one

2005-06-10 Thread Jan Sabee
Dear all, How to read a row dataset one by one and then print it. x1 x2 x3 x4 x5 y a b a c cM1 c b b c cM4 c c a c cM2 c a c a aM2 c c a a aM1 c a b c aM3 c c a b cM3 c a c a bM2 c c a b aM1 I need a result like read row

Re: [R] Robustness of Segmented Regression Contributed by Muggeo

2005-06-10 Thread vito muggeo
Hi, sorry for my delay.. In addition to valuable Achim's comments. As Achim said, you can try different starting values to assess how the final solution depends on them. Then select one having the best logLik (or the minimum RSS). Everybody dealing with nonlinear models knows that the logLik

[R] Replacing for loop with tapply!?

2005-06-10 Thread Sander Oom
Dear all, We have a large data set with temperature data for weather stations across the globe (15000 stations). For each station, we need to calculate the number of days a certain temperature is exceeded. So far we used the following S code, where mat88 is a matrix containing rows of 365

[R] R installation

2005-06-10 Thread Léo Tranchevent
Hi everybody, I am trying to update R on my computer, and I got the following problem: After tar the downloaded file, when I try ./configure, I got the following error configure: error: --with-readline=yes (default) and headers/libs are not available So I try ./configure

Re: [R] How to read a row dataset one by one

2005-06-10 Thread Clark Allan
use a loop associated with the scan function. for (i in 1:9) { print(scan(file=c:/a.txt,sep=\t,skip=i,nlines=1,fill=T,quiet=T,what=raw)) } this works but there has to be a better solution Jan Sabee wrote: Dear all, How to read a row dataset one by one and then print it. x1 x2 x3 x4

Re: [R] Replacing for loop with tapply!?

2005-06-10 Thread Dimitris Rizopoulos
maybe you are looking for something along these lines: mat - matrix(sample(-15:50, 365 * 15000, TRUE), 365, 15000) temps - c(37, 39, 41) # ind - matrix(0, length(temps), ncol(mat)) for(i in seq(along = temps)) ind[i, ] - colSums(mat temps[i]) ind I hope it helps. Best,

Re: [R] dir() and RegEx and gsub()

2005-06-10 Thread Hans-Peter
2005/6/9, Sarah Goslee [EMAIL PROTECTED]: Under R, for reasons I've never quite understood, \\. evaluates to . Thanks to the answers of B. Ripley and Gabor I think, I understand now: 1) the patter-string \\.csv$ gives the regular expression \.csv$ 2) now the backslash lets the dot to be

Re: [R] How to read a row dataset one by one

2005-06-10 Thread Henrik Bengtsson
Open a connection a read line by line from that one, e.g. myReadPrint - function(pathname, ...) { con - file(pathname, open=r) on.exit(close(con)) # Guarantees to close connection! count - 0; while(TRUE) { line - scan(con, sep=\t, nlines=1, fill=TRUE,

Re: [R] Replacing for loop with tapply!?

2005-06-10 Thread Sander Oom
Thanks Dimitris, Very impressive! Much faster than before. Thanks to new found R.basic, I can simply rotate the result with rotate270{R.basic}: mat - matrix(sample(-15:50, 365 * 15000, TRUE), 365, 15000) temps - c(37, 39, 41) # #ind - matrix(0, length(temps), ncol(mat))

[R] Different results in different runs with identical parameters in CLARA

2005-06-10 Thread Amir Safari
Dear All R Friends, When I run my data in any time with the below codes, I receive different results. My data , k , samples, trace are identical in any run. c- clara(mydata,4, metric= euclidean , stand= TRUE, samples=5 , trace=3, keep.data=TRUE , rngR=TRUE) result of first try:

Re: [R] R installation

2005-06-10 Thread Uwe Ligges
Lo Tranchevent wrote: Hi everybody, I am trying to update R on my computer, and I got the following problem: After tar the downloaded file, when I try ./configure, I got the following error configure: error: --with-readline=yes (default) and headers/libs are not available

Re: [R] How to read a row dataset one by one

2005-06-10 Thread Jan Sabee
For Henrik and Clark, thanks for your help. Then If I load to dataframe, MM16 - read.table(G:\\Stuff\\data\\MM16.txt) MM16 x1 x2 x3 x4 x5 y a b a c cM1 c b b c cM4 c c a c cM2 c a c a aM2 c c a a aM1 c a b c aM3 c c a b cM3 c a c a b

Re: [R] dir() and RegEx and gsub()

2005-06-10 Thread Uwe Ligges
Hans-Peter wrote: 2005/6/9, Sarah Goslee [EMAIL PROTECTED]: Under R, for reasons I've never quite understood, \\. evaluates to . Thanks to the answers of B. Ripley and Gabor I think, I understand now: 1) the patter-string \\.csv$ gives the regular expression \.csv$ 2) now the

Re: [R] How to read a row dataset one by one

2005-06-10 Thread Henrik Bengtsson
Jan, I'm not sure what you asking for. First, you cannot make read.table() to output verbose message at every line read. Second, if you have trouble to read you file, which looks like what you are showing, then read ?read.table carefully. Most likely you wish to add argument header=TRUE

Re: [R] Rpy and RSPython

2005-06-10 Thread François Pinard
[Weiwei Shi] I am thinking to use one of them but not sure which one is better. I think Rpy cannot call python from R while the PRPython can in two-directional calling. Am I right? s/PRPython/RSPython/ ? :-) This is also what I understood. Yet, despite the uni-directionality of RPy, this

Re: [R] Different results in different runs with identical parametersin CLARA

2005-06-10 Thread TEMPL Matthias
Dear All R Friends, When I run my data in any time with the below codes, I receive different results. Of course. See in L. Kaufman and P. Rousseeuw. Finding Groups in Data. John Wiley Sons, Inc, 1990. There is a random part in clara. My data , k , samples, trace are identical in

[R] Sum up the values of a function

2005-06-10 Thread Stefan Wagner
Dear R-Users, I have to do a maximum-likelihood estimation and have now a problem concerning how to sum up my function values in a smart way. I don't know how to explain it easyly, so I give you the code and show you where my problem is. I have shorten the code a little bit, so that you only

[R] Replies of the question about robustness of segmented regression

2005-06-10 Thread Park, Kyong H Mr. RDECOM
I appreciate to Roger Koenker, Achim Zeileis and Vito Muggeo for their informative answers. Listed below is unedited replies I got followed by the question I posted. Kyong 1. Roger Koenker: You might try rqss() in the quantreg package. It gives piecewise linear fits for a nonparametric form

Re: [R] Sum up the values of a function

2005-06-10 Thread Petr Pikal
Hi Your example is not reproducible as we do not know x1. Do you by chance seek for something like aggregate? If yes see ?aggregate or ?by BTW, do you have some trouble with your space bar? HTH Petr On 10 Jun 2005 at 13:08, Stefan Wagner wrote: Dear R-Users, I have to do a

Re: [R] Sum up the values of a function

2005-06-10 Thread Stefan Wagner
I am really sorry, that the code appears in a bad way but it is really not my fault (at least I think so). I think my mail-programm is responsible for it because I sent it in a normal way without havin g trouble with my space bar. Here are the missing values (but please remember, this is only a

[R] contours and cutlines

2005-06-10 Thread Robin Hankin
Hi Using contour(), one can plot pretty pictures (of potential flow); my problem is removing the cutline. If I do: x - seq(from=-4,to=4,len=100) z - outer(x,1i*x,+) contour(x,x,Im(log(z))) everything is fine except the discontinuity along the negative real axis, which derives from log()'s

[R] data.frame to character

2005-06-10 Thread Muhammad Subianto
Hi, Excuse me for this simple question. How to convert as.data.frame to as.character? ?data.frame L3 - LETTERS[1:3] L10 - LETTERS[1:10] d - data.frame(cbind(x=c(XYZ), y=L10), fac=sample(L3, 10, repl=TRUE)) d x y fac 1 XYZ A A 2 XYZ B A 3 XYZ C A 4 XYZ D A 5 XYZ E

Re: [R] data.frame to character

2005-06-10 Thread Liaw, Andy
Is this what you want? d[] - lapply(d, as.character) str(d) `data.frame': 10 obs. of 3 variables: $ x : chr XYZ XYZ XYZ XYZ ... $ y : chr A B C D ... $ fac: chr B A C B ... Andy From: Muhammad Subianto Hi, Excuse me for this simple question. How to convert as.data.frame to

Re: [R] data.frame to character

2005-06-10 Thread Gabor Grothendieck
On 6/10/05, Muhammad Subianto [EMAIL PROTECTED] wrote: Hi, Excuse me for this simple question. How to convert as.data.frame to as.character? ?data.frame L3 - LETTERS[1:3] L10 - LETTERS[1:10] d - data.frame(cbind(x=c(XYZ), y=L10), fac=sample(L3, 10, repl=TRUE)) d[] -

Re: [R] logistic regressioin - course ornotes

2005-06-10 Thread Wuming Gong
Hi Stephen, I think the tutorial written by Brett Presnell (http://web.stat.ufl.edu/~presnell/Teaching/sta4504-2000sp/R/) is a good start point for categorical data analysis using R. Wuming On 6/10/05, Stephen Choularton [EMAIL PROTECTED] wrote: Hi I am using R for logistic regression and

Re: [R] How to plot more than 3 sets in Venn Diagrams?

2005-06-10 Thread Jim Lemon
Tan Hui Hui Jenny wrote: I'm trying to plot Venn diagrams with more than 3 sets (5 actually) in order to describe graphically the genetic variation between populations. I tried the limma library but realised it can only plot 3 sets. Is there any solution? Of course I could plot the chart

Re: [R] data.frame to character - thanks

2005-06-10 Thread Muhammad Subianto
Dear all, Thank you very much for your help. I would like to thank Andy Liaw and Gabor Grothendieck for their fast help. Regards, Muhammad Subianto On this day 6/10/2005 2:30 PM, Liaw, Andy wrote: Is this what you want? d[] - lapply(d, as.character) str(d) `data.frame': 10 obs. of

[R] Error with function lda in package MASS (dimnames not equal?)

2005-06-10 Thread Joshua Gilbert
This question appears to have been asked previously, but not answered. the last response I can find to this previous thread is here: http://tolstoy.newcastle.edu.au/R/help/04/07/0126.html. The asnwer was to provide debugging info, not an answer. So the problem is that I'm trying to use lda on my

Re: [R] data.frame to character - thanks

2005-06-10 Thread Gabor Grothendieck
Here is one minor improvement -- it does not overwrite the input: replace(d,,as.matrix(d)) On 6/10/05, Muhammad Subianto [EMAIL PROTECTED] wrote: Dear all, Thank you very much for your help. I would like to thank Andy Liaw and Gabor Grothendieck for their fast help. Regards, Muhammad

Re: [R] Error with function lda in package MASS (dimnames not equal?)

2005-06-10 Thread Gabor Grothendieck
On 6/10/05, Joshua Gilbert [EMAIL PROTECTED] wrote: This question appears to have been asked previously, but not answered. the last response I can find to this previous thread is here: http://tolstoy.newcastle.edu.au/R/help/04/07/0126.html. The asnwer was to provide debugging info, not an

Re: [R] Error with function lda in package MASS (dimnames not equal?)

2005-06-10 Thread Prof Brian Ripley
lda.default - MASS:::lda.default and proceed. Look at the group means of your data: they are identical to machine accuracy. The question has to be `why are you trying to use lda to separate two groups with identical means'? Lda is not protected against that and it is rather unlikely unless you

[R] In connection with the creation of a grid of point with gstat

2005-06-10 Thread alain ponsero
Hi ALL, With the pakage gstat How can one create a grid of points as in the meuse.grid example, from measuring sites as in the Meuse example Thanks in advance! Alain PONSERO Conservateur Réserve Naturelle de la Baie de Saint-Brieuc site de l'Etoile 22120 HILLION tel/fax :

[R] rgl.snapshot failed

2005-06-10 Thread roger koenker
I've installed the rgl package on a Suse x86-64 machine (further details below) and it produces nice screen images. Unfortunately, rgl.snapshot attempts to make png files produces only the response failed. For other graphics png() works fine, and capabilities indicates that it is there.

Re: [R] question on boot

2005-06-10 Thread Frank Samuelson
boot will handle multivariate data as well as multistratum data, though it isn't clear from your question what is multivariate about your problem or what you are bootstrapping. Lisa Wang wrote: Hello there: Dear Dr. Murdoch, I'm a statistician at Princess Margaret Hospital. Could you

[R] What is median in survfit

2005-06-10 Thread Remigijus Lapinskas
Dear All, A very simple question: library(survival) fit - coxph(Surv(time, status) ~ x, data=aml) survfit(fit) Call: survfit.coxph(object = fit) n events median 0.95LCL 0.95UCL 23 18 30 18 45 I believe I know what is median here, but how to extract it?

Re: [R] In connection with the creation of a grid of point with gstat

2005-06-10 Thread Roger Bivand
On Fri, 10 Jun 2005, alain ponsero wrote: Hi ALL, With the pakage gstat How can one create a grid of points as in the meuse.grid example, from measuring sites as in the Meuse example In the meuse example, the meuse.grid locations are provided. They are essentially the same

Re: [R] What is median in survfit

2005-06-10 Thread Thomas Lumley
On Fri, 10 Jun 2005, Remigijus Lapinskas wrote: Dear All, A very simple question: library(survival) fit - coxph(Surv(time, status) ~ x, data=aml) survfit(fit) Call: survfit.coxph(object = fit) n events median 0.95LCL 0.95UCL 23 18 30 18 45 I

Re: [R] rgl.snapshot failed

2005-06-10 Thread Duncan Murdoch
On 6/10/2005 10:17 AM, roger koenker wrote: I've installed the rgl package on a Suse x86-64 machine (further details below) and it produces nice screen images. Unfortunately, rgl.snapshot attempts to make png files produces only the response failed. For other graphics png() works

[R] Abrupt shut down of R session

2005-06-10 Thread Park, Kyong H Mr. RDECOM
Dear R users, I'm using R 2.1.0 with Windows 2000. In the middle of a R session, an error message pops up saying Rgui.exe generated errors and shuts down the R session. How can I correct this problem? Appreciate your help. Kyong [[alternative HTML version deleted]]

Re: [R] How to plot more than 3 sets in Venn Diagrams?

2005-06-10 Thread Charles Geyer
Message: 23 Date: Thu, 9 Jun 2005 09:53:02 -0500 From: [EMAIL PROTECTED] Subject: RE: [R] How to plot more than 3 sets in Venn Diagrams? To: Tan Hui Hui Jenny [EMAIL PROTECTED], r-help@stat.math.ethz.ch Message-ID: [EMAIL PROTECTED] Content-Type: text/plain;

Re: [R] Rpy and RSPython

2005-06-10 Thread Weiwei Shi
Hi, I tried to install RSPython, here is my problem: I load the env variables correctly as below: #-- # R section #-- R_HOME=/usr/lib/R PYTHONPATH=${R_HOME}/library/RSPython/Python:${R_HOME}/library/RSPython/libs

[R] Fortran compilation error

2005-06-10 Thread Brian Beckage
Hello, I'm trying to install a package that requires a Fortran compiler (Hmisc) using R CMD INSTALL. I downloaded the package source onto my Desktop, unzipped it, and then typed: R CMD INSTALL /Users/brianbeckage/Desktop/Hmisc * Installing *source* package 'Hmisc' ... ** libs g77

[R] discovery (was: data.frame to character)

2005-06-10 Thread Robert Citek
How can one discover or list all available built-in objects? On Jun 10, 2005, at 7:23 AM, Muhammad Subianto wrote: L3 - LETTERS[1:3] L10 - LETTERS[1:10] LETTERS is apparently a built-in character vector. ls() and objects () only lists the ones I've created. Is there a function that lists

[R] gc() and gc trigger

2005-06-10 Thread Pavel Khomski
hello, the question concerning to the memory used and g.c. after having removed objects. What is wrong? bevor --- gc() used (Mb) gc trigger (Mb) max used (Mb) Ncells 313142 8.4 1801024 48.1 1835812 49.1 Vcells

Re: [R] Replacing for loop with tapply!?

2005-06-10 Thread Dimitris Rizopoulos
for the maximum you could use something like: ind[, 1] - apply(mat, 2, max) I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/16/336899 Fax:

Re: [R] discovery (was: data.frame to character)w

2005-06-10 Thread Douglas Grove
Help pages are useful, you should try them e.g. ?pi or ?LETTERS How can one discover or list all available built-in objects? On Jun 10, 2005, at 7:23 AM, Muhammad Subianto wrote: L3 - LETTERS[1:3] L10 - LETTERS[1:10] LETTERS is apparently a built-in character vector. ls() and objects

Re: [R] Abrupt shut down of R session

2005-06-10 Thread Duncan Murdoch
On 6/10/2005 12:03 PM, Park, Kyong H Mr. RDECOM wrote: Dear R users, I'm using R 2.1.0 with Windows 2000. In the middle of a R session, an error message pops up saying Rgui.exe generated errors and shuts down the R session. How can I correct this problem? Appreciate your help. Make it

Re: [R] discovery (was: data.frame to character)

2005-06-10 Thread Prof Brian Ripley
On Fri, 10 Jun 2005, Robert Citek wrote: How can one discover or list all available built-in objects? sort(unlist(lapply(search()[-1], ls))) would be a good start, but with over 2000 what are you going to do with the information? (Actually, there are more ...). LETTERS is apparently a

Re: [R] discovery (was: data.frame to character)

2005-06-10 Thread Liaw, Andy
From: Robert Citek How can one discover or list all available built-in objects? On Jun 10, 2005, at 7:23 AM, Muhammad Subianto wrote: L3 - LETTERS[1:3] L10 - LETTERS[1:10] LETTERS is apparently a built-in character vector. ls() and objects () only lists the ones I've created. Is

Re: [R] discovery (was: data.frame to character)

2005-06-10 Thread Achim Zeileis
On Fri, 10 Jun 2005 11:59:07 -0500 Robert Citek wrote: How can one discover or list all available built-in objects? Depends what you mean with `built-in'... On Jun 10, 2005, at 7:23 AM, Muhammad Subianto wrote: L3 - LETTERS[1:3] L10 - LETTERS[1:10] LETTERS is apparently a built-in

Re: [R] How to plot more than 3 sets in Venn Diagrams?

2005-06-10 Thread David Forrest
On Fri, 10 Jun 2005, Jim Lemon wrote: Tan Hui Hui Jenny wrote: I'm trying to plot Venn diagrams with more than 3 sets (5 actually) in order to describe graphically the genetic variation between populations. I tried the limma library but realised it can only plot 3 sets. Is there any

Re: [R] rgl.snapshot failed

2005-06-10 Thread Dirk Eddelbuettel
Roger, roger koenker rkoenker at uiuc.edu writes: I've installed the rgl package on a Suse x86-64 machine (further details below) and it produces nice screen images. Unfortunately, rgl.snapshot attempts to make png files produces only the response failed. For other graphics png()

Re: [R] discovery (was: data.frame to character)

2005-06-10 Thread Gabor Grothendieck
On 6/10/05, Robert Citek [EMAIL PROTECTED] wrote: How can one discover or list all available built-in objects? On Jun 10, 2005, at 7:23 AM, Muhammad Subianto wrote: L3 - LETTERS[1:3] L10 - LETTERS[1:10] LETTERS is apparently a built-in character vector. ls() and objects () only

Re: [R] discovery

2005-06-10 Thread Duncan Murdoch
On 6/10/2005 12:59 PM, Robert Citek wrote: How can one discover or list all available built-in objects? search() tells you all the attached environments. ls(n) gives you a list of what is in the n'th one. The 1300+ page R Reference Manual gives you all the help pages for base packages.

Re: [R] discovery (was: data.frame to character)

2005-06-10 Thread Gabor Grothendieck
On 6/10/05, Gabor Grothendieck [EMAIL PROTECTED] wrote: On 6/10/05, Robert Citek [EMAIL PROTECTED] wrote: How can one discover or list all available built-in objects? On Jun 10, 2005, at 7:23 AM, Muhammad Subianto wrote: L3 - LETTERS[1:3] L10 - LETTERS[1:10] LETTERS is

[R] Solved: Fortran compilation error

2005-06-10 Thread Brian Beckage
I solved this problem by installing a newer version of the fortran (g77 3.4) compiler from hpc.sourceforge.net into usr/local. Best, Brian At 12:49 PM -0400 6/10/05, Brian Beckage wrote: Hello, I'm trying to install a package that requires a Fortran compiler (Hmisc) using R CMD INSTALL. I

[R] ANOVA vs REML approach to variance component estimation

2005-06-10 Thread Adaikalavan Ramasamy
Can anyone verify my calculations below or explain why they are wrong ? I have several animals that were measured thrice. The only blocking variable is the animal itself. I am interested in calculating the between and within object variations in R. An artificial example : y - c( 2.2, -1.4,

[R] Linear regression with variation in x and y

2005-06-10 Thread Johannes Graumann
Hello, According to my archive digging, this issue comes up again and again. I was just wondering whether somebody had code lying around for a Deming Regression or a MLFR method for this, which I might peruse. I came across a post from Brian Ripley stating that something like it was close to be

Re: [R] gc() and gc trigger

2005-06-10 Thread Thomas Lumley
On Fri, 10 Jun 2005, Pavel Khomski wrote: hello, the question concerning to the memory used and g.c. after having removed objects. What is wrong? Nothing is wrong. bevor --- gc() used (Mb) gc trigger (Mb) max used (Mb) Ncells 313142

Re: [R] Plot/manage spatial boundary data

2005-06-10 Thread David Forrest
On Fri, 10 Jun 2005, Roger Bivand wrote: On Thu, 9 Jun 2005, David Forrest wrote: I have some disconnected boundary data from a finite element ocean model and I'd like to make a plot. Maptools looks promising, but since my data is not in a shapefile or a map, I'm unclear on what the

[R] Default Format for Dates?

2005-06-10 Thread khobson
Is there anyway to preset date formats? I have a date from a cover.dbf that is shown as this: cover$FINALREPOR [1] 2003-06-24 The numeric value in cover$FINALREPOR is 12227. I'd rather not create another vector to hold the properly formatted date. When I put this in a WordPerfect merge, I

[R] Estimate of baseline hazard in survival

2005-06-10 Thread Hanke, Alex
Dear All, I'm having just a little terminology problem, relating the language used in the Hosmer and Lemeshow text on Applied Survival Analysis to that of the help that comes with the survival package. I am trying to back out the values for the baseline hazard, h_o(t_i), for each event time or

Re: [R] Replacing for loop with tapply!?

2005-06-10 Thread Sander Oom
Dear all, Dimitris and Andy, thanks for your great help. I have progressed to the following code which runs very fast and effective: mat - matrix(sample(-15:50, 15 * 10, TRUE), 15, 10) mat[mat45] - NA mat-NA mat temps - c(35, 37, 39) ind - rbind( t(sapply(temps, function(temp)

Re: [R] Default Format for Dates?

2005-06-10 Thread Gabor Grothendieck
On 6/10/05, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote: Is there anyway to preset date formats? I have a date from a cover.dbf that is shown as this: cover$FINALREPOR [1] 2003-06-24 The numeric value in cover$FINALREPOR is 12227. I'd rather not create another vector to hold the

Re: [R] Default Format for Dates?

2005-06-10 Thread Roger D. Peng
Is something like this what you want? x - as.Date(2003-06-24) format(x, %B %d %Y) or perhaps as.character(x, %B %d %Y) -roger [EMAIL PROTECTED] wrote: Is there anyway to preset date formats? I have a date from a cover.dbf that is shown as this: cover$FINALREPOR [1] 2003-06-24

Re: [R] Rpy and RSPython

2005-06-10 Thread Seth Falcon
On 10 Jun 2005, [EMAIL PROTECTED] wrote: Hi, I tried to install RSPython, here is my problem: I load the env variables correctly as below: #-- # R section #-- R_HOME=/usr/lib/R

Re: [R] Estimate of baseline hazard in survival

2005-06-10 Thread Thomas Lumley
On Fri, 10 Jun 2005, Hanke, Alex wrote: Dear All, I'm having just a little terminology problem, relating the language used in the Hosmer and Lemeshow text on Applied Survival Analysis to that of the help that comes with the survival package. I am trying to back out the values for the

Re: [R] ANOVA vs REML approach to variance component estimation

2005-06-10 Thread Chuck Cleland
They look fine to me. Also, note varcomp() in the ape package and VarCorr() in the nlme package. I think in this case the ANOVA estimate of the intercept variance component is negative because the true value is close to zero. y - c( 2.2, -1.4, -0.5, # animal 1 +-0.3, -2.1,

[R] 1st derivatives using gam.

2005-06-10 Thread Bill Shipley
Hello. Using the smooth.spline function one can obtain the 1st derivatives by specifiying deriv=1 in the predict function (i.e. predict(fit,fit$x,deriv=1) where fit is the object created by smooth.spline). Can one obtain the first derivatives using gam() and predict.gam? If so, how?

Re: [R] Default Format for Dates?

2005-06-10 Thread Gabor Grothendieck
On 6/10/05, Gabor Grothendieck [EMAIL PROTECTED] wrote: On 6/10/05, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote: Is there anyway to preset date formats? I have a date from a cover.dbf that is shown as this: cover$FINALREPOR [1] 2003-06-24 The numeric value in cover$FINALREPOR

[R] mask a matrix

2005-06-10 Thread Lei Jiang
Hi, there. I have two matrix with identical dimentions. matrix A contains information of 0 and 1, and matrix B contains data. I only want the data in matrix B where matrix A has 1's. the places where matrix A has 0's, I want NA's in matrix B. How do I do that?? Thank you very much. Lei

Re: [R] mask a matrix

2005-06-10 Thread Liaw, Andy
If you just multiply the two together, you get 0 in the right positions. If you want those to be NAs, just do: result - ifelse(A, B, NA) Andy From: Lei Jiang Hi, there. I have two matrix with identical dimentions. matrix A contains information of 0 and 1, and matrix B contains data. I

Re: [R] mask a matrix

2005-06-10 Thread Mwalili, S. M.
Jiang, This example could help: A - matrix(rbinom(20,1,.8),ncol=2) B - cbind(rnorm(20,-1,0),rnorm(20,1,0)) B - ifelse(A==0,NA,B) B A S M Mwalili. Lei Jiang [EMAIL PROTECTED] wrote: Hi, there. I have two matrix with identical dimentions. matrix A contains information of 0 and 1, and matrix

[R] problem with polr ?

2005-06-10 Thread Marc Girondot
I want to fit a multinomial model with logit link. For example let this matrix to be analyzed: male female aborted factor 10 12 1 1.2 14 14 4 1.3 15 12 3 1.4 (this is an example, not the true data which are far more complex...) I

Re: [R] Default Format for Dates?

2005-06-10 Thread Gabor Grothendieck
On 6/10/05, Gabor Grothendieck [EMAIL PROTECTED] wrote: On 6/10/05, Gabor Grothendieck [EMAIL PROTECTED] wrote: On 6/10/05, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote: Is there anyway to preset date formats? I have a date from a cover.dbf that is shown as this:

[R] Performance difference between 32-bit build and 64-bit build on Solaris 8

2005-06-10 Thread Scott Gilpin
Hi everyone - I'm seeing a 32-bit build perform significantly faster (up to 3x) than a 64 bit build on Solaris 8. I'm running R version 2.1.0. Here are some of my system details, and some resulting timings: uname -a SunOS lonetree 5.8 Generic_117350-16 sun4u sparc SUNW,Sun-Fire-V440 lonetree

Re: [R] Replacing for loop with tapply!?

2005-06-10 Thread Kjetil Brinchmann Halvorsen
Sander Oom wrote: Dear all, We have a large data set with temperature data for weather stations across the globe (15000 stations). For each station, we need to calculate the number of days a certain temperature is exceeded. So far we used the following S code, where mat88 is a matrix

[R] RCMD Warnings on src directory.

2005-06-10 Thread Jagarlamudi, Choudary
Hi Group, I performed the following commands to build my package in R 2.0 under Windows XP I got all my tools from Dr. Duncan Mudroch's website. I did a RCMD build dnal and it built a tar file for me. I did a RCMD INSTALL dnal and it installed well. When i do RCMD check dnal i get the

Re: [R] problem with polr ?

2005-06-10 Thread John Fox
Dear Marc, -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Marc Girondot Sent: Friday, June 10, 2005 3:44 PM To: r-help@stat.math.ethz.ch Subject: [R] problem with polr ? I want to fit a multinomial model with logit link. For example let this

Re: [R] RCMD Warnings on src directory.

2005-06-10 Thread Thomas Lumley
On Fri, 10 Jun 2005, Jagarlamudi, Choudary wrote: Hi Group, I performed the following commands to build my package in R 2.0 under Windows XP I got all my tools from Dr. Duncan Mudroch's website. I did a RCMD build dnal and it built a tar file for me. I did a RCMD INSTALL dnal and it

[R] launching RSJava from external application

2005-06-10 Thread Sean
Hi,all, I have a problem about launching SJava in R from external application (Cytoscape). My R version is 2.1.0 and Sjava version is 0.68. R is installed by issuing 1. configure --enable-R-shlib 2. make 3. make check 4. make install SJava is installed by issuing R CMD INSTALL -c

[R] us zipcode data map

2005-06-10 Thread Mike R
i've search the email archives, searched the documention of various map packages and done an R-site search, but have been unable to find direct resources for creating maps of the US that are colored or annotated or ... by zipcode data. For example, create a map of the US and color each

Re: [R] How to read a row dataset one by one

2005-06-10 Thread Francisco J. Zagmutt
You could use by(dat,row.names(dat), print) #It will return a list with results for each row Is this what you want? Cheers Francisco From: Jan Sabee [EMAIL PROTECTED] Reply-To: Jan Sabee [EMAIL PROTECTED] To: R-help@stat.math.ethz.ch Subject: [R] How to read a row dataset one by one Date:

Re: [R] us zipcode data map

2005-06-10 Thread Francisco J. Zagmutt
library(maps) example(match.map) #for coloring If you want to annotate the map look at ?map.text I hope this helps Francisco From: Mike R [EMAIL PROTECTED] Reply-To: r-help@stat.math.ethz.ch To: r-help@stat.math.ethz.ch Subject: [R] us zipcode data map Date: Fri, 10 Jun 2005 15:18:26 -0700

Re: [R] Performance difference between 32-bit build and 64-bit bu ild on Solaris 8

2005-06-10 Thread Liaw, Andy
I'm not familiar with Solaris, so take this with appropriate dose of NaCl... For the 64-bit build, why not have the -O2 for gcc, since you have it for g77 and g++? If you just run vanilla configure for the 32-bit build, I believe it uses -O2 for all three compilers by default. If that's the

Re: [R] us zipcode data map

2005-06-10 Thread Mike R
On 6/10/05, Francisco J. Zagmutt [EMAIL PROTECTED] wrote: library(maps) example(match.map) #for coloring If you want to annotate the map look at ?map.text thanks Francisco, correct me if i am wrong, but maps_2.0-27.tar.gz does many many maps, but not any zipcode maps ?

Re: [R] Replacing for loop with tapply!?

2005-06-10 Thread Adaikalavan Ramasamy
OK, so you want to find some summary statistics for each column, where some columns could be completely missing. Writing a small wrapper should help. When you use apply(), you are actually applying a function to every column (or row). First, let us simulate a dataset with 15 days/rows and 10

Re: [R] us zipcode data map

2005-06-10 Thread Mike R
i've found a 3-column text file of zip, long, lat from the 2000 US census, so i am part way there. http://www.cryptnet.net/fsp/zipdy/ __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] us zipcode data map

2005-06-10 Thread Mike R
i'm likely looking for something that does not exist: http://www.census.gov/cgi-bin/geo/tigerfaq?Q16 ==begin quote Note that the task of creating a lat/long or polygon file of all ZIP Codes is not as easy as it seems since ZIP Codes are not designed to be polygons and can't easily be forced