Re: [R] data.frame size limit

2004-06-09 Thread Prof Brian Ripley
On Tue, 8 Jun 2004, Philip Sobolik wrote: Is there a limit to the number of columns that a data.frame can have? Yes. A data frame is a list, and a list is limited to 2^31-1 items. For example, can I read.csv() a file that has 1000 columns and 10,000 rows, will it break or is it limited

[R] Help with a Lattice plot that fails with an empty unique combination

2004-06-09 Thread Tom Mulholland
While using Lattice I received the following error. Error in if (xx != 0) xx/10 else z/10 : argument is of length zero In addition: Warning messages: 1: is.na() applied to non-(list or vector) in: is.na(x) 2: is.na() applied to non-(list or vector) in: is.na(x) 3: no finite arguments to min;

Re: [R] SJAVA error

2004-06-09 Thread Vincent MUTEAUD
Hi, I have the same error than you two weeks ago. As I have not found how to do, I give up my example. Yunko, I am interesting if you found how to proceed. Thanks Vincent A 01:07 09/06/2004 +0900, =?iso-2022-jp?B?GyRCTHBMbhsoQiAbJEI9ZztSGyhC?= a écrit : Hi I'm trying to use SJava and I have

Re: [R] SJAVA error

2004-06-09 Thread Prof Brian Ripley
SJava is not part of R but rather of Omegahat, so this is the wrong list. But as a hint, I think you have not set R_HOME, or not set it correctly. But please consult the R posting guide for hints as to what information might be needed to help either of you, as neither or you have provided basic

Re: [R] Is there an R-version of rayplot

2004-06-09 Thread Martin Maechler
richard == richard kittler [EMAIL PROTECTED] on Tue, 8 Jun 2004 09:19:51 -0700 writes: richard I need to make plots similar to those produced by richard the s-plus rayplot function but can't seem to find richard it in R. These 'vector maps' plot a ray or vector richard at

[R] nominal data

2004-06-09 Thread Bianca
Hi, I am a new user of R. I have 2 series of nominal data (2 series of answers for the same question) and I want to calculate the correlation between these 2 series. I've tried to use the correlation function (corr, ...) but all are for numeric data... Does anyone know what function should I use

Re: [R] nominal data

2004-06-09 Thread Unung Istopo Hartanto
Hi, try ?cor but. I doubt it's your problem, and I think there's an another way to make correlation with nominal data, not used ?cor. But, i've been explore it yet. ex. : cbind(c(one,one,two,three,two,three),c(one,three,three,three,two,two)) - test test [,1][,2] [1,] one one [2,]

Re: [R] SJAVA error

2004-06-09 Thread Vincent MUTEAUD
Thanks I am using windows NT, R 1.9.0 and SJava 0.65 modified. I think my R_HOME is set correctly A 08:26 09/06/2004 +0100, Prof Brian Ripley a écrit : SJava is not part of R but rather of Omegahat, so this is the wrong list. But as a hint, I think you have not set R_HOME, or not set it

[R] GlmmPQL

2004-06-09 Thread O Tosas Auguet
Dear all, I have two questions concerning model simplification in GlmmPQL, for for random and fixed effects: 1. Fixed effects: I don't know if I can simply specify anova(model) and trust the table that comes up with the p value for each variable in the fixed effects formula. I have read that

[R] attach in R corr. to Spus one

2004-06-09 Thread Utsav Boobna
Hi, I am a newbie to R, just trying to switch from Splus. I am wondering to know if there is any such command in R as attach in Splus which can be used to attach the full directory. Also, any idea/referrence about, how to load c functions in R would be highly appreciated. Thanks, Utsav

Re: [R] SJAVA error

2004-06-09 Thread Wolski
Hi! For Windows and mac os X use of rJava may help. http://stats.math.uni-augsburg.de/iPlots/ From there also as it seems a working older version of SJave can be downloaded. Sincerely Eryk *** REPLY SEPARATOR *** On 6/9/2004 at 9:12 AM Vincent MUTEAUD wrote: Hi, I have the

Re: [R] SJAVA error

2004-06-09 Thread Prof Brian Ripley
Please do read what I said. This is *the wrong list*, so please use teh correct one. On Wed, 9 Jun 2004, Vincent MUTEAUD wrote: Thanks I am using windows NT, R 1.9.0 and SJava 0.65 modified. I think my R_HOME is set correctly It cannot be or the base package would be found. The relevant

Re: [R] GlmmPQL

2004-06-09 Thread Prof Brian Ripley
I think you need to read the references given in ?glmmPQL and its reference to understand what PQL actually does. If you don't know the theory behind a statistical method, you should try to understand it before trying to use it. The same applies to AIC. Even if AIC were computable, do you know

[R] inconsistency on p-value calculation of anova for quasi binomial

2004-06-09 Thread Ulrich Halekoh
Hej, providing the dispersion parameter estimate to the anova function for a quasibinomial fit results in two different ways to calculate the p-value for the same statistic. In the following example I test for the interaction effect. In the versions (a1 and a2) the p-value is based on the

[R] Getting Pr from Summary(lm)

2004-06-09 Thread David J. Netherway
Hello, I am trying to get the P values from the output of a summary for lm. lm - lm(y ~ age + sex) s - summary(lm) I thought that I might be able to get them using a combination of scan, grep and sub. But I got stuck on the first step - being able to process s as a text string. I could perhaps

Re: [R] attach in R corr. to Spus one

2004-06-09 Thread Petr Pikal
Hallo On 9 Jun 2004 at 2:06, Utsav Boobna wrote: Hi, I am a newbie to R, just trying to switch from Splus. I am wondering to know if there is any such command in R as attach in Splus which can be used to attach the full directory. I would recommend you to spent some time going through

Re: [R] Getting Pr from Summary(lm)

2004-06-09 Thread Peter Dalgaard
David J. Netherway [EMAIL PROTECTED] writes: Hello, I am trying to get the P values from the output of a summary for lm. lm - lm(y ~ age + sex) s - summary(lm) I thought that I might be able to get them using a combination of scan, grep and sub. But I got stuck on the first step -

Re: [R] Getting Pr from Summary(lm)

2004-06-09 Thread Prof Brian Ripley
coef(summary(lmfit))[, 4] See ?summary.lm, and note that summary() methods in general produce an object which the print method then prints. On Wed, 9 Jun 2004, David J. Netherway wrote: Hello, I am trying to get the P values from the output of a summary for lm. lm - lm(y ~ age + sex)

Re: [R] Getting Pr from Summary(lm)

2004-06-09 Thread Torsten Hothorn
Hello, I am trying to get the P values from the output of a summary for lm. lm - lm(y ~ age + sex) s - summary(lm) s$coefficients gives you a matrix with the P-values in the fourth column Torsten I thought that I might be able to get them using a combination of scan, grep and sub.

Re: [R] nominal data

2004-06-09 Thread Jonathan Baron
On 06/09/04 10:15, Bianca wrote: Hi, I am a new user of R. I have 2 series of nominal data (2 series of answers for the same question) and I want to calculate the correlation between these 2 series. I've tried to use the correlation function (corr, ...) but all are for numeric data... Does anyone

[R] testing effects of quantitative predictors on a categorical response variable

2004-06-09 Thread Avril Coghlan
Hello, I have a small statistics question, and as I'm quite new to statistics and R, I'm not sure if I'm doing things correctly. I am looking at two quantitative variables (x,y) that are correlated. When I divide the data set according to a categorical variable z, then x and y are more

Re: [R] Getting Pr from Summary(lm)

2004-06-09 Thread jferrer
Hi David, To see what is stored in an object, use names(). In your case you could try this: lm - lm(y ~ age + sex) s - summary(lm) names(s) [1] call terms residuals coefficients [5] aliased sigma dfr.squared [9] adj.r.squared fstatistic

Re: [R] testing effects of quantitative predictors on a categorical response variable

2004-06-09 Thread jferrer
Hi Avril, I'm not sure what you want to show. Do you want to know the effects of each variable? or just predict when you get z=A and when z=B? In the latter case, I think that, if x and y are in the same units, you could simply try w - y-x glm1 - glm(z ~ w,family=binomial(),trace=T) Hope it

Re: [R] poisson regression with robust error variance ('eyestudy')

2004-06-09 Thread Lutz Ph. Breitling
okay, so now the bootcov works fine. aren't the lower bootstrap variances just what Karla is talking about when she writes on the website describing the eyestudy that i was trying to redo in the first place: Using a Poisson model without robust error variances will result in a confidence

Re: [R] Is there an R-version of rayplot

2004-06-09 Thread joerg van den hoff
maybe this qd try helps? #=cut herer= vectorplot - function (field) { #input is a (N x 4 array) of N vectors: # field[,1:2] - x/y position of vectors # field[,3:4] - x/y componnent of vectors # plotted are the 2-D vectors attached to the specified

[R] market-basket analysis in R

2004-06-09 Thread Wayne Jones
Hi there fellow R-users, Does anyone know if there exists a package for associated rules data mining (market basket analysis) in R. I have tried searching CRAN but with no luck. Regards Wayne KSS Ltd Seventh Floor St James's Buildings 79 Oxford Street Manchester M1 6SS England

[R] Greek fonts

2004-06-09 Thread Johannes SCHNITZLER
Dear all, I need Greek text in titles, axis labels, etc. for graphs. The Graphs are produced in png-format - to use them later on in HTML pages. As the output has to be multilingual they are produced in a loop (for each language). The strings for the titles, labels ... are stored in a

RE: [R] market-basket analysis in R

2004-06-09 Thread RINNER Heinrich
Hi, I have asked a similar question quite a while ago, maybe you find the replys given in that thread useful: http://finzi.psych.upenn.edu/R/Rhelp02/archive/3977.html. Anyway, this is two years old now, so maybe there's some news in the meantime. Regards Heinrich. -Ursprüngliche

Re: [R] attach in R corr. to Spus one

2004-06-09 Thread Frank E Harrell Jr
Petr Pikal wrote: Hallo On 9 Jun 2004 at 2:06, Utsav Boobna wrote: Hi, I am a newbie to R, just trying to switch from Splus. I am wondering to know if there is any such command in R as attach in Splus which can be used to attach the full directory. I would recommend you to spent some time

RE: [R] more obvious contribution mechanism?

2004-06-09 Thread Liaw, Andy
There's no paypal button, but it isn't that much harder: Download the pdf form, fill it out and fax it back. How hard can that be? Andy From: Jason Turner On Wed, 2004-06-09 at 13:22, ivo welch wrote: can we put a how to donate money to R on the R webpage? perhaps with a paypal

RE: [R] more obvious contribution mechanism?

2004-06-09 Thread Peter Flom
It certainly isn't hard to do it this way, but a PayPal button is even easier. I think adding one would bring in extra money. Just my opinion, of course Peter Peter L. Flom, PhD Assistant Director, Statistics and Data Analysis Core Center for Drug Use and HIV Research National Development

Re: [R] e1071, R1.9.0, Solaris 2.9, should I be worried?

2004-06-09 Thread Orit Harel (Gdalyahu)
Hi How are you ? Did you get any answer regarding String literal converted to char* ? Thanks. Thanks and regards, Orit Harel * 972-9-7767513 * [EMAIL PROTECTED] The information contained in this message is proprietary of Amdocs, protected from disclosure, and may be privileged. The

Re: [R] market-basket analysis in R

2004-06-09 Thread Christian Schulz
I don't know any package!? but you could try: http://fuzzy.cs.uni-magdeburg.de/~borgelt/software.html http://www.cs.umb.edu/~laur/ARtool/ http://www.cs.waikato.ac.nz/ml/weka/ regards, Christian Am Mittwoch, 9. Juni 2004 13:27 schrieb Wayne Jones: Hi there fellow R-users, Does anyone know

[R] Re: more obvious contribution mechanism

2004-06-09 Thread ivo welch
hi: thanks everyone for pointing me to the contribution page.. [a] it was not obvious to me how to find this page. if you like contributions, please make it a bit more obvious. I would think a button on the home page would be a good idea. [b] any kind of trinket sale (a CD?) for money would

RE: [R] attach in R corr. to Spus one

2004-06-09 Thread Liaw, Andy
From: Frank E Harrell Jr Petr Pikal wrote: Hallo On 9 Jun 2004 at 2:06, Utsav Boobna wrote: Hi, I am a newbie to R, just trying to switch from Splus. I am wondering to know if there is any such command in R as attach in Splus which can be used to attach the full

RE: [R] attach in R corr. to Spus one

2004-06-09 Thread Liaw, Andy
Now that I acutally read the original question: If there are .RData files in other directories that you want to attach(), you can certainly do that. One difference between R and S-plus is that in S-plus you can attach() a directory at position 1, effectively making that the place that new

Re: [R] GLMM(..., family=binomial(link=cloglog))?

2004-06-09 Thread Göran Broström
On Tue, Jun 08, 2004 at 08:32:24AM -0700, Spencer Graves wrote: Hi, Doug: Thanks. I'll try the things you suggests. The observed proportions ranged from roughly 0.2 to 0.8 in 100 binomial random samples where sigma is at most 0.05. Jim Lindsey's glmm does Gauss-Hermite

RE: [R] Specifying xlevels in effects library

2004-06-09 Thread Liaw, Andy
Prof. Fox will be able to give the definitive answer, but from my reading of ?effect, xlevels refers to the values of the factor whose effect you're interested in, not the ones being `marginalized'. I believe you need to play with the `typical' argument. HTH, Andy From: David J. Netherway

Re: [R] SJAVA error

2004-06-09 Thread $BLpLn(B $B=g;R(B
I am using windows xp, R 1.9.0 and SJava 0.65 modified. (Bperhaps I set correctly. (B (Bcould you send me your compiled SJava package with the (Bmodified REmbed.c because in Windows i'm not able to recompile!!! (Bcould you provide me your solution. (B (B--- (BCLASSPATH (B

RE: [R] Specifying xlevels in effects library

2004-06-09 Thread Liaw, Andy
Here's an example: library(effects) age - round(rnorm(100, mean=30, sd=5)) sex - factor(sample(c(M,F), 100, replace=TRUE)) y - rnorm(100) fit - lm(y ~ age + sex) eff - effect(sex, fit, typical=function(...) 25) eff sex effect sex F M -0.23387685 0.07063834 Now we

RE: [R] fast mkChar

2004-06-09 Thread Paulo Nuin
Hello everyone This is my first message to the list and I believe the question I am including is a simple one. I have a matrix where I need to calculate ANOVA for the rows as the columns represent a different treatment. I would like to know if there is a command or a series of commans that I can

RE: [R] attach in R corr. to Spus one

2004-06-09 Thread Petr Pikal
Hallo Andy, Well, my original respond was to Utsav's question if it is possible to attach a directory in R as in S+. I personally prefer to start R session from a directory where I have all data and results from certain project (txt, doc, xls, png, pdf). Therefore I use only .RData and I

[R] Giovanna Jonalasinio è fuori ufficio- away from the office

2004-06-09 Thread Giovanna . Jonalasinio
Sarò fuori ufficio a partire dal 09/06/04 fino al 12/06/04 Sono al congresso SIS2004 a Bari. Non so se avrò accesso alla posta elettronica. I'll be away from the 09th to the 12th of June (SIS Conference in Bari) and I'm not sure I'll be able to read my mail [[alternative HTML

[R] Re: fighting with ps.options and xlim/ylim

2004-06-09 Thread ivo welch
Thanks again for all the messages. Is the 4% in par('usr') hardcoded? if so, may I suggest making this a user-changeable parameter for x and y axis? I looked at psfrag, and it seems like a great package. alas, I have switched to pdflatex, and pdffrag does not exist. :-( I also discovered

RE: [R] attach in R corr. to Spus one

2004-06-09 Thread Liaw, Andy
Hi Petr, I don't think I've said anything that is inconsistent with what you said, have I? The strategy that you mentioned is the one that is recommended, and I use it myself. However, life is not always that simple. There are times when the task at hand involves data/code/whatever that live

Re: [R] SJAVA error

2004-06-09 Thread Uwe Ligges
$BLpLn(B $B=g;R(B wrote: (B (B I am using windows xp, R 1.9.0 and SJava 0.65 modified. perhaps I set (B correctly. (B (B could you send me your compiled SJava package with the modified (B REmbed.c because in Windows i'm not able to recompile!!! (B could you provide me your solution.

Re: [R] Re: fighting with ps.options and xlim/ylim

2004-06-09 Thread Uwe Ligges
ivo welch wrote: Thanks again for all the messages. Is the 4% in par('usr') hardcoded? if so, may I suggest making this a user-changeable parameter for x and y axis? See ?par and its argumets xaxp, yaxp which can be set to i. I looked at psfrag, and it seems like a great package. alas, I have

[R] Dot chart question

2004-06-09 Thread Peter Flom
Running R 1.8.1 on a Windows machine In dotchart, I would like to shrink the labels on the tick marks (that is, the numbers) without shrinking anything else. I could not find this in either the Rhelp archives or in ?dotchart, which recmmends cex to avoid 'label overlap', but cex shrinks all the

Re: [R] Dot chart question

2004-06-09 Thread Chuck Cleland
See cex.axis and other cex.* arguments to par() in ?par. For example: data(VADeaths) dotchart(VADeaths, main = Death Rates in Virginia - 1940) par(cex.axis = .6) dotchart(VADeaths, main = Death Rates in Virginia - 1940) hope this helps, Chuck Cleland Peter Flom wrote: Running R

Re: [R] GLMM(..., family=binomial(link=cloglog))?

2004-06-09 Thread Spencer Graves
Hi, Go"ran: (B (BThanks for the analysis. Unfortunately, it still leaves me with 2 (Bproblems. First, I'm dealing with extremely small defect rates involving (Bthousands and millions of Bernoulli trials, so creating bigDF would (Brequire computers with much more memory and processing speed

[R] X-12-ARIMA

2004-06-09 Thread Paul Gilbert
(subject changed from [R] a doubt ) The fortran source code for X-12-ARIMA seems to be available at http://www.census.gov/srd/www/x12a/x12down_unix.html. Incorporating this into a standard R package would be a much better approach than trying to hack an OS specific version with system calls to

Re: [R] Dot chart question

2004-06-09 Thread Peter Flom
Thanks, this worked. I had tried fx.dotchart2(freqrelig[order(freqrelig)], label = levels(religfact)[order(freqrelig)], main = 'Religions among young adults in Bushwick \n (log scale)', xlab = 'Frequency', log = 'x', cex.axis = .7) which did nothing... Peter Chuck Cleland [EMAIL

[R] Anova question

2004-06-09 Thread Paulo Nuin
Hello everyone This is my first message to the list and I believe the question I am including is a simple one. I have a matrix where I need to calculate ANOVA for the rows as the columns represent a different treatment. I would like to know if there is a command or a series of commans that I can

[R] robust correlation in R

2004-06-09 Thread Xavier . Abulker
Dear R user group, I'm looking for a robust mesure of correlation in R. I found a very interesting article by Dr Rich Herrington on http://www.unt.edu/benchmarks/archives/2001/december01/rss.htm and I'd like to implement exaclty this method but my problem is that everything is here developped

Reshaping matrix (was: [R] fast mkChar)

2004-06-09 Thread Duncan Murdoch
On 09 Jun 2004 09:52:27 -0400, Paulo Nuin [EMAIL PROTECTED] wrote : Hello everyone This is my first message to the list and I believe the question I am including is a simple one. I have a matrix where I need to calculate ANOVA for the rows as the columns represent a different treatment. I would

[R] robust correlation in R

2004-06-09 Thread Xavier . Abulker
Dear R user group, I'm looking for a robust mesure of correlation in R. I found a very interesting article by Dr Rich Herrington on http://www.unt.edu/benchmarks/archives/2001/december01/rss.htm and I'd like to implement exaclty this method but my problem is that everything is here

[R] ask for data manipulation!

2004-06-09 Thread vinkwai wong
i have such a data: V1 V2 V3 V4 abc e afd NA bde f hdNA NA adf e (V1,V2,V3,V4 can be any one of a, b, c, d, e, f, g, h. that is to say ,it has 8 possible values.) and i want to change it the the following form: abcdef

Re: [R] X-12-ARIMA

2004-06-09 Thread Dirk Eddelbuettel
On Wed, Jun 09, 2004 at 10:41:58AM -0400, Paul Gilbert wrote: (subject changed from [R] a doubt ) The fortran source code for X-12-ARIMA seems to be available at http://www.census.gov/srd/www/x12a/x12down_unix.html. Incorporating this into a standard R package would be a much better

Wrong question [Wasn't: Re: [R] fast mkChar]

2004-06-09 Thread Matej Cepl
On Wednesday 09 of June 2004 09:52, you wrote: This is my first message to the list and I believe the question I am including is a simple one. http://www.r-project.org/posting-guide.html -- Matej Cepl, http://www.ceplovi.cz/matej GPG Finger: 89EF 4BC6 288A BF43 1BAB 25C3 E09F EF25 D964 84AC

[R] Multiple regression

2004-06-09 Thread Manoj - Hachibushu Capital
Hi, I am trying to do multiple regression on a set of data using backward stepwise regressionhowever backward stepwise regression is critised for overfitting data. To actually observe the bias and to come up with a better method to use..Could you all stats experts kindly give me

RE: [R] Is there an R-version of rayplot

2004-06-09 Thread richard . kittler
Thanks for this example. I have since found the 'arrows' function which is not as comprehensive as rayplot but seems to work well for basic vector field plots, e.g. y - rep(1:5,rep(5,5)) x - rep(1:5,5) x1 - x-sqrt(x*x + y*y)/25 y1 - y-2*x/125 plot(x,y,,main='Sample Vector

Re: [R] Multiple regression

2004-06-09 Thread Dimitris Rizopoulos
Have a look at Prof. Harrell's book (Chapters 4 and 5) @Book{harrell:01, author= {F. E. Harrell, Jr.}, title = {Regression Modeling Strategies: With Applications to Linear Models, Logistic Regression and Survival Analysis}, year = {2001}, address = {New

Re: [R] market-basket analysis in R

2004-06-09 Thread Christian Schulz
Hi, when i want start more than one year ago with a fuzzy-association rule algorithm to learn programming in r i have to work in business and it's too bad that here until now nobody need a tool for market basket analysis in r :-). regards,Christian Am Mittwoch, 9. Juni 2004 15:38 schrieb

Re: [R] Dot chart question

2004-06-09 Thread Frank E Harrell Jr
Peter Flom wrote: Running R 1.8.1 on a Windows machine In dotchart, I would like to shrink the labels on the tick marks (that is, the numbers) without shrinking anything else. I could not find this in either the Rhelp archives or in ?dotchart, which recmmends cex to avoid 'label overlap', but

[R] About dll from c++ routine

2004-06-09 Thread Rui
Hi folks, My system is Windows98 + R1.9.0. The path for my system is c:\perl\bin; c:\mingw\bin; c:\rtools; c:\windows; c:\windows\command; c:\rw1090\bin. I created three files followed the examples in “Writing R extensions” in the directory c:\temp: // X.hh class X { public: X (); ~X (); };

RE: [R] Comparing two pairs of non-normal datasets in R?

2004-06-09 Thread Lutz Prechelt
Here are the boxplots if that helps: http://www.ps.masny.dk/guests/misc/A1.png http://www.ps.masny.dk/guests/misc/A2.png http://www.ps.masny.dk/guests/misc/C1.png http://www.ps.masny.dk/guests/misc/C2.png Here is how I would do it: It looks like your distributions can be characterized by

Re: [R] Re: fighting with ps.options and xlim/ylim

2004-06-09 Thread Marc Schwartz
On Wed, 2004-06-09 at 09:30, Uwe Ligges wrote: ivo welch wrote: Thanks again for all the messages. Is the 4% in par('usr') hardcoded? if so, may I suggest making this a user-changeable parameter for x and y axis? See ?par and its argumets xaxp, yaxp which can be set to i.

Re: [R] Anova question

2004-06-09 Thread Adaikalavan Ramasamy
On Wed, 2004-06-09 at 15:36, Paulo Nuin wrote: Hello everyone This is my first message to the list and I believe the question I am including is a simple one. I have a matrix where I need to calculate ANOVA for the rows as the columns represent a different treatment. I would like to know

Re: [R] About dll from c++ routine

2004-06-09 Thread Prof Brian Ripley
On Wed, 9 Jun 2004, Rui wrote: Hi folks, My system is Windows98 + R1.9.0. The path for my system is c:\perl\bin; c:\mingw\bin; c:\rtools; c:\windows; c:\windows\command; c:\rw1090\bin. I created three files followed the examples in “Writing R extensions” in the directory c:\temp: //

[R] lsoda with arbitrary zero thresholds

2004-06-09 Thread Martin Henry H. Stevens
using R 2.0.0 I am trying to do some population modeling with lsoda, where I set arbitrary zero population sizes when values get close to zero, but am having no luck. As an example of what I have tried, I use code below from the help page on lsoda in which I include my modification bordered by

Re: [R] Specifying xlevels in effects library

2004-06-09 Thread John Fox
Dear David and Andy, First, my apologies for somehow missing the original message. Andy has it right: xlevels refers to the predictors in the term corresponding to the effect to be computed. His solution is clever -- I wouldn't have thought of it -- but it could get you into trouble in a more

[R] question related to S-Plus

2004-06-09 Thread Rick Picard
Dear r-help, Having used S-Plus for many years, it has been suggested to me that I could benefit from learning R. A question to this end, though. Would all of my existing S-Plus arrays, functions, and so on have to be re-created from scratch in R, or is there a way to copy them into the

Re: [R] lsoda with arbitrary zero thresholds (with psuedo-solution)

2004-06-09 Thread Martin Henry H. Stevens
I have a new and less distressing, but potentially more interesting, problem. I realized the major flaw my old solution and now have a solution that kind of works but is rather inelegant and I think may be problematic in difficult systems. Borrowing from the lsoda example again I once again

[R] [R-pkgs] gregmisc 1.11.2 including read.xls()

2004-06-09 Thread Warnes, Gregory R
Release 1.11.2 of the gregmisc() package of functions is now available on CRAN at http://cran.r-project.org/src/contrib/Descriptions/gregmisc.html for both Unix and Windows systems. The most notable enhancement provided by this release is: - read.xls(), a function to read Microsoft Excel files

[R] Re: R equivalent of Splus rowVars function

2004-06-09 Thread David Brahm
Mark Leeds [EMAIL PROTECTED] wrote (to S-News): does anyone know the R equivalent of the SPlus rowVars function ? Andy Liaw [EMAIL PROTECTED] replied: More seriously, I seem to recall David Brahms at one time had created an R package with these dimensional summary statistics, using C code.

[R] Building package on Windows: No rule to make target '-llapack'

2004-06-09 Thread Zhu Wang
Dear all, I have a problem to build a package on Windows XP while there is no problem on Linux. The Makefile is something like: ### LIBNAME=cts PKG_LIBS = $(LAPACK_LIBS) $(BLAS_LIBS) $(FLIBS) OBJS=file1.o ... file20.o -llapack -lblas $(LIBNAME)$(SHLIB_EXT): $(OBJS) $(SHLIB_LD)

Re: [R] Building package on Windows: No rule to make target '-llapack'

2004-06-09 Thread A.J. Rossini
Zhu Wang [EMAIL PROTECTED] writes: Now I think maybe two problems: one is that maybe I do not have Lapack and Blas installed I think this definitely needs to be solved. -- [EMAIL PROTECTED]http://www.analytics.washington.edu/ Biomedical and Health Informatics University of

RE: [R] fast mkChar

2004-06-09 Thread Vadim Ogranovich
Thank you for the lead, Peter. It may be useful for other packages I write. As to the strings, I think I have to take what is already there. I agree that strings would be better managed in malloc-style fashion (probably with reference counter) and not by gc(). However I don't want to have a

[R] how to initialize random seed properly ?

2004-06-09 Thread ryszard . czerminski
I want to start R processes on multiple processors from single shell script and I want all of them to have different random seeds. One way of doing this is sleep 2 # (with 'sleep 1' I am often getting the same number) ... set.seed(unclass(Sys.time())) Is there a

[R] direct data frame entry

2004-06-09 Thread ivo welch
hi: I searched the last 2 hours for a way to enter a data frame directly in my program. (I know how to read from a file.) that is, I would like to say something like d - this.is.a.data.frame( c(obs1name, 0.2, 0.3), c(obs2name, 0.4, 1.0),

RE: [R] direct data frame entry

2004-06-09 Thread Liaw, Andy
?data.frame says: Usage: data.frame(..., row.names = NULL, check.rows = FALSE, check.names = TRUE) Arguments: ...: these arguments are of either the form 'value' or 'tag=value'. Component names are created based on the tag (if present) or the deparsed argument

Re: [R] direct data frame entry

2004-06-09 Thread Tony Plate
easy to do it by column: d - data.frame(name=c(obs1name,obs2name,obs3name),val1=c(0.2,0.4,0.6),val2=c(0.3,1.0,2.0),row.names=c(r1,r2,r3)) d name val1 val2 r1 obs1name 0.2 0.3 r2 obs2name 0.4 1.0 r3 obs3name 0.6 2.0 (when you do it by row, you get the numbers as factors because

Re: [R] how to initialize random seed properly ?

2004-06-09 Thread A.J. Rossini
For reproducibility, you probably want to specify the starting seeds. I'd worry about using systime. (for selecting the seeds, you might consider random draws from a uniform). best, -tony [EMAIL PROTECTED] writes: I want to start R processes on multiple processors from single shell script

[R] moving data and output?

2004-06-09 Thread Martin Klaffenboeck
Hello, I have a few questions now: 1. How can I move data the following way: I have 2 variables: one two 1 5 ^ 3 4 | 1 3 | 4 4 | Now I want to move the two one arround (sorry I don't know how to say that in english). That means: I want to move the first item at

Re: [R] direct data frame entry

2004-06-09 Thread ivo welch
thank you, chaps. ok, so this is not as straightforward as I had thought. perhaps the read.table() function should have the ability to read inline (terminated, e.g., by two newlines, or a usersettable string), rather than just from a file. this would be a nice feature. regards, /iaw

Re: [R] Comparing two pairs of non-normal datasets in R?

2004-06-09 Thread Peter Sebastian Masny
On Wednesday 09 June 2004 08:58 am, Lutz Prechelt wrote: ... Does that make sense? Lutz Thanks, everyone, very much for your replies. I think I'm getting a little out of my league with quantile regession estimates and bootstrap tests of distribution parameters. It makes sense, but I'm

Re: [R] direct data frame entry

2004-06-09 Thread Wolski
Hi Ivo! https://www.stat.math.ethz.ch/pipermail/r-help/2004-June/050601.html Sincerely Eryk *** REPLY SEPARATOR *** On 6/9/2004 at 3:29 PM ivo welch wrote: hi: I searched the last 2 hours for a way to enter a data frame directly in my program. (I know how to read from a

[R] Two-dimensional Kolmogorov-Smirnov test

2004-06-09 Thread Rich Grenyer
Hi - is the 2-D Kolmogorov-Smirnov test of Fasano Francheschini (1987) implemented in any of the R-packages at present? The algorithm and code exist in the Numerical Recipes series (e.g. Press et al. 1992) but I thought I should check before playing with .C() and dyn.load() for the first

Re: [R] direct data frame entry

2004-06-09 Thread Gabor Grothendieck
ivo welch ivo.welch at yale.edu writes: thank you, chaps. ok, so this is not as straightforward as I had thought. perhaps the read.table() function should have the ability to read inline (terminated, e.g., by two newlines, or a usersettable string), rather than just from a file. this

Re: [R] Help with a Lattice plot that fails with an empty unique combination

2004-06-09 Thread Deepayan Sarkar
On Wednesday 09 June 2004 01:58, Tom Mulholland wrote: While using Lattice I received the following error. Error in if (xx != 0) xx/10 else z/10 : argument is of length zero In addition: Warning messages: 1: is.na() applied to non-(list or vector) in: is.na(x) 2: is.na() applied to non-(list

Re: [R] bar plot patterns

2004-06-09 Thread Paul Murrell
Hi Osman wrote: Thank you for your answer. I have about 7 stacks it is not very appealing to have just stripes with changing angles. I was wondering if there is a way to have varying patters in black and white. R graphics is inherently vector-based; there is no native support for bitmaps so

Re: [R] Specifying xlevels in effects library

2004-06-09 Thread David J. Netherway
Andy and John, I looked at typical when xlevels did not work but when I saw that it was a function I went no further. Setting the function to a constant was a good idea. John's method seems to require that I change the model: eff -effect(sex*age,mod,xlevel=(Age=c(120,120))) Error in

Re: [R] Getting Pr from Summary(lm)

2004-06-09 Thread David J. Netherway
Thanks for all the relies. I recently discovered names and applied it to lm objects but did not think to apply it to the summary object. Cheers, David __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do

RE: [R] Specifying xlevels in effects library

2004-06-09 Thread John Fox
Dear David, You don't have to change the model, but you do have to name the variable age (or Age) consistently, and give a list as the xlevels argument. Hence, if the variable names are really lower-case: eff - effect(sex*age, mod, xlevels=list(age=c(120,120))). Regards, John -Original

RE: [R] Specifying xlevels in effects library [2]

2004-06-09 Thread John Fox
Dear David, An addendum: Looking at my original posting, I see that failing to specify a list as the xlevels argument was my error. Sorry, John -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of David J. Netherway Sent: Wednesday, June 09, 2004

Re: [R] X-12-ARIMA

2004-06-09 Thread Alistair Gray
Dear All, I've used the X-12-ARIMA or its earlier versions from S+ and R under both Unix and Windows platforms for many years using the klugey approach of calling an executable using in R the system function. I've found this serviceable for the following reasons. 1) Paul Gilbert's hunch is

Re: [R] X-12-ARIMA

2004-06-09 Thread Dirk Eddelbuettel
Alistair, That look like a very good way of getting things done. I had at one point in the past thought about undertaking a similar approach of dealing with the input/output files to x12arima from R. However, my job priorities changed and I didn't have the need for x12arima any more. For the

Re: [R] moving data and output?

2004-06-09 Thread Jason Turner
Hello, I have a few questions now: Yes. Since the mail is archived to help other people, in future please send a small mail for each question, with a descriptive subject line for each. This makes it easier for everyone. 1. How can I move data the following way: I have 2 variables: one