[R] Rmetrics - R-Forge - Workshop

2008-03-18 Thread chalabi
Dear Members of the R-Core Team, Rmetrics Developers, and Rmetrics Users ... The repository of the development version of the Rmetrics software environment has been moved to R-forge. The new R-Forge framework for R-project developers based on GForge

Re: [R] problems with rgl in Ubuntu 'gutsy'

2008-03-18 Thread Hans Ekbrand
On Mon, Mar 17, 2008 at 09:11:05PM -, Foadi, J (James) wrote: On Mon, Mar 17, 2008 at 05:43:58PM -, Foadi, J (James) wrote: [...] | all rgl windows device do not show the outer frame, so I cannot move or resize the window. | In addition to that, the image is like frozen, it can't

[R] Partition data into clusters

2008-03-18 Thread Suhaila Zainudin
Greetings R-users, I have been using the fpc package in R to cluster my data. Speficically I am using kmeansruns clustering. I would like to know how I use R to partition data into clusters. What I am doing is as follows. # Use csv file as input # wholeset =

Re: [R] Saving Chinese text in Quartz

2008-03-18 Thread Prof Brian Ripley
Which version of R is this? -- there are two completely different quartz() devices. (Please remember to supply the 'at a minimum' information the posting guide asked for.) You haven't even told us your OS, but we can deduce that this is some version of Mac OS X. You would do better to ask

[R] re ad .csv separated with coma

2008-03-18 Thread mysimbaa
Dear R users, I try to read .csv composed of two rows which are separated with coma in the followin manner: 1.ignore the three first lines 2. read datas separated with coma. I tried the following : pathname - file.choose() z -read.csv(pathname,

Re: [R] Table of basic descriptive statistics like SPSS

2008-03-18 Thread Philippe Grosjean
Jim Lemon wrote: [EMAIL PROTECTED] wrote: Dear list readers, I want to: 1. Get a table of basic descriptive statistics for my variables with the variable names one below the other like SPSS descriptive statistics: Varname N Min Max Mean SD x x xx x xxx x x xx

[R] two cols in a data frame are the same factor

2008-03-18 Thread Andres Legarra
Dear all, I have a data set (QTL detection) where I have two cols of factors in the data frame that correspond logically (in my model) to the same factor. In fact these are haplotype classes. Another real-life example would be family gas consumption as a function of car company (e.g. Ford, GM, and

[R] R plot as image data and another question

2008-03-18 Thread Bio7
Dear R developers, i would like to get the plots of R in the form of matrix data which i want to transfer to an external image software. Is it generally possible to get the plot information in form of data values? Another question concerns about the command line execution of scripts. From an

Re: [R] re ad .csv separated with coma

2008-03-18 Thread Eik Vettorazzi
You should skip the first 3 rows by using option skip=3 and set header=FALSE, which defaults to TRUE. But you could also use the 3rd row as header, doing so by setting skip=2. z -read.csv(pathname, skip=3, header=F) or z -read.csv(pathname, skip=2) (By the way, it is obvious not necessary to

Re: [R] R plot as image data and another question

2008-03-18 Thread Prof Brian Ripley
On Tue, 18 Mar 2008, Bio7 wrote: Dear R developers, i would like to get the plots of R in the form of matrix data which i want to transfer to an external image software. Is it generally possible to get the plot information in form of data values? Another question concerns about the

Re: [R] odbcQuery , memory.size

2008-03-18 Thread Prof Brian Ripley
On Tue, 18 Mar 2008, Lukas Rohrer wrote: Dear R cracks I am trying to fetch 38 Tables from 38 ESRI Geodatabases through an ODBC connection. On Windows, using RODBC? I stored the 38 channels in a list and the 38 tablenames in an other list. Why? It looks like you only need the channel

Re: [R] empty array

2008-03-18 Thread Robin Hankin
Hello everyone I know other, more knowledgeable, people have replied to Christophe's question, but perhaps the List would be interested to know that zero-extent arrays are useful (to me at least) because although such an array has no content, the dimname are nevertheless retained: a -

Re: [R] R plot as image data and another question

2008-03-18 Thread Bio7
In my case i use R by means of Rserve from my Java application. Since i`ve integrated the well known image analysis tool ImageJ i can send images to R and vice versa. This works well and is quite fast. My idea now was to get the plot information from R in the form of matrix data which i can

[R] Gaining access to complex LAPACK functions

2008-03-18 Thread Rishi
Hello, I am using R version 2.6.2 on Windows XP. I installed R using the installation file provided through one of the CRAN mirrors; I did not build it myself. as a result, the LAPACK functionality is limited. I'm trying to install a package that requires access to complex LAPACK functions

Re: [R] grouped colSums without for loops?

2008-03-18 Thread Albert Greinoecker
try: aggregate(d[,2:3], by=list(d$foo), FUN=sum) cheers, Albert Am Dienstag, den 18.03.2008, 11:41 +0100 schrieb Hans Ekbrand: Consider the following data.frame, d. d - data.frame(c(a,b,c,d,b,a), c(1,4,85,3,4,1), c(7,6,2,4,15,33)) names(d) - c(foo, bar, baz) To get the colsum of d[[2]]

Re: [R] grouped colSums without for loops?

2008-03-18 Thread Hans Ekbrand
On Tue, Mar 18, 2008 at 05:57:59AM -0500, jim holtman wrote: Is this what you want? lapply(split(d, d$foo), function(x) colSums(x[,-1])) Yes! Thank you! the *apply functions seem very powerful, thanks again for giving me a hint on how to use them. -- Hans Ekbrand (http://sociologi.cjb.net)

Re: [R] rgl build warnings and loading error on Linux

2008-03-18 Thread Liviu Andronic
On Mon, Mar 17, 2008 at 11:07 PM, Liviu Andronic [EMAIL PROTECTED] wrote: This is what I have on my system: localhost liviu # locate libGL.so /usr/lib/opengl/xorg-x11/lib/libGL.so.1 /usr/lib/opengl/xorg-x11/lib/libGL.so.1.2 /usr/lib/opengl/xorg-x11/lib/libGL.so /usr/lib/libGL.so

[R] How to change symbols in a trellis graph (UNCLASSIFIED)

2008-03-18 Thread Park, Kyong H Mr ECBC
Classification: UNCLASSIFIED Caveats: NONE Hello, R Users, I created a trellis graph following methods suggest in Mixed-Effects Models in S and S-Plus by Bates and Pinheiro using 'groupedData' and 'outer' commands as follows: 1. data2-groupedData(LCU~ST1|Blcok, data=data1) 2. plot(data2,

[R] How to get week of the year, scale 1-52?

2008-03-18 Thread Lauri Nikkinen
R users, I have a vector of dates days - seq(as.Date(2007/1/1), as.Date(2008/1/31), days) and I would like to have week numbers from 1 to 52 for each year. How do I do that? Now I get 00-53 using format(days, %W) range(format(days, %W)) [1] 00 53 I have read Date and Time Classes in R (R

Re: [R] Table of basic descriptive statistics like SPSS

2008-03-18 Thread Frank E Harrell Jr
William Revelle wrote: At 7:05 PM +0100 3/17/08, [EMAIL PROTECTED] wrote: Dear list readers, I want to: 1. Get a table of basic descriptive statistics for my variables with the variable names one below the other like SPSS descriptive statistics: Varname N Min Max Mean SD x x x

[R] two cols in a data frame are the same factor

2008-03-18 Thread Andres Legarra
Dear all, I have a data set (QTL detection) where I have two cols of factors in the data frame that correspond logically (in my model) to the same factor. In fact these are haplotype classes. Another real-life example would be family gas consumption as a function of car company (e.g. Ford,

[R] UNSOLITED E_MAILS: Integrate R data-analysis projects with Microsoft Office for free

2008-03-18 Thread Gorden T Jemwa
Dear R Admins, I received an unsolicited e-mail from BlueInference as an R user. Does it mean that R that our e-mails (and names) is sharing it's user database with third parties without our consent? Or perhaps the BlueInference guys are using an e-mail address miner to get our contact

Re: [R] UNSOLITED E_MAILS: Integrate R data-analysis projects wi

2008-03-18 Thread Ted Harding
On 18-Mar-08 12:08:44, Gorden T Jemwa wrote: Dear R Admins, I received an unsolicited e-mail from BlueInference as an R user. Does it mean that R that our e-mails (and names) is sharing it's user database with third parties without our consent? Or perhaps the BlueInference guys are

Re: [R] UNSOLITED E_MAILS: Integrate R data-analysis projects wi

2008-03-18 Thread Doran, Harold
Can a CAPTCHA be implemented as a prevenative measure -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of [EMAIL PROTECTED] Sent: Tuesday, March 18, 2008 7:33 AM To: Gorden T Jemwa Cc: r-help@r-project.org Subject: Re: [R] UNSOLITED E_MAILS:

Re: [R] How to get week of the year, scale 1-52?

2008-03-18 Thread Gabor Grothendieck
On Tue, Mar 18, 2008 at 7:55 AM, Lauri Nikkinen [EMAIL PROTECTED] wrote: R users, I have a vector of dates days - seq(as.Date(2007/1/1), as.Date(2008/1/31), days) and I would like to have week numbers from 1 to 52 for each year. How do I do that? Now I get 00-53 using format(days, %W)

[R] Re move a number from a vector

2008-03-18 Thread Josh Roofchop
This should be an easy one, but I haven't been able to figure it out, x-c(2, 4, 6, 9, 10) #how do I remove a single value x-x[-2] # will remove the second element of the vector not the number 2 This is part of a loop that looks a little like this... x-c(2, 4, 6, 9, 10) remove.value-4 #then I

[R] Time scale x-axis, Reformat data.table

2008-03-18 Thread Stephan K.
Hello, i'm pretty new on R and i'm sure the problem i have is easy to fix. my data table looks like 223344 556677.3 123455 665545.3 . . and so one i wanna plot that table but the firts column is a time with hhmmss. So i think i have to reformat the column (and maybe give

Re: [R] Re move a number from a vector

2008-03-18 Thread John Kane
?which x - x[which(x==remove.value)] --- Josh Roofchop [EMAIL PROTECTED] wrote: This should be an easy one, but I haven't been able to figure it out, x-c(2, 4, 6, 9, 10) #how do I remove a single value x-x[-2] # will remove the second element of the vector not the number 2 This

Re: [R] How to get week of the year, scale 1-52?

2008-03-18 Thread Liaw, Andy
Looks like you can convert to Julian date with the right origin and then divide by 7, and take modulo 52. Something like: julian(days, origin=days[1]-1) %/% 7 %% 52 That gets you 0-51. You just need to add one to get what you want. (I subtract one from days[1] because it's a Monday. I'm

Re: [R] Re move a number from a vector

2008-03-18 Thread Albert Greinoecker
does x[!x == remove.value] do what you need? [1] 2 6 9 10 regards, Albert Am Dienstag, den 18.03.2008, 05:21 -0700 schrieb Josh Roofchop: This should be an easy one, but I haven't been able to figure it out, x-c(2, 4, 6, 9, 10) #how do I remove a single value x-x[-2] # will remove

Re: [R] grouped colSums without for loops?

2008-03-18 Thread Liaw, Andy
See ?rowsum. The other two suggestions are both shown in the Example section as well (with comment). Andy From: Hans Ekbrand * PGP Signed by an unknown key: 03/18/2008 at 06:22:03 AM On Tue, Mar 18, 2008 at 05:57:59AM -0500, jim holtman wrote: Is this what you want? lapply(split(d,

Re: [R] How to get week of the year, scale 1-52?

2008-03-18 Thread John Kane
365/7 != 52. You have more than 52 weeks in each year. --- Lauri Nikkinen [EMAIL PROTECTED] wrote: R users, I have a vector of dates days - seq(as.Date(2007/1/1), as.Date(2008/1/31), days) and I would like to have week numbers from 1 to 52 for each year. How do I do that? Now I

Re: [R] Windows Vista 64-bit support

2008-03-18 Thread Daniel Gatti
OK, so there is a 32-bit version of R that runs on Vista, but no 64-bit version, due to the lack of 64-bit compilers. Is there a forecast for when a 64-bit version of R might compile in Vista? Dan From: chaogai [EMAIL PROTECTED] The standard windows build works in Vista 64. A 64 bits build

Re: [R] How to get week of the year, scale 1-52?

2008-03-18 Thread Lauri Nikkinen
Yeah, that's correct. I can accept that there is more than 52 weeks in a year (:-)) but it's hard to accept that the first week e.g in this example is 00 (year 2008). Should the first week be the first one, 1? -Lauri 2008/3/18, John Kane [EMAIL PROTECTED]: 365/7 != 52. You have more than 52

[R] Walds Z in a mixed model

2008-03-18 Thread Daniel Bolnick
I recently reviewed a paper that used Wald's Z in a mixed model analysis. The dataset (see end of message for complete dataset) has individual animals (SAMPLE) from within social groups, from within geographic regions. Y is the response variable. The paper ran a mixed model with year and

Re: [R] How to get week of the year, scale 1-52?

2008-03-18 Thread John Kane
--- Lauri Nikkinen [EMAIL PROTECTED] wrote: Yeah, that's correct. I can accept that there is more than 52 weeks in a year (:-)) but it's hard to accept that the first week e.g in this example is 00 (year 2008). Should the first week be the first one, 1? I'd think so but I don't know much

Re: [R] how to plot a map on a non-rectilinear grid

2008-03-18 Thread jgarcia
There is another solution, which depending on the resolution and amount of your data may be useful. You can find through interpolation the corners of four coordinates around any data point that define a four sides polygon. This grid so created does not need to be regular. Then, you can simple call

Re: [R] problems with rgl in Ubuntu 'gutsy'

2008-03-18 Thread Ben Bolker
Hans Ekbrand hans at sociologi.cjb.net writes: [snip snip snip] On my system that window can be moved and even resized without problems (keeping a high rendering rate). Sorry, but I my test didn't give any conclusive new facts about the cause of the problem. If you have fancy 3D

Re: [R] UNSOLITED E_MAILS: Integrate R data-analysis projects wi

2008-03-18 Thread Douglas Bates
Usually a captcha is used to prevent creation of email accounts for use by spammers. (There was an interesting article recently on whether the Gmail captcha scheme had been broken so that spammers could create masses of gmail accounts. The general conclusion is that the capcha scheme is intact

[R] ask for help how to get accumulate contingency table?

2008-03-18 Thread alex chin
Hi, I post the question again. x-rep(c(2,3,4),20) y-sample(rep(c(0,1),30)) tr-rep(c(1:5),6) data1-data.frame(x,y,tr) data1-data1[order(data1$tr),] data1 x y tr 2 1 1 4 0 1 3 1 1 2 1 1 4 1 1 3 1 1 2 0 1 4 0 1 3 0 1 2 0 1 4 0 1 3 1 1 3 0 2 2 0 2 table(data1) tr = 1 y x 0 1

[R] Splitting data

2008-03-18 Thread Kathy Maher
Does anyone know a loop that will allow me to split rows of data to get sum of squares; and then to move rows around to minimize the sum of squares? I am stuck on how to split and move rows around. [[alternative HTML version deleted]] __

Re: [R] How to get week of the year, scale 1-52?

2008-03-18 Thread Eik Vettorazzi
But usually weeks start with a Monday. So what should happen when new year isn't on a Monday? format(days, %W) puts the first Monday of the year in the first week, which makes sense. hth. Lauri Nikkinen schrieb: Yeah, that's correct. I can accept that there is more than 52 weeks in a year

Re: [R] Re move a number from a vector

2008-03-18 Thread Josh Roofchop
Thanks to both John and Arthur for their replies, they both worked great. Josh Josh Roofchop wrote: This should be an easy one, but I haven't been able to figure it out, x-c(2, 4, 6, 9, 10) #how do I remove a single value x-x[-2] # will remove the second element of the vector not the

[R] read.dta for files from stata 9.0

2008-03-18 Thread Albrecht Kauffmann
Dear R-helpers, if I want to read a .dta-file generated by stata 9.0 with read.dta (foreign), I get the message not a stata version 5-8 .dta-file. I'm using R-2.6.2 and the latest version of the foreign package. Has someone any hint? With many thanks, Albrecht

[R] how to reset slogic.f file

2008-03-18 Thread coldeyes Liao
Hi there: recently i try to use LogicReg package for a tree model(logistics fit ) . i list my code and error below: dim(model.dat) [1] 48000 745 fit1 - logreg(resp = model.dat[,745], bin=model.dat[, 9:700], type = 3, select = 3, ntrees = c(1,2), nleaves=c(1,7), ) Insufficient declaration

Re: [R] read.dta for files from stata 9.0

2008-03-18 Thread K. Elo
Hi, Albrecht Kauffmann wrote (18.03.2008): Dear R-helpers, if I want to read a .dta-file generated by stata 9.0 with read.dta (foreign), I get the message not a stata version 5-8 .dta-file. I'm using R-2.6.2 and the latest version of the foreign package. Has someone any hint? Seems quite

Re: [R] problems with rgl in Ubuntu 'gutsy'

2008-03-18 Thread Foadi, J (James)
If you have fancy 3D graphics effects turned on, try turning them off. (System - Preferences - Appearance - Visual Effects - None (Provides a simple desktop appearance ...) unfortunately rgl and compiz/etc. both try to use the same OpenGL interface, so you can't use both at the same

Re: [R] UNSOLITED E_MAILS: Integrate R data-analysis projects wi

2008-03-18 Thread Doran, Harold
I think on the python list, when you review the archives, the poster address is viewed like a CAPTCHA. So, it makes it slightly more difficult (though not impossible) to pullout poster emails addresses and replace john.doe at domainname.com -Original Message- From: [EMAIL PROTECTED]

[R] Puzzled at generating combinations

2008-03-18 Thread Ajay Shah
I have two data frames. Suppose the first has rows r1 r2 r3 and the second has rows R1 R2 R3 I'd like to generate the data frame: r1 R1 r1 R2 r1 R3 r2 R1 r2 R2 r2 R3 r3 R1 r3 R2 r3 R3 How would I go about doing this? I'm sure there's a

Re: [R] problem with merge

2008-03-18 Thread Don MacQueen
Do any of your Entrez.Gene values appear in more than one row in either of the dataframes? For example: dfa - data.frame( a = c('A','A','B','C') , x= 1:4) dfb - data.frame( a = c('A','B','B','D'), y=1:4) merge(dfa,dfb) a x y 1 A 1 1 2 A 2 1 3 B 3 2 4 B 3 3 length(intersect(dfa$a,

[R] Error when producing multiple graphes

2008-03-18 Thread Uli Kleinwechter
Dear all, Within a larger script I use the following code to produce barplots for all variables contained in a dataframe type3m: for(xn in names(type3m)) barplot(table(type3m[[xn]]), main = paste (Barplot of, xn), width=0.75,legend.text=barnamesm, xlim=c(0,4),ylim=c(0,400)) The

Re: [R] read.dta for files from stata 9.0

2008-03-18 Thread Stefan Grosse
On Tuesday 18 March 2008 03:33:02 pm K. Elo wrote: KE Seems quite obvious to me: not a stata version 5-8 means, that the 9.0 KE version is not supported. So You should save the file in stata ver. 8 KE format and then give a new try. KE KE Hope this helps, No. I have stata 9.2 and the data import

[R] download.file()

2008-03-18 Thread Paul Evans
Hi, I wanted to download a file and did the following: - fileLink - 'ftp://ftp.ncbi.nih.gov/pub/geo/DATA/supplementary/series/GSE1000/GSE1000_RAW.tar' download.file(fileLink,'/geoDat') trying URL

Re: [R] download.file()

2008-03-18 Thread Henrik Bengtsson
On Tue, Mar 18, 2008 at 8:46 AM, Paul Evans [EMAIL PROTECTED] wrote: Hi, I wanted to download a file and did the following: - fileLink - 'ftp://ftp.ncbi.nih.gov/pub/geo/DATA/supplementary/series/GSE1000/GSE1000_RAW.tar'

Re: [R] download.file()

2008-03-18 Thread Stefan Grosse
On Tuesday 18 March 2008 04:46:12 pm Paul Evans wrote: PE download.file(fileLink,'/geoDat') I would have expected download.file(fileLink,/geoDat/yourfilename) note that you need your complete path (or use setwd) -- Microeconomics University of Erfurt signature.asc Description: This is a

Re: [R] download.file()

2008-03-18 Thread Henrik Bengtsson
On Tue, Mar 18, 2008 at 8:57 AM, Henrik Bengtsson [EMAIL PROTECTED] wrote: On Tue, Mar 18, 2008 at 8:46 AM, Paul Evans [EMAIL PROTECTED] wrote: Hi, I wanted to download a file and did the following: - fileLink -

Re: [R] download.file()

2008-03-18 Thread Paul Evans
Aaah...I see now. Thanks !! - Original Message From: Henrik Bengtsson [EMAIL PROTECTED] To: Paul Evans [EMAIL PROTECTED] Cc: [EMAIL PROTECTED] Sent: Tuesday, March 18, 2008 11:59:45 AM Subject: Re: [R] download.file() On Tue, Mar 18, 2008 at 8:57 AM, Henrik Bengtsson [EMAIL PROTECTED]

Re: [R] download.file()

2008-03-18 Thread Prof Brian Ripley
/geodat *is* an absoute path! The second argument of download.file is called 'destfile', so what makes you think it is a desination *directory*? Your command works for me, using an account that has write permission in / , and not otherwise. (It seems very unsafe that you would be using such

Re: [R] download.file()

2008-03-18 Thread Stefan Grosse
On Tuesday 18 March 2008 04:59:45 pm Henrik Bengtsson wrote: HB Also, you want to download the file in a binary fashion, i.e. use HB argument mode=wb, otherwise your binary tar file will be corrupt. HB HB R-core: I'd suggest to replace the default to mode=wb for file HB transfers. I had no

Re: [R] Puzzled at generating combinations

2008-03-18 Thread Richard . Cotton
I have two data frames. Suppose the first has rows r1 r2 r3 and the second has rows R1 R2 R3 I'd like to generate the data frame: r1 R1 r1 R2 r1 R3 r2 R1 r2 R2 r2 R3 r3 R1 r3 R2 r3 R3 Try: col1 - paste(r, 1:3, sep=)

Re: [R] download.file()

2008-03-18 Thread Henrik Bengtsson
On Tue, Mar 18, 2008 at 9:10 AM, Prof Brian Ripley [EMAIL PROTECTED] wrote: /geodat *is* an absoute path! The second argument of download.file is called 'destfile', so what makes you think it is a desination *directory*? The help file might be the source of confusion: destfile: A character

[R] lme library

2008-03-18 Thread roberta_varriale
Dear authors, I’m an Italian PhD student and I’m dealing with linear mixed models. I’d like to use your lme library, but I have a problem. I’m able to estimate a null model. For example, using SAS data, I can estimate the model: lme(y ~ 1, data = Mississippi, random = ~ 1|influent, method=ML)

Re: [R] read.dta for files from stata 9.0

2008-03-18 Thread Prof Brian Ripley
On Tue, 18 Mar 2008, Stefan Grosse wrote: On Tuesday 18 March 2008 03:33:02 pm K. Elo wrote: KE Seems quite obvious to me: not a stata version 5-8 means, that the 9.0 KE version is not supported. So You should save the file in stata ver. 8 KE format and then give a new try. KE KE Hope this

Re: [R] Puzzled at generating combinations

2008-03-18 Thread Henrique Dallazuanna
Try: expand.grid(col2, col1) On 18/03/2008, Ajay Shah [EMAIL PROTECTED] wrote: I have two data frames. Suppose the first has rows r1 r2 r3 and the second has rows R1 R2 R3 I'd like to generate the data frame: r1 R1 r1 R2 r1 R3 r2 R1 r2 R2 r2 R3

Re: [R] Puzzled at generating combinations

2008-03-18 Thread John Kane
Will something like this work? x - data.frame(xx=c(r1, r2, r3), yy=c(1,2,3)) y - data.frame(aa=c(R1, R2, R3), bb=c(10,12, 14)) merge(x,y) --- Ajay Shah [EMAIL PROTECTED] wrote: I have two data frames. Suppose the first has rows r1 r2 r3 and the second has rows R1 R2 R3 I'd

[R] Automatic command line options

2008-03-18 Thread Aaron MacNeil
Hello users, Does anyone know how to turn off the automatic double quoting and bracketing on the command line that appeared in R 2.6.x (OS X). It's driving me nuts! Many thanks, Aaron. M. Aaron MacNeil __

Re: [R] Table of basic descriptive statistics like SPSS

2008-03-18 Thread ukoenig
Dear all, thank you very much for the provided support and all the hints, which were very helpful! Now I know how to continue. See ?pdf Please, note first search commands: apropos(pdf) RSiteSearch(latex) Using these, you would have found by yourself the various corresponding functions. Philippe

Re: [R] read.dta for files from stata 9.0

2008-03-18 Thread Stefan Grosse
On Tuesday 18 March 2008 05:39:10 pm Prof Brian Ripley wrote: PB Did you check the sources? The source for read.dta explicitly checks for PB VERSION_5 ... VERSION_8 and otherwise says PB PB not a Stata version 5-8 .dta file No, I did not check the sources, but I checked the help file and there

[R] Trouble fitting models with nonNULL datadist in Design Library

2008-03-18 Thread Holder, Daniel J.
Having trouble with datadist in the Design library. When I set the datadist option, all of the fitting functions give error messages similar to Error in exists(x, envir, mode, inherits) : invalid first argument. Code below was run in R version 2.5.0 under windows, but get similar results with R

Re: [R] read.dta for files from stata 9.0

2008-03-18 Thread Stefan Grosse
On Tuesday 18 March 2008 05:39:10 pm Prof Brian Ripley wrote: PB According to http://www.stata.com/help.cgi?dta, Stata now writes in PB 'format-114'. read.dta expects only versions 113, 111 and earlier. PB There is very little difference between formats 113 and 114, so it should PB be easy to

Re: [R] rgl build warnings and loading error on Linux

2008-03-18 Thread Charles C. Berry
On Tue, 18 Mar 2008, Liviu Andronic wrote: On Mon, Mar 17, 2008 at 11:07 PM, Liviu Andronic [EMAIL PROTECTED] wrote: This is what I have on my system: localhost liviu # locate libGL.so /usr/lib/opengl/xorg-x11/lib/libGL.so.1 /usr/lib/opengl/xorg-x11/lib/libGL.so.1.2

[R] Course: R/Splus Fundamentals and Programming Techniques.***April 2008*** in San Francisco and New York City

2008-03-18 Thread [EMAIL PROTECTED]
We've added 2 additional courses to our April Schedule: (++) R/S Fundamentals and Programming Techniques *** San Francisco / April 24-25, 2008 *** *** New York City / April 28-29, 2008 *** *** Seattle / April 21-22 Regards - Sue Original Message Subject: [BioC]

Re: [R] download.file()

2008-03-18 Thread Paul Evans
Thanks Henrik. Also, is there a method/package (in R) with which I can unpack a tar file programatically? thanks again! - Original Message From: Henrik Bengtsson [EMAIL PROTECTED] To: Prof Brian Ripley [EMAIL PROTECTED] Cc: Paul Evans [EMAIL PROTECTED]; [EMAIL PROTECTED] Sent:

[R] Symetric matrix as random effect in lmer

2008-03-18 Thread Marc
UseRs, I am trying to fit the following model aa.lmer - lmer(y ~ x + (1 | s), na.action=na.omit, family=binomial) where s is a length(x) by length(x) matrix and I get this error: Error in eval(expr, envir, enclos) : Ztl[[1]] must have 28392 columns I am wondering if it is possible to fit

Re: [R] download.file()

2008-03-18 Thread Stefan Grosse
On Tuesday 18 March 2008 06:32:18 pm Paul Evans wrote: PE Also, is there a method/package (in R) with which I can unpack a tar file PE programatically? I guesses that question when I saw your tar file. Not that you not only need to un-tar but also to un-zip then several other files in that

Re: [R] download.file()

2008-03-18 Thread Paul Evans
Stefan, Correct. However I have a large number of tar files to download (with several zip files in each) - so doing it manually is going to be somewhat difficult. I may have to package this code for other users and if possible, would only like to use R for these operations. thanks for the

Re: [R] [PS] Re: Table of basic descriptive statistics like SPSS

2008-03-18 Thread Frank E Harrell Jr
Ben Fairbank wrote: tiny high-resolution histograms Could you say a bit more on that topic? How achieved? Where latex(describe()) uses LaTeX's picture environment after binning into 101 bins. documented? Similar to Tufte's sparklines? yes, a bit See the link from

Re: [R] lme library

2008-03-18 Thread Ben Bolker
roberta_varriale at ds.unifi.it writes: I’m able to estimate a null model. For example, using SAS data, I can estimate the model: lme(y ~ 1, data = Mississippi, random = ~ 1|influent, method=ML) As suggested in the literature I want to “test” the significance of the second level variance.

[R] MENU TCL-TK

2008-03-18 Thread ermimi
Hello! I need help in a menu in TCL-TK. I have 3 functions, one function (function1) that sum 2 numbers and return the result. Other function(function2) that print a number. And the other (function3), is a menu that has several options. I need use the result that returns the option that calls to

Re: [R] MENU TCL-TK

2008-03-18 Thread Gabor Grothendieck
Make sure d has been created (see first line marked ##) and then use d - function(a, b) as shown in second line marked ##. See ?- Another solution is to use the proto package. The home page http://rproto.googlecode.com has pointers to examples of using gWidgets with proto and the same idea

[R] S4 : dumpMethod

2008-03-18 Thread Christophe Genolini
Hi the list, I almost finished the green book (chapter 7 and 8 on S4), this is probably one of my last question :-) I hope you will not miss them to much ! Here is my today nighmare: I do not manage to use dumpMethod : setClass(A,representation(a=numeric))

[R] useR! 2008 submission deadline

2008-03-18 Thread Achim Zeileis
Dear useRs, this is just a quick reminder that the submission deadline for useR! 2008 at TU Dortmund is approaching: Submission deadline: 2008-03-31 (and along with that the early registration deadline). Both, online submission and registration is available on the conference Web page at

[R] invalid \uxxxx sequence while trying to create a package

2008-03-18 Thread Charles Annis, P.E.
Greetings, R-Helpers: I am trying to build a package on WinXP, something that I have done before, but infrequently, and have recently updated both the Rtools and the version of R being used. My test case fails after successfully creating the compiled html files with this error message: .

Re: [R] MENU TCL-TK

2008-03-18 Thread ermimi
My doubts were resolved, thank you Gabor ermimi wrote: Hello! I need help in a menu in TCL-TK. I have 3 functions, one function (function1) that sum 2 numbers and return the result. Other function(function2) that print a number. And the other (function3), is a menu that has

[R] How to reverse colors in filled.contour?

2008-03-18 Thread Ravi Varadhan
Hi, I am plotting values of log(hazards ratio) as a function of two predictors, using the plotting function filled.contour(). Here is a simple simulated example of this: x - seq(0,1, length=20) y - seq(0,1, length=20) z - outer(x,y, function(x,y) x^2 + y^2 ) zmat - matrix(

Re: [R] invalid \uxxxx sequence while trying to create a package

2008-03-18 Thread Duncan Murdoch
On 18/03/2008 5:17 PM, Charles Annis, P.E. wrote: Greetings, R-Helpers: I am trying to build a package on WinXP, something that I have done before, but infrequently, and have recently updated both the Rtools and the version of R being used. My test case fails after successfully

Re: [R] How to reverse colors in filled.contour?

2008-03-18 Thread Bert Gunter
from the Help file: color.palette: a color palette function to be used to assign colors in the plot. so color = function(x)rev(heat.colors(x)) should do it. -- Bert Gunter Genentech -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Ravi Varadhan Sent:

Re: [R] grouped colSums without for loops?

2008-03-18 Thread Hans Ekbrand
On Tue, Mar 18, 2008 at 12:05:23PM +0100, Albert Greinoecker wrote: try: aggregate(d[,2:3], by=list(d$foo), FUN=sum) Great! Now I can get a data.frame as well as a list, thanks! -- Hans Ekbrand (http://sociologi.cjb.net) [EMAIL PROTECTED] Signature generated by Signify v1.14. For

[R] rmultinomial() function

2008-03-18 Thread Mary Black
After scouring the online R resources and help pages, I still need clarification on the function rmultinomial(). I would like to create a vector, say of 100 elements, where every element in the vector can take on the value of 0, 1 or 2, and where each of those values have a specific

[R] ROracle installation problem

2008-03-18 Thread rmail ye
Hi there, I tried to install ROracle package (Windows XP pro), got the following errors; *C:\Program Files\R\R-2.6.2\binrcmd install c:\ROracle_0.5-9.tar.gz* * installing to 'c:/PROGRA~1/R/R-26~1.2/library'* *-- Making package roracle * *

Re: [R] rmultinomial() function

2008-03-18 Thread Erik Iverson
Mary Black wrote: After scouring the online R resources and help pages, I still need clarification on the function rmultinomial(). I would like to create a vector, say of 100 elements, where every element in the vector can take on the value of 0, 1 or 2, and where each of those values have a

[R] Tcl/tk question solved

2008-03-18 Thread Erin Hodgess
Whoops! I found it...James Wettenhall. sorry for the trouble. erin -- Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown mailto: [EMAIL PROTECTED] __ R-help@r-project.org mailing list

Re: [R] Tcl/tk question

2008-03-18 Thread Ben Tupper
Hi, I'm just a newbie, but I don't mind stealing the gurus' thunder. Try RSiteSearch(Tcl/tk James) - it brings up a number of posts authored by a James W. Cheers, Ben On Mar 18, 2008, at 7:27 PM, Erin Hodgess wrote: Dear R Gurus: What is the name of the person who has all of the Tcl/tk

Re: [R] Tcl/tk question

2008-03-18 Thread John Fox
Dear Erin, The name you're looking for is James Wettenhall, and the examples are at http://bioinf.wehi.edu.au/~wettenhall/RTclTkExamples/. Regards, John John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario, Canada L8S 4M4

[R] comparing length-weight regressions

2008-03-18 Thread saduria
I'd like to compare length-weight regressions among years. Any information would be appreciated. a. gray fisheries consultant _ 08 [[alternative HTML version deleted]] __

Re: [R] Tcl/tk question

2008-03-18 Thread Gabor Grothendieck
Philip Grosjean maintains the examples. You can find them on his site now: http://www.sciviews.org/_rgui/tcltk On Tue, Mar 18, 2008 at 7:27 PM, Erin Hodgess [EMAIL PROTECTED] wrote: Dear R Gurus: What is the name of the person who has all of the Tcl/tk stuff, please? I know it's James W,

[R] Tcl/tk question

2008-03-18 Thread Erin Hodgess
Dear R Gurus: What is the name of the person who has all of the Tcl/tk stuff, please? I know it's James W, but can't remember his last name. I wanted to look at some of his examples. Thanks in advance! Sincerely, Erin -- Erin Hodgess Associate Professor Department of Computer and

Re: [R] invalid \uxxxx sequence while trying to create a package

2008-03-18 Thread Charles Annis, P.E.
Thank you Duncan: I found the culprit - my menu contains lines like this one: winMenuAdd(menuname = mh1823/\u00e2 vs a/2. Build \u00e2 vs a Linear Model) This does produce a menu item that looks like â vs a and since R doesn't seem to like the a-hat character I tried to get around the

[R] ls() and classes

2008-03-18 Thread Erin Hodgess
Dear R People: I want to get the class of all of the objects in my directory. I was trying: do.call(class,list=ls()) but got an unused argument error. I'm sure it's simple, but I'm just not seeing it. Any help would be much appreciated. Sincerely, Erin -- Erin Hodgess Associate Professor

Re: [R] ls() and classes

2008-03-18 Thread Rolf Turner
On 19/03/2008, at 4:39 PM, Erin Hodgess wrote: Dear R People: I want to get the class of all of the objects in my directory. I was trying: do.call(class,list=ls()) but got an unused argument error. I'm sure it's simple, but I'm just not seeing it. Any help would be much appreciated.

  1   2   >