Re: [R] Kmeans centers

Thank you very much Gavin, The set.seed is the correct function i need. Now the kmeans is permanent and doesn't change results every time i run. 2007/3/30, Gavin Simpson <[EMAIL PROTECTED]>: > > On Fri, 2007-03-30 at 09:07 +0200, Sergio Della Franca wrote: > > My simple problem is that when

Re: [R] R equivalent of S+SeqTrial?

Thanks again; It's ok I'm using the ldBands function now, the ld98.exe is just a small file to download into the path. I can calculate the alpha spending function and determine (given the Z statistics of my observations) if at the interim analysis the null hypothesis is to be rejected or not. But

Re: [R] Hmisc summary.formula.reverse export problem

Thanks Brian, I was trying to get that object printed out with columns separated by semicolon. The reason for this is that I'm trying to get it to MS Excel. If this is not possible, I probably should turn into using Latex. -Lauri 2007/3/30, Prof Brian Ripley <[EMAIL PROTECTED]>: > > Well, wri

Re: [R] Wikibooks

Bert Gunter wrote: > Question: > > Many (perhaps most?) questions on the list are easily answerable simply by > checking existing R Docs (Help file/man pages, Intro to R, etc.). Why would > a Wiki be more effective in deflecting such questions from the mailing list > than them? Why would too hel

Re: [R] Wikibooks

..<°}))>< ) ) ) ) ) ( ( ( ( (Prof. Philippe Grosjean ) ) ) ) ) ( ( ( ( (Numerical Ecology of Aquatic Systems ) ) ) ) ) Mons-Hainaut University, Belgium ( ( ( ( ( .. Diet

[R] substitute NA values

Dear R-Helpers, I have the following data set(y): Test_Result #_Test t 10 f 14 f 25 f NA f 40 t45 t44 47 tNA I want to

Re: [R] Wikibooks

On 3/30/2007 5:27 AM, Philippe Grosjean wrote: > > Bert Gunter wrote: >> Question: >> >> Many (perhaps most?) questions on the list are easily answerable simply by >> checking existing R Docs (Help file/man pages, Intro to R, etc.). Why would >> a Wiki be more effective in deflecting such question

Re: [R] Wikibooks

Philippe Grosjean wrote: > > As other have pointed out, the main reason for the lack of success > of the R Wiki is that the mailing lists, particularly R-Help, are > sooo successful. However, I continue to consider that the mailing > list is suboptimal in two cases: (1) when text is not enough

Re: [R] Regarding Vista

--- Dieter Menne <[EMAIL PROTECTED]> wrote: > ccilindia.co.in> writes: > > > > > I am facing the same problem in my case. R 2.4.1 > have installed > > successfully, but when i try to install the > packages from a local zip > > file. It gives the following error message. > > +++

Re: [R] Regarding Vista

John Kane wrote: > > As a somewhat desperate workaround try installing R on > a USB and see if you can run if from there. I have > 2.4.1 on a USB and it seems to work fine albeit a bit > more slowly than from the hard drive. > I love this quote, in http://zoonek2.free.fr/UNIX/48_R/02.html, fro

[R] Model comparison

Dear all, I would like to know if I can compare by a significance test 2 models with different kind of parameters. Perhaps I am wrong but I think that we can only compare 2 models if one is a sub model of the other. Med venlig hilsen / Regards João Fadista Ph.d. studerende / Ph.d. student

Re: [R] Wikibooks

On 3/30/2007 7:34 AM, Alberto Monteiro wrote: > Duncan Murdoch wrote: >> >> But the wiki doesn't offer a way to ask questions. I'd be just as >> happy to answer questions there as here, but there are none there to >> answer >> (and the advice there is to ask questions here). >> >> I don't kno

Re: [R] Tail area of sum of Chi-square variables

Hi, thanks everyone! pchisqsum() in the "survey" package does exactly what I was looking for! Best wishes, Klaus Original-Nachricht Datum: Thu, 29 Mar 2007 07:45:15 -0700 (PDT) Von: Thomas Lumley <[EMAIL PROTECTED]> An: S Ellison <[EMAIL PROTECTED]> CC: [EMAIL PROTECTED], [

[R] Inquiry

Good morning, I have a question about R, I would like to know how it is possible not to chop the result of an operation. How many decimals it is possible to obtain? Thank you in advance, I. López - Inmaculada López García Dpto. Estadística y Matemátic

Re: [R] Wikibooks

Duncan Murdoch wrote: > > But the wiki doesn't offer a way to ask questions. I'd be just as > happy to answer questions there as here, but there are none there to > answer > (and the advice there is to ask questions here). > > I don't know how to organize a wiki to make it easy to ask and >

Re: [R] Wikibooks

> > > > I once tried: > > > > http://wiki.r-project.org/rwiki/doku.php?id=guides:lmer-tests > > > > but I don't think I will do this again on the existing Wiki. I am a frequent > > Wikipedia-Writer, so I know how it works, but this was discouraging. > > > > 1) The structure of the Wiki was and is s

Re: [R] substitute NA values

I assume you are referring to na.roughfix in randomForest. I don't think it works for logical vectors or for factors outside of data frames: > library(randomForest) > DF <- data.frame(a = c(T, F, T, NA, T), b = c(1:3, NA, 5)) > na.roughfix(DF) Error in na.roughfix.data.frame(DF) : na.roughfix onl

Re: [R] Wikibooks

Romain Francois wrote: > >> Say I don't know (and I can't understand the help) how to >> use the rnorm function. If I do RSiteSearch("rnorm"), I >> will get too much useless information. OTOH, an ideal wikipedia >> would have a page http://www.r-wiki.org/rnorm, where I could >> find examples, learn

Re: [R] Wikibooks

On 3/30/2007 9:16 AM, Alberto Monteiro wrote: > Romain Francois wrote: >> >>> Say I don't know (and I can't understand the help) how to >>> use the rnorm function. If I do RSiteSearch("rnorm"), I >>> will get too much useless information. OTOH, an ideal wikipedia >>> would have a page http://www.r-

[R] math-operations

Hallo R-experts, for a function I need to work with the commands "div" and "mod" known from Pascal and Ruby. The only help I know is "ceiling()" and "floor()" in R. Do "div" and "mod" exist? When yes please send me a little example. In Pascal-Syntax I want: 513 div 100 = 5 513 mod 100 = 13 How

Re: [R] math-operations

513 %/% 100 513 %% 100 check ?"%/%" for more info. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/(0)16/336899 Fax: +32/(0)16/337015 Web: http://med.kuleuven

Re: [R] CLUSTER Package

It seems nobody else was willing to help here (when the original poster did not at all follow the posting guide). In the mean time, someone else has asked me about part of this, so let me answer in public : > "MM" == Martin Maechler <[EMAIL PROTECTED]> > on Mon, 12 Mar 2007 17:23:30 +

Re: [R] math-operations

On Friday 30 March 2007 14:34, Schmitt, Corinna wrote: > Hallo R-experts, > > for a function I need to work with the commands "div" and "mod" known > from Pascal and Ruby. The only help I know is "ceiling()" and "floor()" > in R. Do "div" and "mod" exist? When yes please send me a little > example.

[R] Inquiry

Petr Pikal [EMAIL PROTECTED] - Postoupil do Petr PIKAL/CTCAP dne 30.03.2007 15:47 - Petr PIKAL/CTCAP napsal dne 30.03.2007 15:01:58: > Hi > > Do you mean rounding. If yes you can consult > ?round, ?floor, ?ceiling. > > If you want just to print different amount of digits you can look a

[R] R: math-operations

I guess you need "%%" and ""%/%" try > 513 %/% 100 [1] 5 > 513 %% 100 [1] 13 ?"%%" Stefano -Messaggio originale- Da: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] conto di Schmitt, Corinna Inviato: venerdì 30 marzo 2007 15.34 A: r-help@stat.math.ethz.ch Oggetto: [R] math-operations

Re: [R] math-operations

On 3/30/2007 9:34 AM, Schmitt, Corinna wrote: > Hallo R-experts, > > for a function I need to work with the commands "div" and "mod" known > from Pascal and Ruby. The only help I know is "ceiling()" and "floor()" > in R. Do "div" and "mod" exist? When yes please send me a little > example. > > In

Re: [R] math-operations

On Friday 30 March 2007 14:46, Alberto Murta wrote: > On Friday 30 March 2007 14:34, Schmitt, Corinna wrote: > > Hallo R-experts, > > > > for a function I need to work with the commands "div" and "mod" known > > from Pascal and Ruby. The only help I know is "ceiling()" and "floor()" > > in R. Do "d

[R] faster computation of cumulative multinomial distribution

Dear list members, I have a series of /unequal/ probabilities [p1,p2,...,pk], describing mutually exclusive events, and a "remainder" class with a probability p0=1-p1-p2--pk, and need to calculate, for a given number of trials t>=k, the combined probability that each of the classes 1...k conta

Re: [R] math-operations

Thanks, it works. Corinna -Ursprüngliche Nachricht- Von: Guazzetti Stefano [mailto:[EMAIL PROTECTED] Gesendet: Freitag, 30. März 2007 15:41 An: Schmitt, Corinna; r-help@stat.math.ethz.ch Betreff: R: [R] math-operations I guess you need "%%" and ""%/%" try > 513 %/% 100 [1] 5 > 513

Re: [R] math-operations and the R-Wiki

Dimitris Rizopoulos wrote: > > 513 %/% 100 > > 513 %% 100 > Now this is a great opportunity to improve the R-Wiki. What about a "Pascal" page in the R-Wiki, where a list of Pascal-to-R translations would be available? Alberto Monteiro __ R-help@stat.m

Re: [R] substitute NA values

This is that i obtained. There isn't a method to replace the NA values only for character variable? 2007/3/30, Gabor Grothendieck <[EMAIL PROTECTED]>: > > I assume you are referring to na.roughfix in randomForest. I don't think > it > works for logical vectors or for factors outside of data

Re: [R] Wikibooks

Duncan Murdoch wrote: > >> This works when there's a decent documentation for the function. >> The functions in the tcltk package, for example, are horribly >> undocumented, and asking for help only loops to a general >> help about all (and none) of the functions. > > I don't remember if you've s

[R] ANOVA and confidence intervals plot

Dear *, I would like to obtain for each factor of my anova model the "response variable vs factor" plot with means and 95% Tukey HSD intervals. I would appreciate any information on how to do that. Cheers

Re: [R] Inquiry

I am not sure I understand your question but have a look at ?round and the signif command on that page. --- Inmaculada López García <[EMAIL PROTECTED]> wrote: > > > Good morning, > > I have a question about R, I would like to know how > it is possible not to > chop the result of an operation

Re: [R] substitute NA values

Not as part of na.roughfix. You could convert your character strings to factors and back again: library(randomForest) DF <- data.frame(a = c(T, F, T, NA, T), b = c(1:3, NA, 5), c = c("b", "b", NA, "d", "e"), d = factor(c("a", "a", NA, "d", "e")), stringsAsFactors = FALSE) DF$a <- factor(D

Re: [R] Wikibooks

Duncan Murdoch wrote: > >> This works when there's a decent documentation for the function. >> The functions in the tcltk package, for example, are horribly >> undocumented, and asking for help only loops to a general >> help about all (and none) of the functions. > > I don't remember if you've s

Re: [R] substitute NA values

On Fri, 2007-03-30 at 16:25 +0200, Sergio Della Franca wrote: > This is that i obtained. > > There isn't a method to replace the NA values only for character variable? This is R, there is always a way (paraphrasing an R-Helper the name of whom I forget just now). If you mean a canned function, no

[R] Testing random effect in logistic mixed model

Hi When using linear mixed model, I could test for the effect of the random part of the model using a likelihood ratio test comparing two model with and without the random part. model.lmer=lmer(y~x+(1|r)) model.lm = lm(y~x) anova(model.lmer,model.lm) However, this does not work with a mixed mode

[R] Using split() several times in a row?

Hi, fellow R users. I have a question about sapply and split combination. I have a big dataframe (4 observations, 21 variables). First variable (factor) is "date" and it is in format "8.29.97", that is, I have monthly data. Second variable (also factor) has levels 1 to 6 (fractiles 1 to 5 and

Re: [R] Inquiry

My interpretation of the question (which of course may be wrong) has an answer that is the opposite of 'round': When a result is printed, it is rounded to a certain number of digits (controlled by the 'digits' argument of 'print'). Just because it is printed like that, doesn't mean the actual val

[R] Minimum valid number of observations for rpart

Hi, I wonder if anyone knows a study dealing with the minimum valid number of observations when using CART?. On top of that, when using RandomForest, is it possible to obtained a interpretable tree model as the graphical output of the analysis, just like in "rpart"? Thanks a lot in advance Javier

[R] Makeconf on Windows - where?

Read in the Windows release notes for 2.4.0 that: There is a new file etc/Makeconf that provides an approximation to the Unix version and may help with src/Makefile's. Yet when I look at the CRAN binary I have for 2.4.1patched, there are five files in the $R_HOME/etc directory, but no Makec

Re: [R] Wikibooks

Alberto Monteiro wrote: > Philippe Grosjean wrote: > >> As other have pointed out, the main reason for the lack of success >> of the R Wiki is that the mailing lists, particularly R-Help, are >> sooo successful. However, I continue to consider that the mailing >> list is suboptimal in two cas

Re: [R] Wikibooks

On 3/30/07, Alberto Monteiro <[EMAIL PROTECTED]> wrote: > Philippe Grosjean wrote: > > > > As other have pointed out, the main reason for the lack of success > > of the R Wiki is that the mailing lists, particularly R-Help, are > > sooo successful. However, I continue to consider that the mailing >

Re: [R] Wikibooks

> If entering a new page is really the way to ask a question, then you > should write this on the front page, and as a possible way to contribute > on the getting-started page. It would also be a good idea to tell > people like me how to find those questions. ("Recent Edits" seems a > little too b

Re: [R] Makeconf on Windows - where?

On 3/30/2007 11:53 AM, Paul Roebuck wrote: > Read in the Windows release notes for 2.4.0 that: > > There is a new file etc/Makeconf that provides an > approximation to the Unix version and may help with > src/Makefile's. > > Yet when I look at the CRAN binary I have for 2.4.1patched, > ther

Re: [R] Using functions in LAPACK in a C program

Many thanks to Brian for his very useful information and help. A quick look at mgcv package already gives me some directions to try. I remember that I tried to use Rlapack.dll instead of R.dll without success, perhaps because I don't know how to build an import library for Rlapack.dll and link a

Re: [R] Minimum valid number of observations for rpart

On Fri, 30 Mar 2007, [EMAIL PROTECTED] wrote: Hi, I wonder if anyone knows a study dealing with the minimum valid number of observations when using CART?. I have no idea what you mean by 'valid' here. Could you answer the question for logistic regression to indicate to us what form of answer

Re: [R] Using functions in LAPACK in a C program

I think all these points are answered in the README.packages for R 2.5.0 alpha (and probably for 2.4.1 too). You can use pexports to create Rlapack.exp, for example. On Fri, 30 Mar 2007, [EMAIL PROTECTED] wrote: Many thanks to Brian for his very useful information and help. A quick look at mg

Re: [R] Model comparison

On Fri, 30 Mar 2007, Jo?o Fadista wrote: Dear all, I would like to know if I can compare by a significance test 2 models with different kind of parameters. Perhaps I am wrong but I think that we can only compare 2 models if one is a sub model of the other. The literature you seek is on '

Re: [R] Minimum valid number of observations for rpart

Regarding interpretable output, I assume you have looked at the mds plot? regards, Farrar [EMAIL PROTECTED] wrote: Hi, I wonder if anyone knows a study dealing with the minimum valid number of observations when using CART?. On top of that, when using RandomForest, is it possible

[R] Trouble installing package 'sp' in R 2.4.1

Hi all, I just installed R 2.4.1 and tried to install the package 'sp' working down from the 'Packages' menu in the R GUI. The package would not install and I got the following message: "trying URL 'http://cran.us.r-project.org/bin/windows/contrib/2.4/sp_0.9-10.zip' Error in download.file(url

Re: [R] Trouble installing package 'sp' in R 2.4.1

On 3/30/2007 2:39 PM, Gaurav Ghosh wrote: > Hi all, > > I just installed R 2.4.1 and tried to install the package 'sp' working down > from the 'Packages' menu in the R GUI. The package would not install and I > got the following message: > > "trying URL > 'http://cran.us.r-project.org/bin/win

Re: [R] Wikibooks

On 3/30/07, Dieter Menne <[EMAIL PROTECTED]> wrote: > Ben Bolker zoo.ufl.edu> writes: > > > Well, we do have an R wiki -- http://wiki.r-project.org/rwiki/doku.php -- > > although it is not as active as I'd like. (We got stuck halfway through > > porting Paul Johnson's "R Tips" to it ...) Plea

[R] RWiki, tcltk and plot

I think I - almost - got the knack for GUI programming using the tcltk library. Maybe I will update the RWiki with this: # # # library(tcltk) # # Create some matrix - nothing about tcltk here # matrix <- cbind(rnorm(100), rpois(100, lambda=10),

Re: [R] Wikibooks

Deepayan Sarkar wrote: > > I was just looking at this page, and it makes me curious: what gives > anyone the right to take someone else's mailing list post and include > that in a Wiki? > Thinks there were posted to public mailing lists are freely copied and distributed. It's a scary thought; I m

Re: [R] Trouble installing package 'sp' in R 2.4.1

On Fri, 30 Mar 2007, Duncan Murdoch wrote: > On 3/30/2007 2:39 PM, Gaurav Ghosh wrote: > > Hi all, > > > > I just installed R 2.4.1 and tried to install the package 'sp' working down > > from the 'Packages' menu in the R GUI. The package would not install and I > > got the following message: >

Re: [R] Wikibooks

On 3/30/07, Alberto Monteiro <[EMAIL PROTECTED]> wrote: > Deepayan Sarkar wrote: > > > > I was just looking at this page, and it makes me curious: what gives > > anyone the right to take someone else's mailing list post and include > > that in a Wiki? > > > Thinks there were posted to public mailin

Re: [R] Wikibooks

> Under US law (the only one I'm familiar with), the author of a mailing list > post or any other written work _automatically holds copyright_ to that > post (although not to the ideas contained therein, but to that particular > description of the ideas). That's true in almost any country - see th

Re: [R] RWiki, tcltk and plot

On Fri, Mar 30, 2007 at 05:28:14PM -0200, Alberto Monteiro wrote: > My question: is there any way to integrate the plot part into > a tcltk window? Are you aware of the tkrplot package on CRAN ? Dirk -- Hell, there are no rules here - we're trying to accomplish something.

Re: [R] Wikibooks

On 3/30/07, Sarah Goslee <[EMAIL PROTECTED]> wrote: > On 3/30/07, Alberto Monteiro <[EMAIL PROTECTED]> wrote: > > Deepayan Sarkar wrote: > > > > > > I was just looking at this page, and it makes me curious: what gives > > > anyone the right to take someone else's mailing list post and include > > >

Re: [R] RWiki, tcltk and plot

Look at the tkrplot package. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] (801) 408-8111 > -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of > Alberto Monteiro > Sent: Friday, March 30, 2007 1

[R] zip.file.extract on Windows

Hello R-help! I just wanted to share a tip might help others of us stuck on R. The documentation for zip.file.extract is rather scant, and a few examples would. When using zip.file.extract, the first argument should be the full (or relative) path and the second the name of the zip file, it s

Re: [R] RWiki, tcltk and plot

Dirk Eddelbuettel wrote: > >> My question: is there any way to integrate the plot part into >> a tcltk window? > > Are you aware of the tkrplot package on CRAN ? > No. Alberto Monteiro __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mail

Re: [R] Wikibooks

On a related note, one might be interested in checking out citizendium which is spin off wikipedia but 1) has more stringent identity verification and 2) uses a two-tier system of editors and authors. See http://www.citizendium.org/cfa.html. Deepayan Sarkar wrote: > On 3/30/07, Sarah Goslee <

Re: [R] Wikibooks

Deepayan Sarkar wrote: > On 3/30/07, Sarah Goslee <[EMAIL PROTECTED]> wrote: > >> On 3/30/07, Alberto Monteiro <[EMAIL PROTECTED]> wrote: >> >>> Deepayan Sarkar wrote: >>> I was just looking at this page, and it makes me curious: what gives anyone the right to take someone

Re: [R] Copyright of messages (was Wikibooks)

What you agree to _is_ covered in the R posting guide and its references, as well as what is expected of posters in respecting the rights of others. Perhaps time yet again to remind people: > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html On Fri, 30 Mar 2007, Pete

Re: [R] Wikibooks

On Fri, 30 Mar 2007, Peter Dalgaard wrote: > Deepayan Sarkar wrote: > > On 3/30/07, Sarah Goslee <[EMAIL PROTECTED]> wrote: > > > >> On 3/30/07, Alberto Monteiro <[EMAIL PROTECTED]> wrote: > >> > >>> Deepayan Sarkar wrote: > >>> > I was just looking at this page, and it makes me curious: what

Re: [R] faster computation of cumulative multinomial distribution

Theo Borm wrote: > > I have a series of /unequal/ probabilities [p1,p2,...,pk], describing > mutually exclusive events, and a "remainder" class with a probability > p0=1-p1-p2--pk, and need to calculate, for a given number of trials > t>=k, the combined probability that each of the classes 1..

Re: [R] substitute NA values

On Mar 30, 2007, at 10:56 AM, Gavin Simpson wrote: > This is R, there is always a way (paraphrasing an R-Helper the name of > whom I forget just now). Can't resist, it's one of my favorite fortunes ;) That would be Simon 'Yoda' Blomberg: library(fortunes) fortune(109) Haris Skiadas Department

Re: [R] faster computation of cumulative multinomial distribution

Alberto Monteiro wrote: > Theo Borm wrote: > >> I have a series of /unequal/ probabilities [p1,p2,...,pk], describing >> mutually exclusive events, and a "remainder" class with a probability >> p0=1-p1-p2--pk, and need to calculate, for a given number of trials >> t>=k, the combined probabil

[R] X11 and linux and plotting and the vertical axis

I use R on linux and I go through exceed from a windows machine. Depending on the amount of plots I do on a screen, sometimes the numbers on the vertical axis of the plots don't show up. I tried infinite of combinations of mar and cex but nothing helps. Yet, if instead of shooting the graphic upto

Re: [R] Wikibooks

On 3/30/2007 5:05 PM, Peter Dalgaard wrote: > Deepayan Sarkar wrote: >> On 3/30/07, Sarah Goslee <[EMAIL PROTECTED]> wrote: >> >>> On 3/30/07, Alberto Monteiro <[EMAIL PROTECTED]> wrote: >>> Deepayan Sarkar wrote: > I was just looking at this page, and it makes me curio

Re: [R] Using split() several times in a row?

Hi Sergey, I believe the code below should get you close to want you want. For dates, I usually store them as "POSIXct" classes in data frames, but according to Gabor Grothendieck and Thomas Petzoldt's R Help Desk article , I should probably b

Re: [R] Trouble installing package 'sp' in R 2.4.1

On 3/30/2007 3:37 PM, Roger Bivand wrote: > On Fri, 30 Mar 2007, Duncan Murdoch wrote: > >> On 3/30/2007 2:39 PM, Gaurav Ghosh wrote: >>> Hi all, >>> >>> I just installed R 2.4.1 and tried to install the package 'sp' working down >>> from the 'Packages' menu in the R GUI. The package would not i

Re: [R] faster computation of cumulative multinomial distribution

Theo Borm wrote: > > Think of a sort of "power roulette", played with 58 balls > simultaneously, with a wheel containing 36 red/black slots of unequal > size, and 1 green slot. I need to calculate the probability that each of > the 36 red/black slots contains at least one ball. > Ah, now we come to

[R] Linear model and time series

Dear all, I have three timeseries Uts, Vts, Wts. The relation between the time series can be expressed as Uts = x Vts + y Wts + residuals How would I feed this to lm() to evaluate the unknowns x and y? Thanks, andre __ R-help@stat.math.ethz.ch mailin