[R] RE: do0i4grjj40j09gjijgpd

2005-04-06 Thread domino . admin
Privitak je zaraen virusom. Kontaktirajte sistem administratora. An attachment is infected by virus. Contact administrator. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: [R] French Curve

2005-04-06 Thread Martin Maechler
Michael == Michael A Miller [EMAIL PROTECTED] on Tue, 05 Apr 2005 10:28:21 -0500 writes: dream == dream home [EMAIL PROTECTED] writes: Does it sound like spline work do the job? It would be hard to persuave him to use some modern math technique but he did ask me to help him

Re: [R] two methods for regression, two different results

2005-04-06 Thread Jari Oksanen
On Tue, 2005-04-05 at 22:54 -0400, John Sorkin wrote: Please forgive a straight stats question, and the informal notation. let us say we wish to perform a liner regression: y=b0 + b1*x + b2*z There are two ways this can be done, the usual way, as a single regression, fit1-lm(y~x+z)

Re: [R] error messages on R CMD check

2005-04-06 Thread Johannes Hüsing
Hallo Uwe, Uwe Ligges schrieb: [...] In Johannes' case, the problem is different, because the error message is not that clear. Johannes, can you install and load the package? Is the DESCRIPTION file correct? If so, you might want to send the package in a private message... I have found the

RE: [R] Help with three-way anova

2005-04-06 Thread michael watson \(IAH-C\)
OK, so I tried using lm() instead of aov() and they give similar results: My.aov - aov(IL.4 ~ Infected + Vaccinated + Lesions, data) My.lm - lm(IL.4 ~ Infected + Vaccinated + Lesions, data) If I do summary(My.lm) and summary(My.aov), I get similar results, but not identical. If I do

Re: [R] Fitdistr and likelihood

2005-04-06 Thread Martin Maechler
{BCC'ed to VR's maintainer} Carsten == Carsten Steinhoff [EMAIL PROTECTED] on Tue, 5 Apr 2005 17:31:04 +0200 writes: Carsten Hi all, I'm using the function fitdistr (library Carsten MASS) to fit a distribution to given data. What I Carsten have to do further, is getting

Re: [R] step error

2005-04-06 Thread Dimitris Rizopoulos
an indirect solution is to use the following: pro - function(indep, dep){ d - data.frame(indep) form - formula(lm(dep~., data=d)) assign(d, d, envir=.GlobalEnv) res - step(lm(dep~X1, data=d), scope=form, trace=0, direction=f) rm(d, envir=.GlobalEnv) res } # m -

RE: [R] Help with three-way anova

2005-04-06 Thread Federico Calboli
On Wed, 2005-04-06 at 09:11 +0100, michael watson (IAH-C) wrote: OK, so I tried using lm() instead of aov() and they give similar results: My.aov - aov(IL.4 ~ Infected + Vaccinated + Lesions, data) My.lm - lm(IL.4 ~ Infected + Vaccinated + Lesions, data) Incidentally, if you want

RE: [R] Help with three-way anova

2005-04-06 Thread michael watson \(IAH-C\)
OK, now I am lost. I went from using aov(), which I fully understand, to lm() which I probably don't. I didn't specify a contrasts matrix in my call to lm() Basically I want to find out if Infected/Uninfected affects the level of IL.4, and if Vaccinated/Unvaccinated affects the level of

[R] conditional selection with Factors

2005-04-06 Thread H J Gerber
Hi All, I want to select a level of a factor variable from a dataset. I have the folowing data: dataset: use factor: month (use$month, levels=February, July) use[1:5,] month registration use department size 1 February KKG151GP Y Safety 1.6 To select February I tried:

Re: [R] R can not show plots (in Mac OS X terminal)

2005-04-06 Thread David Ruau
Hi, You should use X11. It doesn't work in Terminal. You can use the basic Xterm in X11 or like I do Aterm. David Ruau On Apr 5, 2005, at 20:12, Minyu Chen wrote: Dear all: I am a newbie in Mac. Just installed R and found R did not react on my command plot (I use command line in terminal). It did

Re: [R] Using kmeans given cluster centroids and data with NAs

2005-04-06 Thread Ales Ziberna
Hello! I would suggest using some form of imputations, such as MICE package (http://web.inter.nl.net/users/S.van.Buuren/mi/hmtl/mice.htm) or similar (I heard that this can be also done with aregImpute function in the Hmisc package, although I have not tried it) to fill in the NA's. Then you can

[R] nls.control

2005-04-06 Thread Anthony Landrevie
Hello everyone, I'm trying to test the accurracy of R on the Eckerle4 dataset from NIST and I don't understand how the control option of the nls function works. I tought nls(...) was equivalent to nls(...control=nls.control()) i.e nls.control() was the default value of control, but here is the

Re: [R] conditional selection with Factors

2005-04-06 Thread Petr Pikal
Hi On 6 Apr 2005 at 11:34, H J Gerber wrote: Hi All, I want to select a level of a factor variable from a dataset. I have the folowing data: dataset: use factor: month (use$month, levels=February, July) use[1:5,] month registration use department size 1 February KKG151GP Y

RE: [R] Help with three-way anova

2005-04-06 Thread John Fox
Dear Mick, For a three-way ANOVA, the difference between aov() and lm() is mostly in the print and summary methods -- aov() calls lm() but in its summary prints an ANOVA table rather than coefficient estimates, etc. You can get the same ANOVA table from the object returned by lm via the anova()

[R] Precision

2005-04-06 Thread Josef Eschgfaeller
How precise is R numerically? For example I wrote the following function for calculating the volume of the ball inscribed in the unit cube in m dimensions. In order to see what happens in 40 dimensions, I created an output of 24 digits. But how many are precise? Thanks Josef Eschgfäller Ferrara

RE: [R] Help with three-way anova

2005-04-06 Thread michael watson \(IAH-C\)
Hi John Thanks for your help, that was a very clear answer. It looks as though, due to my design, the best way forward is: contrasts(il4$Infected) [,1] I-1 UI1 contrasts(il4$Vaccinated) [,1] UV -1 V 1 summary(lm(IL.4 ~ Infected * Vaccinated, il4)) Thanks Mick

Re: [R] nls.control

2005-04-06 Thread Martin Maechler
Anthony == Anthony Landrevie [EMAIL PROTECTED] on Wed, 6 Apr 2005 02:54:50 -0700 (PDT) writes: Anthony Hello everyone, I'm trying to test the accurracy of Anthony R on the Eckerle4 dataset from NIST Do you know that there's an R package 'NISTnls' (on CRAN) exactly for this

Re: [R] French Curve

2005-04-06 Thread roger koenker
On Apr 6, 2005, at 1:48 AM, Martin Maechler wrote: Median filtering aka running medians has one distinctive advantage {over smooth.spline() or other so called linear smoothers}: It is robust i.e. not distorted by gross outliers. Running medians is implemented in runmed() {standard stats

Re: [R] R can not show plots (in Mac OS X terminal)

2005-04-06 Thread Tony Han Bao
Hi On 6 Apr 2005, at 10:47, David Ruau wrote: Hi, You should use X11. It doesn't work in Terminal. To make Apple's Terminal use X11 first one should set the DISPLAY environment variable if you are using bash, put the following line in .bash_profile [[ -z $DISPLAY ]] export DISPLAY=:0.0 You can

[R] how to estimate Type I, Type III SS

2005-04-06 Thread Kum-Hoe Hwang
Howdy, R gurus I 'd like to know hwo to calculate or estimate SS of Type I and Type III in ANOVA or other anaysis in R. Thanks, -- Kum-Hoe Hwang, Ph.D. Kyonggi Research Institute, Korea (ROK) (Urban Planning and GIS) Phone : 82-31-250-3283 Email : [EMAIL PROTECTED] [[alternative

Re: [R] Precision

2005-04-06 Thread Jan T. Kim
On Wed, Apr 06, 2005 at 02:00:58PM +0200, Josef Eschgfaeller wrote: How precise is R numerically? For example I wrote the following function for calculating the volume of the ball inscribed in the unit cube in m dimensions. In order to see what happens in 40 dimensions, I created an output

Re: [R] Precision

2005-04-06 Thread Duncan Murdoch
Josef Eschgfaeller wrote: How precise is R numerically? For example I wrote the following function for calculating the volume of the ball inscribed in the unit cube in m dimensions. In order to see what happens in 40 dimensions, I created an output of 24 digits. But how many are precise? For most

Re: [R] Precision

2005-04-06 Thread Duncan Murdoch
Duncan Murdoch wrote: For most floating point operations R uses double precision, which gives about 18-19 significant digit precision. Leading zeros don't count. Oops, Jan Kim is right: double precision is only 15-16 digit precision. Sorry. Duncan Murdoch

Re: [R] how to estimate Type I, Type III SS

2005-04-06 Thread Federico Calboli
On Wed, 2005-04-06 at 21:40 +0900, Kum-Hoe Hwang wrote: Howdy, R gurus I 'd like to know hwo to calculate or estimate SS of Type I and Type III in ANOVA or other anaysis in R. Thanks, If memory seves me well, try Anova in the car package F -- Federico C. F. Calboli Department of

RE: [R] how to estimate Type I, Type III SS

2005-04-06 Thread John Fox
Dear Federico and Kum-Hoe, The Anova function in the car package will compute Type II or Type III tests (with the former as the default); anova() computes Type I (sequential) tests. Be careful with the contrast coding for Type III tests. Regards, John John Fox

[R] Courses- May 2005***R/S-plus Programming Techniques

2005-04-06 Thread elvis
XLSolutions Corporation (www.xlsolutions-corp.com) is proud to announce our May-2005 2-day R/S-plus Fundamentals and Programming Techniques at 4 US locations: www.xlsolutions-corp.com/training.htm Philadelphia --- May 19-20 Raleigh

[R] Introduce a new function in a package?

2005-04-06 Thread Luis Ridao Cruz
R-help, Sometimes I define functions I wish to have in any R session. The obvious thing to do is copy-paste the code The thing is that sometimes I don't know where I have the function code. My question is if somehow I could define a function and introduce it (let's say 'base' package ) so that

Re: [R] Introduce a new function in a package?

2005-04-06 Thread Roger D. Peng
I think the usual way is to create an R package for yourself and load it when you need it for whatever project. -roger Luis Ridao Cruz wrote: R-help, Sometimes I define functions I wish to have in any R session. The obvious thing to do is copy-paste the code The thing is that sometimes I don't

Re: [R] Introduce a new function in a package?

2005-04-06 Thread Dimitris Rizopoulos
or you could create a package containing all these functions and edit .Rprofile to load it at start-up (see also ?.Startup). Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven,

Re: [R] Introduce a new function in a package?

2005-04-06 Thread TEMPL Matthias
See at http://cran.r-project.org/doc/manuals/R-intro.pdf at page 54. .First() can help you. Or create an own package (see http://cran.r-project.org/doc/manuals/exts.pdf ) and load the package, when needed. Best, Matthias -Ursprüngliche Nachricht- Von: [EMAIL PROTECTED]

RE: [R] Introduce a new function in a package?

2005-04-06 Thread Liaw, Andy
Another approach, if making a package is a bit more than what you want to do, is: 1. Save those functions/objects in an image using save(). 2. attach() that image every time you start R. There are a few ways that you can do #2 above. See ?Startup. Andy From: Roger D. Peng I think the

[R] par(mfcol=2, mfrow=3) equivalent for trellis

2005-04-06 Thread Dieter Menne
Dear friends of lattice, I know how to position trellis plots with print(...,split,more=T) or (...position). Sometimes I wish I had something like the old par(mfcol=2, mfrow=3) mechanism, where the next free viewport is automatically chosen. I tried fiddling with grid-viewports, but could not

Re: [R] Introduce a new function in a package?

2005-04-06 Thread Jan T. Kim
On Wed, Apr 06, 2005 at 09:57:00AM -0400, Roger D. Peng wrote: I think the usual way is to create an R package for yourself and load it when you need it for whatever project. -roger Alternatively, one can also write the function in question into one's ~/.Rprofile; then, it's automatically

Re: [R] par(mfcol=2, mfrow=3) equivalent for trellis

2005-04-06 Thread Dimitris Rizopoulos
probably you want to use the layout argument, see its description in ?xyplot(). 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] the number of cluster

2005-04-06 Thread msck9
Are there any package that is using v-fold cross-validation algorithm to test the number the clusters? Thanks, Ming On Sat, Mar 19, 2005 at 12:44:42PM +0100, Uwe Ligges wrote: XP Sun wrote: hi, all, how to decide the number of cluster before you use kmeans and hclust? thank

Re: [R] Introduce a new function in a package?

2005-04-06 Thread Don MacQueen
Expressions in .Rprofile are executed *before* any previously saved global environment is loaded (i.e., before the .RData file in the current working directory is loaded, causing the message [Previously saved workspace restored] to a appear). If you define a function in .Rprofile, and then

[R] Error in hist.default(A) : `x' must be numeric

2005-04-06 Thread Mag. Ferri Leberl
Dear everybody! I have load a list A of numbers and want a histogram to be drawn. on hist(Y) the Machine returns: Error in hist.default(A) : `x' must be numeric I found out, that the list is of type data.frame. Y-as.numeric(Y) returns Error in as.double.default(A) : (list) object cannot be coerced

Re: [R] R can not show plots (in Mac OS X terminal)

2005-04-06 Thread Thomas Lumley
On Wed, 6 Apr 2005, David Ruau wrote: Hi, You should use X11. It doesn't work in Terminal. You can use the basic Xterm in X11 or like I do Aterm. This is not actually true. It does work in Terminal, you just have to specify the DISPLAY, either in the shell before entering R % setenv DISPLAY :0

Re: [R] par(mfcol=2, mfrow=3) equivalent for trellis

2005-04-06 Thread Deepayan Sarkar
On Wednesday 06 April 2005 09:06, Dieter Menne wrote: Dear friends of lattice, I know how to position trellis plots with print(...,split,more=T) or (...position). Sometimes I wish I had something like the old par(mfcol=2, mfrow=3) mechanism, where the next free viewport is automatically

Re: [R] Error in hist.default(A) : `x' must be numeric

2005-04-06 Thread Uwe Ligges
Mag. Ferri Leberl wrote: Dear everybody! I have load a list A of numbers and want a histogram to be drawn. on hist(Y) the Machine returns: Error in hist.default(A) : `x' must be numeric I found out, that the list is of type data.frame. Y-as.numeric(Y) returns Error in as.double.default(A) : (list)

Re: [R] Error in hist.default(A) : `x' must be numeric

2005-04-06 Thread Sundar Dorai-Raj
Mag. Ferri Leberl wrote on 4/6/2005 9:57 AM: Dear everybody! I have load a list A of numbers and want a histogram to be drawn. on hist(Y) the Machine returns: Error in hist.default(A) : `x' must be numeric I found out, that the list is of type data.frame. Y-as.numeric(Y) returns Error in

Re: [R] R can not show plots (in Mac OS X terminal)

2005-04-06 Thread Don MacQueen
Did you install R from source code, or did you install the binary? If you installed the binary, then you can start R by double-clicking on the R application icon. Then your default graphics device will not require X windows, and will be fully interactive (in the R sense). If you installed from

Re: [R] Introduce a new function in a package?

2005-04-06 Thread Gabor Grothendieck
Some other advantages of making your own package are: - you can use help.search to search for your own functions even if you don't load the package - if you can't even remember where your functions are (and I often can't) then you may not remember what they do either and packaging them

Re: [R] Introduce a new function in a package?

2005-04-06 Thread Josef Eschgfaeller
~/.Rprofile You could also write in .Rprofile soemthing like this: for (x in dir(Mylibrary,full.names=T,recursive=T)) source(x) where Mylibrary is a directory which contains your functions without making a package. Josef Eschgfäller__

RE: [R] Error in hist.default(A) : `x' must be numeric

2005-04-06 Thread Liaw, Andy
From: Uwe Ligges Mag. Ferri Leberl wrote: Dear everybody! I have load a list A of numbers and want a histogram to be drawn. on hist(Y) the Machine returns: Error in hist.default(A) : `x' must be numeric I found out, that the list is of type data.frame. Y-as.numeric(Y) returns

Re: [R] R can not show plots (in Mac OS X terminal)

2005-04-06 Thread Minyu Chen
Thank you very much. This is very informative and I already save it for future reference. Now I got the double clicking icon (quite mysteriously, since I tried several ways recommended by others, so I don't know which one make it works). Thanks, Minyu Chen On 6 Apr 2005, at 16:11, Don MacQueen

[R] off-topic question: Latex and R in industries

2005-04-06 Thread Wensui Liu
Latex and R are really cool stuff. I am just wondering how they are used in industry. But based on my own experience, very rare. Why? How about the opinion of other listers? Thanks. __ R-help@stat.math.ethz.ch mailing list

Re: [R] Introduce a new function in a package?

2005-04-06 Thread roger bos
I tried making a package on windows and got a make error, so I was happy I was able to get source(mystuff.R) to work in .First(). Since my utility functions are pretty simple and few in number, this is good enough for me for now. But I got a curious error. I can submit the command

[R] newman-keuls

2005-04-06 Thread Faouzi LYAZRHI
Bonjour Je cherche à faire un test de newman-keuls sous R. Est-ce que quelqu'un peut m'aider ? En vous remerciant d'avance fawtzy __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: [R] Introduce a new function in a package?

2005-04-06 Thread Prof Brian Ripley
Better to attach a new environment and source() the files into that. local({ for (.x in dir(Mylibrary,full.names=T,recursive=T)) source(.x, local=T); rm(.x)}, env = attach(NULL, name=myenv)) for all the reasons why a package() is a good idea. Leaving an object `x' around from .Rprofile is

Re: [R] off-topic question: Latex and R in industries

2005-04-06 Thread roger koenker
my favorite answer to this question is because there is no one to sue. url:www.econ.uiuc.edu/~rogerRoger Koenker email [EMAIL PROTECTED] Department of Economics vox:217-333-4558University of Illinois fax:217-244-6678

Re: [R] Introduce a new function in a package?

2005-04-06 Thread Prof Brian Ripley
On Wed, 6 Apr 2005, roger bos wrote: I tried making a package on windows and got a make error, so I was happy I was able to get source(mystuff.R) to work in .First(). Since my utility functions are pretty simple and few in number, this is good enough for me for now. But I got a curious error. I

Re: [R] Introduce a new function in a package?

2005-04-06 Thread Josef Eschgfaeller
Leaving an object `x' around from .Rprofile is not at all a good idea, Actually I thought to put it inside a function. Josef Eschgfäller __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] bootstrap vs. resampleing

2005-04-06 Thread array chip
Hi, I understand bootstrap can be used to estimate 95% confidence interval for some statistics, e.g. variance, median, etc. I have someone suggesting that by resampling certain proportion of the total samples (e.g. 80%) without replacement, we can also get the estimate of confidence intervals.

Re: [R] off-topic question: Latex and R in industries

2005-04-06 Thread A.J. Rossini
Ha -- that's a good one, Roger. Which demonstrates that most industrial people don't bother to read EULA's :-). (of course, it depends on which industry, and for some industries, which segment you are in ). On Apr 6, 2005 6:05 PM, roger koenker [EMAIL PROTECTED] wrote: my favorite answer to

Re: [R] off-topic question: Latex and R in industries

2005-04-06 Thread Marc Schwartz
On Wed, 2005-04-06 at 11:38 -0400, Wensui Liu wrote: Latex and R are really cool stuff. I am just wondering how they are used in industry. But based on my own experience, very rare. Why? How about the opinion of other listers? Thanks. As Tony has referenced, the answer will depend upon what

RE: [R] bootstrap vs. resampleing

2005-04-06 Thread Berton Gunter
I understand bootstrap can be used to estimate 95% confidence interval for some statistics, e.g. ^^ There's no such thing. You can estimate 95% CI's on population **parameters**, which is, I assume, what you mean. If you don't know what the difference is,

Re: [R] bootstrap vs. resampleing

2005-04-06 Thread Roger D. Peng
What you're describing sounds like subsampling, about which John Hartigan has written a few papers. -roger array chip wrote: Hi, I understand bootstrap can be used to estimate 95% confidence interval for some statistics, e.g. variance, median, etc. I have someone suggesting that by resampling

Re: [R] 2d plotting and colours

2005-04-06 Thread Earl F. Glynn
Mulholland, Tom [EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED] Since I was only concentrating on colour issues and not on your specific problem I was just showing the possibilities. Does this code help n - 5 par(mfrow = c(2,2)) palette(default) barplot(1:25,col = 1:25)

Re: [R] bootstrap vs. resampleing

2005-04-06 Thread Thomas Lumley
On Wed, 6 Apr 2005, array chip wrote: Hi, I understand bootstrap can be used to estimate 95% confidence interval for some statistics, e.g. variance, median, etc. I have someone suggesting that by resampling certain proportion of the total samples (e.g. 80%) without replacement, we can also get the

Re: [R] off-topic question: Latex and R in industries

2005-04-06 Thread roger bos
We have S+ at our company, but I choose to use R because I like it. There are two observations I have. One is that many people in IT don't seem to like open source software that much because either they don't trust it or they say there is no one who stands behind it. Second, equally important

[R] bootstrap vs. resampleing

2005-04-06 Thread Jens Oehlschlägel
Confidence intervals depend on the sample size - the bigger the sample the smaller the interval. Subsampling (resampling without replacement) gives smaller samples and underestimates confidence (overestimates confidence interval size) of parameters calculated on the original sample. Best Jens

[R] insignificant factors in regression model

2005-04-06 Thread Weijie Cai
Hi list, I am building a regression model with categorical predictor variable coded by treatment contrasts. The summary of the regression model shows that some levels are significant while others are not. The significant ones show that they are statistically significant from the basis factor

Re: [R] off-topic question: Latex and R in industries

2005-04-06 Thread Wensui Liu
Thank you all for the replies. I've used R and latex in graduate school and absolultely love them. After getting in the industry, everyone is using MS products or SPSS/SAS. But in term of quality, there is no comparison between MS word and Latex or between SAS/SPSS and R. On Apr 6, 2005 2:35

[R] using command line flags with TINN-R

2005-04-06 Thread roger bos
This is a TINN-R editor question rather than an R question, but can anyone tell me how to use command line flags with TINN-R. There is a space to fill in the path to Rgui, and I have C:\Program Files\R\rw2001pat\bin\Rgui.exe. If I try to add a command line flag after that, such as --no-save or

RE: [R] off-topic question: Latex and R in industries

2005-04-06 Thread Huntsinger, Reid
It appears that more than a few people posting on this list are from industry... I use R for research into signal/image analysis and tomography. Most people doing this use Matlab. I prefer R for several reasons. First, I'm a long-time S user. Second, I'm involved in the administration and in this

Re: [R] off-topic question: Latex and R in industries

2005-04-06 Thread Thomas Lumley
On Wed, 6 Apr 2005, roger bos wrote: Insightful will come in and give a company presentation. Who wants to volunteer to come into my company and demo R for my manager? I only learned about R a year ago when a friend of mine told me about it. The real question is, how to get more exposuRe? The

[R] newbie file size question

2005-04-06 Thread Christine Krisky
I've just recently started using R and have a very basic question. Sorry all, I've searched the postings and couldn't find the answer. I'm trying to read a very large ascii file 1,461,363 rows(data points) x 200 columns(time points) for a few hours now and R seems to be 'frozen'. I'm

RE: [R] off-topic question: Latex and R in industries

2005-04-06 Thread Ben Fairbank
I think there may be a bit of an us vs. them perception in business as it views academia (and R is a product of academia). I discussed the use of R with a businessman not long ago and he raised two objections to its use. First, if they give it away for free, how good can it be? After all, you

Re: [R] insignificant factors in regression model

2005-04-06 Thread Bob Wheeler
Weijie Cai wrote: Hi list, I am building a regression model with categorical predictor variable coded by treatment contrasts. The summary of the regression model shows that some levels are significant while others are not. The significant ones show that they are statistically significant from

[R] nnet classification using unbalanced classes

2005-04-06 Thread Tarca Adi Laurentiu
Hi everybody, I want to obtain a classification model using the nnet function for a simple two class problem. My problem is that number of samples in the first class (n1) is about twice higher than the one in class two (n2). I would like to use the weights argument in the nnet function to

Re: [R] off-topic question: Latex and R in industries

2005-04-06 Thread Jan T. Kim
On Wed, Apr 06, 2005 at 11:56:59AM -0700, Thomas Lumley wrote: On Wed, 6 Apr 2005, roger bos wrote: Insightful will come in and give a company presentation. Who wants to volunteer to come into my company and demo R for my manager? I only learned about R a year ago when a friend of mine told

[R] Polynomiographic function in R :-)

2005-04-06 Thread François Pinard
Hi, people. Nothing too serious in this message. Nevertheless, all criticism or advice is welcome :-). Yesterday, I went to a conference by Bahman Kalantari (Rutgers University) about Polynomiography (the Fine Art and Science of Visualizing Polynomials). Since I'm starting my R learning, I

[R] looking for a plot function

2005-04-06 Thread bogdan romocea
Dear useRs, I have a data frame and I want to plot all rows. Each row is represented as a line that links the values in each column. The plot looks like this: dfr - data.frame(A=sample(1:50,10),B=sample(1:50,10), C=sample(1:50,10),D=sample(1:50,10)) xa - 10*1:4 plot(c(10,40),c(0,50)) for

Re: [R] looking for a plot function

2005-04-06 Thread Deepayan Sarkar
On Wednesday 06 April 2005 15:18, bogdan romocea wrote: Dear useRs, I have a data frame and I want to plot all rows. Each row is represented as a line that links the values in each column. The plot looks like this: dfr - data.frame(A=sample(1:50,10),B=sample(1:50,10),

RE: [R] looking for a plot function

2005-04-06 Thread Doran, Harold
If your data were in the long format you could use interaction.plot. But, I think trellis plots in lattice are much better than this approach. -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of bogdan romocea Sent: Wednesday, April 06, 2005 4:19 PM To:

RE: [R] looking for a plot function

2005-04-06 Thread apjaworski
I am not sure about the scaling, but doing simply matplot(xa, t(dfr), type=b) does most of what you want. Andy __ Andy Jaworski 518-1-01 Process Laboratory 3M Corporate Research Laboratory - E-mail: [EMAIL PROTECTED] Tel: (651) 733-6092 Fax: (651)

Re: [R] Is a .R script file name available inside the script?

2005-04-06 Thread Darren Weber
That is useful, when calling the script like this: file - Rscript.R source(file) However, it does not work if we do this from the shell prompt: $ R --vanilla Rscript.R because the eval.parent statement attempts to access a base workspace that does not contain the file object/variable, as

Re: [R] off-topic question: Latex and R in industries

2005-04-06 Thread Patrick Burns
Thomas Lumley wrote: The other real question is Why?. I can see the motivation of people who want to use R and need to convince their management that it is safe, but inflicting R on people who haven't heard of it and are perfectly happy that way seems unnecessary. What would be the benefit?

RE: [R] bootstrap vs. resampleing

2005-04-06 Thread Huntsinger, Reid
I may be misunderstanding the question, but I believe you want a pointwise confidence band for the conditional odds function. The issue here is less bootstrap versus some other resampling plan, and more how to do it at all. For example, if no matter what training data you feed in, you always get

[R] read.table with header and text data

2005-04-06 Thread Laura Holt
Hi R! I am reading in a text file which has one column of alpha data and 5 columns of numeric data. There is a header row. I would like the alpha data column to just be character rather than factor. Is there a way to do this, please? I'm thinking that it might be I() but can't figure out

Re: [R] read.table with header and text data

2005-04-06 Thread Rich FitzJohn
See ?read.table, especially the argument as.is. Cheers, Rich On Apr 7, 2005 9:55 AM, Laura Holt [EMAIL PROTECTED] wrote: Hi R! I am reading in a text file which has one column of alpha data and 5 columns of numeric data. There is a header row. I would like the alpha data column to

Re: [R] par(mfcol=2, mfrow=3) equivalent for trellis

2005-04-06 Thread Paul Murrell
Hi Deepayan Sarkar wrote: On Wednesday 06 April 2005 09:06, Dieter Menne wrote: Dear friends of lattice, I know how to position trellis plots with print(...,split,more=T) or (...position). Sometimes I wish I had something like the old par(mfcol=2, mfrow=3) mechanism, where the next free viewport

[R] Importing SAS transport data

2005-04-06 Thread Dennis Fisher
I have encountered a problem reading SAS transport files on both a Mac (OS X) and Linux (RedHat 9), both using R2.0.1. After loading foreign, the command: read.xport(V1622101_050304.xpt) yields: Error in lookup.xport(file) : File not in SAS transfer format In Linux, I can cat the file. The

Re: [R] 2d plotting and colours

2005-04-06 Thread Paul Murrell
Hi Earl F. Glynn wrote: Mulholland, Tom [EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED] Since I was only concentrating on colour issues and not on your specific problem I was just showing the possibilities. Does this code help n - 5 par(mfrow = c(2,2)) palette(default) barplot(1:25,col =

Re: [R] R can not show plots (in Mac OS X terminal)

2005-04-06 Thread George W. Gilchrist
This is incorrect. x11(display=0:0) opens an x11 graphics device from the terminal assuming (1) that you have installed X11 from Apple's website and (2) that x11 is running. Cheers, George On 4/6/05 5:47 AM, David Ruau [EMAIL PROTECTED] wrote: Hi, You should use X11. It doesn't work in

Re: [R] Is a .R script file name available inside the script?

2005-04-06 Thread Gabor Grothendieck
It works for me. Suppose in.txt is a two line file with these two lines: file - Rscript.R source(file) and Rscript.R is a two line file with these two lines: script.description - function() eval.parent(quote(file), n = 3) print(basename(script.description())) Then here is the output on

[R] newline in lattice axis label

2005-04-06 Thread Sebastian Luque
Hi, I have a 3 panel xyplot with different variables in the y axis. I'm trying to insert a newline after Width (cm), in the ylab argument as in the example below. My goal is to have the y axis label broken into two lines, split after the string just mentioned. plotfun - function() { fakedf -

Re: [R] newline in lattice axis label

2005-04-06 Thread Sebastian Luque
Sebastian Luque [EMAIL PROTECTED] wrote: I have a 3 panel xyplot with different variables in the y axis. I'm trying to insert a newline after Width (cm), Sorry, that should be VarC (cm), in the hypothetical example below! -- Sebastian P. Luque __

Re: [R] How to do aggregate operations with non-scalar functions

2005-04-06 Thread Itay Furman
On Tue, 5 Apr 2005, Gabor Grothendieck wrote: On Apr 5, 2005 6:59 PM, Itay Furman [EMAIL PROTECTED] wrote: Hi, I have a data set, the structure of which is something like this: a - rep(c(a, b), c(6,6)) x - rep(c(x, y, z), c(4,4,4)) df - data.frame(a=a, x=x, r=rnorm(12)) The true data set has 1

Re: [R] How to do aggregate operations with non-scalar functions

2005-04-06 Thread Itay Furman
On Wed, 6 Apr 2005, Rich FitzJohn wrote: [snip] ## This does the hard work of calculating the statistics over your ## combinations, and over the values in `p' y - lapply(p, function(y) tapply(df$r, list(a=a, x=x), quantile, probs=y)) Rich, thank you for your reply. Gabor G has

Re: [R] newline in lattice axis label

2005-04-06 Thread Renaud Lancelot
Sebastian Luque a écrit : Hi, I have a 3 panel xyplot with different variables in the y axis. I'm trying to insert a newline after Width (cm), in the ylab argument as in the example below. My goal is to have the y axis label broken into two lines, split after the string just mentioned. plotfun -

Re: [R] How to do aggregate operations with non-scalar functions

2005-04-06 Thread Itay Furman
On Wed, 6 Apr 2005 [EMAIL PROTECTED] wrote: Here is a method that I use in this situation. I work with the indices of the rows so that copies are not made and it is fast. Result - lapply(split(seq(nrow(df)), df$a), function(.a){ # partition on the first variable lapply(split(.a, df$z[.a]),

Re: [R] using command line flags with TINN-R

2005-04-06 Thread Philippe Grosjean
roger bos wrote: This is a TINN-R editor question rather than an R question, but can anyone tell me how to use command line flags with TINN-R. There is a space to fill in the path to Rgui, and I have C:\Program Files\R\rw2001pat\bin\Rgui.exe. If I try to add a command line flag after that, such