Re: [R] ED50 from logistic model with interactions

2007-05-02 Thread Berwin A Turlach
On Wed, 02 May 2007 11:37:22 +1000 Kate Stark [EMAIL PROTECTED] wrote: [...] My model is: fit - glm(Mature ~ Season * Size - 1, family = binomial, data=dat) where Mature is a binary response, 0 for immature, 1 for mature. There are 3 Seasons. I would use: fit - glm(Mature

Re: [R] Simulation using parts of density function

2007-05-02 Thread Thür Brigitte
Thanks for your code! It is not exactly what I really want - but it is my fault, because my description was wrong... It is not sim but rhater exp(rgamma(...)) that should not exceed 500. So I tried to modify your code but it doesn't really work. sim.test returns just 1 value and not 999.

[R] Log-likelihood function

2007-05-02 Thread Doxastic
I've computed a loglinear model on a categorical dataset. I would like to test whether an interaction can be dropped by comparing the log-likelihoods from two models(the model with the interaction vs. the model without). Since R does not immediately print the log-likelihood when I use the glm

[R] Get the difference of values to their own median value

2007-05-02 Thread Felix Wave
Hello, I've got a matrix (mail end) with the colnames x, y, z. In this matrix are different measurements. x and y are risign coordinates. With the following line I got the median value of z for all x AND y witch are the same (not every measurment in my list hast the same number of x and y

Re: [R] Simulation using parts of density function

2007-05-02 Thread Prof Brian Ripley
Please do not send everything twice: you are using R-help in both the To: and Cc: fields. I disagree with Ted: it _is_ much easier to create a generator for this purpose. Consider rtgamma - function(n, ..., tr = log(500)) { p - pgamma(tr, ...) qgamma(p*runif(n), ...) } as

[R] missing package

2007-05-02 Thread raymond chiruka
l'm trying to find the survdiff package/function but it seems i cant get it how do l instal it if its not there thanks - [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing

Re: [R] Log-likelihood function

2007-05-02 Thread Prof Brian Ripley
I think you need to learn about deviances, which R does print. Log-likelihoods are only defined up to additive constants. In this case the conventional constant differs if you view this as a Poisson or as a product-multinomial log-linear model, and R gives you the log-likelihood for a Poisson

Re: [R] ? R 2.5.0 alpha bug

2007-05-02 Thread Stefan Grosse
The version 2.5.0 has left Alpha status long time ago and its final version has been released so please try the new version. Inman, Brant A. M.D. wrote: This email is intended to highlight 2 problems that I encountered running R 2.5.0 alpha on a Windows XP machine. #1 - Open script error

Re: [R] Simulation using parts of density function

2007-05-02 Thread Ted Harding
On 02-May-07 07:45:48, Prof Brian Ripley wrote: Please do not send everything twice: you are using R-help in both the To: and Cc: fields. I disagree with Ted: it _is_ much easier to create a generator for this purpose. Consider rtgamma - function(n, ..., tr = log(500)) { p

[R] hello

2007-05-02 Thread elyakhlifi mustapha
hello, I have a problem with a R program I don't understand my errors my program looks like this for(i in LE) { +for(j in LEC[[i]]) { +for(k in LR) { + donRep[[k]] - subset(don2, Id_Essai == 1006961 Id_Cara == j Id_Rep == k Id_Geno != 65125,

Re: [R] Log-likelihood function

2007-05-02 Thread Doxastic
You're right. I do need to learn more. I never learned null/residual deviance. I know the deviance is equivalent to an anova decompostion. But I've never dealt with it seperated like this. I understand deviance as the difference between two model's log-likelihood difference between them

Re: [R] missing package

2007-05-02 Thread John Kane
--- raymond chiruka [EMAIL PROTECTED] wrote: l'm trying to find the survdiff package/function but it seems i cant get it how do l instal it if its not there thanks I think you want the survival package. __ R-help@stat.math.ethz.ch mailing list

[R] missing values

2007-05-02 Thread elyakhlifi mustapha
hello, I need your help for this example for(k in LR) { + donGeno[[k]] - as.numeric(levels(factor(subset(don2, Id_Essai == 1006961 Id_Cara == LC[1] Id_Rep == k, select = Id_Geno)[,1]))) + print(donGeno[[k]])} [1] 65125 65126 65127 65128 65129 65130 65131 65132 65133 65134 65135 65136 65137

Re: [R] missing package

2007-05-02 Thread Petr Klasterecky
help.search(survdiff) tells you it is in package survival. Most probably this package is not missing, nevertheless if it is, you can install it in a standard way using install.packages() or from the menu (Windows). Petr raymond chiruka napsal(a): l'm trying to find the survdiff

[R] length of the object is not a multiple of... [was: Re: hello]

2007-05-02 Thread Petr Klasterecky
elyakhlifi mustapha napsal(a): hello, I have a problem with a R program I don't understand my errors my program looks like this for(i in LE) { +for(j in LEC[[i]]) { +for(k in LR) { + donRep[[k]] - subset(don2, Id_Essai == 1006961 Id_Cara ==

Re: [R] to draw a smooth arc

2007-05-02 Thread Jim Lemon
(Ted Harding) wrote: This thread prompts me to ask about something I've been pondering for a while, as to whether there's an implementation somewhere ticked away in the R resources. So far, people have been responding to the original query in terms of increasing the numbers of points, and

Re: [R] Polar graph of time and tide

2007-05-02 Thread Jim Lemon
Alan E. Davis wrote: I have been trying to visualize times of lowest tides, month by month. I have tide predictions with times either in unix time or a text format, and heights in feet or meters. I had been able to derive the clock times of each prediction. I would now like to graph this

Re: [R] Log-likelihood function

2007-05-02 Thread Ross Darnell
Alternatively generate the log-likelihood using the sum(dpois(y, fitted(model), log = TRUE)) Regards Ross Darnell Doxastic wrote: You're right. I do need to learn more. I never learned null/residual deviance. I know the deviance is equivalent to an anova decompostion. But I've

Re: [R] Concepts question: environment, frame, search path

2007-05-02 Thread Duncan Murdoch
On 01/05/2007 11:34 AM, Prof Brian Ripley wrote: On Tue, 1 May 2007, Duncan Murdoch wrote: On 01/05/2007 12:29 AM, Graham Wideman wrote: [...] Refman p122: Environments consist of a frame, or collection of named objects, and a pointer to an enclosing environment. Is the or here explaining

[R] Degrees of freedom in repeated measures glmmPQL

2007-05-02 Thread Charlotte Burn
Hello, I've just carried out my first good-looking model using glmmPQL, and the output makes perfect sense in terms of how it fits with our hypothesis and the graphical representation of the data. However, please could you clarify whether my degrees of freedom are appropriate? I had 106

Re: [R] Log-likelihood function

2007-05-02 Thread Doxastic
Thanks. I used this and it gave me the same result as the logLik function. The reason I ask is the SAS output gives me a loglik = 1089. R gives me -298.09583. Both for my reduced model. For the saturated (or complex) model, SAS gives me an loglik = 1143. R gives me -298.1993. The problem

[R] I need help

2007-05-02 Thread elyakhlifi mustapha
hello, I need help because I don't understand the syntaxe else how can I write it for example I writed a script to cut missings values and I have errors if(na==length(C)){ + pos=match(0,match(donGeno[[na-1]],donGeno[[na]],nomatch=0)) + for(k in 1:(na-1)) { +

Re: [R] I need help

2007-05-02 Thread Duncan Murdoch
On 02/05/2007 8:11 AM, elyakhlifi mustapha wrote: hello, I need help because I don't understand the syntaxe else how can I write it for example I writed a script to cut missings values and I have errors if(na==length(C)){ +

[R] Odp: I need help

2007-05-02 Thread Petr PIKAL
Hi [EMAIL PROTECTED] napsal dne 02.05.2007 14:11:51: hello, I need help because I don't understand the syntaxe else how can I write it for example I writed a script to cut missings values and I have errors if(na==length(C)){ + pos=match(0,match(donGeno[[na-1]],donGeno[[na]],nomatch=0))

Re: [R] Get the difference of values to their own median value

2007-05-02 Thread Gabor Grothendieck
Use ave. Also its easier to use read.table to read it in and then convert it to a matrix if that's what you want. Input - 29 4.5 1.505713 29 4.6 1.580402 29 4.7 1.656875 29 4.8 1.735054 30 0 0 30 0.1 0.00096108 30 0.2 0.00323831 29 4.5 1.495148 29 4.6 1.568961 29 4.7 1.644467 30 0 0 30 0.1

[R] repeated measures and recurrent events data

2007-05-02 Thread Luis Guillermo Diaz Monroy
Hello, I need a data set on repeated measures and recurrent events data which had been registered on the same subject. This is very important to my thesis work. Thanks a lot, Luis Guillermo Díaz Monroy Profesor Asociado Universidad Nacional de Colombia

Re: [R] ED50 from logistic model with interactions

2007-05-02 Thread Christian Ritz
Hi Kate, try looking at the package 'drc' on CRAN and in particular look at the example in the help page for the dataset 'daphnids' (?daphnids). You can obtain arbitrary ED values with approximate standard errors using the function 'ED'. Christian

[R] how to concatinate the elements of some text vectors cat() or print() ?

2007-05-02 Thread John Kane
I have some comment text taken from a SAS data file. It is stored in two vectors and is difficult to read. I would like to simply concatentate the individual entries and end up with a character vector that give me one line of text per comment. I cannot see how to do this, yet it must be very

[R] NAs introduced by coercion in dist()

2007-05-02 Thread slomascolo
I work with Windows and use R version 2.4.1. I am JUST starting to learn this program... I get this warning message 'NAs introduced by coercion' while trying to build a distance matrix (to be analyzed with NMDS later) from a 336 x 100 data matrix. The original matrix has lots of zeros and no

[R] KS test pvalue estimation using mctest (library truncgof)

2007-05-02 Thread D ANIELLO CLAUDIA \(MPS - 05966\)
Hi, I'm trying to evaluate a Monte Carlo p-value (using truncgof package) on a left truncated sample. From an empirical sample I've estimated a generalized pareto distribution parameters (xi, beta, threshold) (I've used fExtremes pkg). I'm in doubt on what of the following command is the most

[R] Few contributed binaries for Mac OS X?

2007-05-02 Thread Ista Zahn
Dear all, After upgrading to R 2.5.0 on my Intel-based Mac I find that there are only a few binary packages available. I've installed most of the packages I need from source, but am having trouble compiling RGtk2. Will Mac binary packages be available for R 2.5.0? Thanks, Ista Zahn

Re: [R] Degrees of freedom in repeated measures glmmPQL

2007-05-02 Thread Charlotte Burn
Apologies, I made a mistake with my maths. The degrees of freedom look correct, assuming they are the denominator and that glmms work this way. I just under-estimated the number of data points I had. Sorry. Charlotte On 02/05/07, Charlotte Burn [EMAIL PROTECTED] wrote: Hello, I've just

Re: [R] Polar graph of time and tide

2007-05-02 Thread Gabor Grothendieck
On 5/2/07, Jim Lemon [EMAIL PROTECTED] wrote: Alan E. Davis wrote: I have been trying to visualize times of lowest tides, month by month. I have tide predictions with times either in unix time or a text format, and heights in feet or meters. I had been able to derive the clock times of

Re: [R] how to concatinate the elements of some text vectors cat() or print() ?

2007-05-02 Thread Liaw, Andy
Is paste() what you're looking for? Andy From: John Kane I have some comment text taken from a SAS data file. It is stored in two vectors and is difficult to read. I would like to simply concatentate the individual entries and end up with a character vector that give me one line of text

Re: [R] ED50 from logistic model with interactions

2007-05-02 Thread Rubén Roa
Kate Stark wrote: Hi, I was wondering if someone could please help me. I am doing a logistic regression to compare size at maturity between 3 seasons. My model is: fit - glm(Mature ~ Season * Size - 1, family = binomial, data=dat) where Mature is a binary response, 0 for immature, 1 for

Re: [R] NAs introduced by coercion in dist()

2007-05-02 Thread Peter Dalgaard
slomascolo wrote: I work with Windows and use R version 2.4.1. I am JUST starting to learn this program... I get this warning message 'NAs introduced by coercion' while trying to build a distance matrix (to be analyzed with NMDS later) from a 336 x 100 data matrix. The original matrix has

Re: [R] ROracle issues

2007-05-02 Thread Don MacQueen
You may have a confusion between 64 bit and 32 bit versions. (I'm not certain; I haven't studied your transcript in complete detail, but rather just noticed a couple of things.) There are a lot of references in your transcript to 64 bit libraries, e.g. gcc -I/usr/lib64 ... but I think that

[R] RE : add arrows to barchart with groups

2007-05-02 Thread GOUACHE David
Hello and thank you for this first response I'll rephrase what I mean by more general: in a case where I have only 2 levels in groups, but in which I have more than 2 levels for variable s, and in a case where the second bar is not necessarily higher than the first. Thanks in advance. Best

Re: [R] NAs introduced by coercion in dist()

2007-05-02 Thread Silvia Lomascolo
It was suggested that the 'NAs introduced by coercion' message might be warning me that my data are not what they should be. I checked this using str(PeaksMatrix), as suggested, and the data seem to be what I thought they were: 'data.frame': 335 obs. of 127 variables: $ Code : Factor w/

Re: [R] Log-likelihood function

2007-05-02 Thread Robert A LaBudde
At 07:30 AM 5/2/2007, Doxastic wrote: Thanks. I used this and it gave me the same result as the logLik function. The reason I ask is the SAS output gives me a loglik = 1089. R gives me -298.09583. Both for my reduced model. For the saturated (or complex) model, SAS gives me an loglik = 1143.

[R] stringification magic in subset?

2007-05-02 Thread ivo welch
dear R wizards: I am trying to replace subset() with my own version that first checks that each name in the select statement has a corresponding name in the data set. preferably, it would have the same syntax and semantics as subset() otherwise. alas, subset works in interesting ways:

[R] upgrade to 2.5

2007-05-02 Thread Iasonas Lamprianou
Hi I am using R version 2.4.1. How can I upgrade to version 2.5 without having to install all the packages again? Thanks Jason Dr. Iasonas Lamprianou Department of Education The University of Manchester Oxford Road, Manchester M13 9PL, UK Tel. 0044 161 275 3485 [EMAIL PROTECTED]

[R] A Question about knnFinder

2007-05-02 Thread Saptarshi Guha
Hello, Thanks to the author for writing knnFinder. I seem to have problem with a data set found here (http://www.stat.purdue.edu/~sguha/random/ mydata.csv) This is the R-Code ma-read.csv(~/mydata.csv) ma-as.matrix(ma) rownames(ma)-NULL colnames(ma)-NULL #ma is a 390x2 matrix

[R] Query about RODBC to access MySQL from Windows

2007-05-02 Thread lalitha viswanath
Hi I am trying to use RODBC in R installed on Windows to access MySQL database (on a linux box). I set up a DSN and specified this DSN in R as follows library(RODBC); channel - odbcConnect(mysqldsn); RODB Connection 5 Details: case=nochange PORT=3306 Although this seems to connect

Re: [R] stringification magic in subset?

2007-05-02 Thread Sundar Dorai-Raj
ivo welch said the following on 5/2/2007 8:13 AM: dear R wizards: I am trying to replace subset() with my own version that first checks that each name in the select statement has a corresponding name in the data set. preferably, it would have the same syntax and semantics as subset()

Re: [R] NAs introduced by coercion in dist()

2007-05-02 Thread Peter Dalgaard
Silvia Lomascolo wrote: It was suggested that the 'NAs introduced by coercion' message might be warning me that my data are not what they should be. I checked this using str(PeaksMatrix), as suggested, and the data seem to be what I thought they were: 'data.frame': 335 obs. of 127

[R] reference in article

2007-05-02 Thread Tomas Mikoviny
Hi all R positive, does anyone know how to refer R in article? thanks tomas __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and

Re: [R] Read every second line from ASCII file

2007-05-02 Thread Roland Rau
Dear all, I just realized that I forgotten to write some kind of final email for this thread and to thank you for your help. It seems that the recommeneded procedure in such circumstances has three steps: 1) readLines() 2) select the desired lines 3) strsplit() Thanks Ferdinand, Jim, and Paul!

[R] choose.dir

2007-05-02 Thread Antje
Hi all, I have written a R-script under Windows using choose.dir. Now, I have seen that this function is missing at MacOS. Does anybody know an alternative? Antje __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] sorting in barplot

2007-05-02 Thread Romain . Mayor
Thank you for your answer, But I think it doesn't work, because the function bargraph.CI (I use it because barplot doesn't have a standard error option) use two internally function(one for mean and the other for standard error) which use something like tapply (which sort the factors by alphabet).

Re: [R] reference in article

2007-05-02 Thread Doran, Harold
citation() -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Tomas Mikoviny Sent: Wednesday, May 02, 2007 11:44 AM To: r-help@stat.math.ethz.ch Subject: [R] reference in article Hi all R positive, does anyone know how to refer R in article?

Re: [R] reference in article

2007-05-02 Thread Dimitris Rizopoulos
try citation() 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.be/biostat/

[R] Query about finding correlations

2007-05-02 Thread lalitha viswanath
Hi I have a dataframe which has 3 columns of numeric data A,B,C each of which has been obtained independent of the other. We are trying to find out, which of A or B cause C i.e. We are hypothesising that C is the effect and either A or B, not both is the cause. i.e. A causes C and this

Re: [R] stringification magic in subset?

2007-05-02 Thread ivo welch
Beautiful. thank you. /ivo Try: methods(subset) which will point you to subset.data.frame. The latter code will answer your questions. HTH, --sundar __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] reference in article

2007-05-02 Thread Thomas Lumley
On Wed, 2 May 2007, Tomas Mikoviny wrote: Hi all R positive, does anyone know how to refer R in article? Every time you start R it says (in part) R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages

[R] choose.dir

2007-05-02 Thread niederlein-rstat
Hi all, I have written a R-script under Windows using choose.dir. Now, I have seen that this function is missing at MacOS. Does anybody know an alternative? Antje - Kennt man wirklich jeden über 3 Ecken? Die Antworten gibt's bei Yahoo! Clever.

[R] Reenviar: repeated measures and recurrent events data

2007-05-02 Thread Luis Guillermo Diaz Monroy
Luis Guillermo Díaz Monroy Profesor Asociado Universidad Nacional de Colombia---BeginMessage--- Hello, I need a data set on repeated measures and recurrent events data which had been registered on the same subject. This is very important to my thesis work. Thanks a lot, Luis Guillermo Díaz

Re: [R] Read every second line from ASCII file

2007-05-02 Thread Gabor Grothendieck
I typically wouldn't use strsplit but would reread it using read.table and textConnection as in: http://www.mail-archive.com/r-help@stat.math.ethz.ch/msg84752.html On 5/2/07, Roland Rau [EMAIL PROTECTED] wrote: Dear all, I just realized that I forgotten to write some kind of final email for

Re: [R] choose.dir

2007-05-02 Thread Gabor Grothendieck
I don't have a Mac but perhaps this works: library(tcltk) tclvalue(tkchooseDirectory()) On 5/2/07, Antje [EMAIL PROTECTED] wrote: Hi all, I have written a R-script under Windows using choose.dir. Now, I have seen that this function is missing at MacOS. Does anybody know an alternative?

Re: [R] sorting in barplot

2007-05-02 Thread Marc Schwartz
Romain, Try this before calling bargraph.CI(): Fst2$REGION - reorder(Fst2$REGION, Fst2$MOY_FST, mean) This will reorder the factor levels for REGION, based upon the mean value of MOY_FST for each region. Thus, your REGIONS should then be sorted by increasing mean values in the plot. For

Re: [R] Query about finding correlations

2007-05-02 Thread Lalitha Viswanath
Hi This is not a homework assignment :) Me and my manager are trying to understand the problem better. In the meanwhile, we thought we would post the problem on this forum to seek some input from statisticians who possibly do this kind of analyses everyday and hence are possibly more proficient

Re: [R] Query about finding correlations

2007-05-02 Thread Stefan Grosse
I am sorry if I misinterpreted this but this questions looks/looked very like a very basic one about statistical inference... (have a look at the Manuals/contributed documentation e.g. Faraway) Lalitha Viswanath wrote: Hi This is not a homework assignment :) Me and my manager are trying to

Re: [R] reference in article

2007-05-02 Thread Peter Dalgaard
Tomas Mikoviny wrote: Hi all R positive, does anyone know how to refer R in article? thanks tomas Like this (well, almost! This is from R 2.4.1 - 2007 on the most recent version) citation() To cite R in publications use: R Development Core Team (2006). R: A language and

[R] Multiple scatterplots

2007-05-02 Thread Kostadin Cholakov
Hi, I have to plot three Ziph distributions for three languages where the x value represents the rank of a given word and the y value represents the relative frequency of this word in the corpus. Is there some way so that I can plot all three distributions on a single scatterplot, preferably with

Re: [R] Query about finding correlations

2007-05-02 Thread Stefan Grosse
How about making your homeworks yourselfes? lalitha viswanath wrote: Hi I have a dataframe which has 3 columns of numeric data A,B,C each of which has been obtained independent of the other. We are trying to find out, which of A or B cause C i.e. We are hypothesising that C is the effect

[R] Multiple scatterplots

2007-05-02 Thread Kostadin Cholakov
Hi, I have to plot three Ziph distributions for three languages where the x value represents the rank of a given word and the y value represents the relative frequency of this word in the corpus. Is there some way so that I can plot all three distributions on a single scatterplot, preferably with

Re: [R] upgrade to 2.5

2007-05-02 Thread Sundar Dorai-Raj
Iasonas Lamprianou said the following on 5/2/2007 8:25 AM: Hi I am using R version 2.4.1. How can I upgrade to version 2.5 without having to install all the packages again? Thanks Jason You may find the following link relevant.

Re: [R] choose.dir

2007-05-02 Thread Antje
Thank you. For Windows it works, I'll check whether it also works for Mac :) Gabor Grothendieck schrieb: I don't have a Mac but perhaps this works: library(tcltk) tclvalue(tkchooseDirectory()) On 5/2/07, Antje [EMAIL PROTECTED] wrote: Hi all, I have written a R-script under Windows

Re: [R] Query about finding correlations

2007-05-02 Thread Alberto Monteiro
Lalitha Viswanath wrote: We are trying to find out, which of A or B cause C i.e. We are hypothesising that C is the effect and either A or B, not both is the cause. (...) I would greatly appreciate any inputs on the best statistcal approach to tackle this problem. I am thinking that we

Re: [R] upgrade to 2.5

2007-05-02 Thread Gabor Grothendieck
On 5/2/07, Sundar Dorai-Raj [EMAIL PROTECTED] wrote: Iasonas Lamprianou said the following on 5/2/2007 8:25 AM: Hi I am using R version 2.4.1. How can I upgrade to version 2.5 without having to install all the packages again? Thanks Jason You may find the following link relevant.

Re: [R] Query about finding correlations

2007-05-02 Thread Lalitha Viswanath
Hi Thanks for your input. I stand corrected. The causation is not linear. We wish to find out which is the cause and under what circumstances. i.e. at what points along a scale for C, is A the cause and when does B become the cause, if at all. As a crude analyses, we assumed that the above is not

Re: [R] Multiple scatterplots

2007-05-02 Thread John Kane
Your title and your posting do not say the same thing. Assuming you want all three distributions on one scatter plot does this help? aa - 1:10 bb - 11:2 cc - bb^2 dd - c(3,4,7,9,11,32,11,14,5,9) plot(aa,cc, col=red) points(aa,bb, col=blue) points(aa,dd, col=green) Also in plotting it is a

Re: [R] Free Webinar: Vendor Neutral Intro to Data Mining for AbsoluteBeginners, May 23, 2007

2007-05-02 Thread Brian Koch
Lisa: Can we expect to see R used [exclusively, I would hope] during your demonstration? Learning how data mining models work: the inputs, the outputs, and the nature of the predictive mechanism only makes sense for me if I can follow/retrace your steps on my systems. Thank you. Brian J. Koch

[R] How to install previous packages after upgrading to R 2.5.0?

2007-05-02 Thread Jeffrey Wood
Hello, I have just upgraded from R-2.4.1 to R-2.5.0 for Windows. I had installed a large number of add-on packages under 2.4.1. Is there an easy way to install (or load, if that's the easier way) those packages under 2.5.0, without having to install each package by hand? Thanks, Jeff

Re: [R] R GUI in Ubuntu Feisty?

2007-05-02 Thread duggula
Hello Andy, I had the same problem in Feisty and now I fixed it. I updated to Java 6 through the repositories from Feisty. Then I choosed java 6 with sudo update-alternatives --config java. After that I started R with sudo R. If you had already installed JGR, but it doesn't run you have to

Re: [R] upgrade to 2.5

2007-05-02 Thread Robert A LaBudde
At 01:41 PM 5/2/2007, you wrote: On 5/2/07, Sundar Dorai-Raj [EMAIL PROTECTED] wrote: Iasonas Lamprianou said the following on 5/2/2007 8:25 AM: Hi I am using R version 2.4.1. How can I upgrade to version 2.5 without having to install all the packages again? Thanks Jason

Re: [R] How to install previous packages after upgrading to R 2.5.0?

2007-05-02 Thread Duncan Murdoch
On 02/05/2007 5:12 PM, Jeffrey Wood wrote: Hello, I have just upgraded from R-2.4.1 to R-2.5.0 for Windows. I had installed a large number of add-on packages under 2.4.1. Is there an easy way to install (or load, if that's the easier way) those packages under 2.5.0, without having

Re: [R] upgrade to 2.5

2007-05-02 Thread Sundar Dorai-Raj
Robert A LaBudde said the following on 5/2/2007 2:39 PM: At 01:41 PM 5/2/2007, you wrote: On 5/2/07, Sundar Dorai-Raj [EMAIL PROTECTED] wrote: Iasonas Lamprianou said the following on 5/2/2007 8:25 AM: Hi I am using R version 2.4.1. How can I upgrade to version 2.5 without having to

Re: [R] How to install previous packages after upgrading to R 2.5.0?

2007-05-02 Thread John C Frain
Rename the library subdirectory in your new installation to something like library.new. Copy the library subdirectory in the old distribution to the new installation. Copy the contents of the renamed library.new to the the directory copied from the old distribution. You should then start your

Re: [R] Multiple scatterplots

2007-05-02 Thread Bert Gunter
Please note: in R you can specify (some of the) graphics parameters as the appropriate length vectors. So your plot example below can also be done as, for example: plot( rep.int(aa,3),c(cc,bb,dd),col=rep(c(red,blue,green),e=length(aa))) However, this doesn't seem to fit the posted request, where

[R] Is R's fast fourier transform function different from fft2 in Matlab?

2007-05-02 Thread Li Li
Hi All, I found mvfft in R and fft2 in Matlab give different result and can't figure out why. My example is: In R: matrix(c(1,4,2,20), nrow=2) [,1] [,2] [1,]12 [2,]4 20 mvfft(matrix(c(1,4,2,20), nrow=2)) [,1] [,2] [1,] 5+0i 22+0i [2,] -3+0i -18+0i In Matlab:

[R] how to reproduce the same sampled units?

2007-05-02 Thread Santanu Pramanik
Hi all, Is it possible to generate the same sample number of times in R? In SAS, using the option seed it is possible to reproduce exactly the same sample. Is there any such feature in R which I can use? For further clarity, for (i in 1:2) { samp = sample(1:1000,100,replace = FALSE)

Re: [R] how to reproduce the same sampled units?

2007-05-02 Thread Ted Harding
On 02-May-07 23:25:18, Santanu Pramanik wrote: Hi all, Is it possible to generate the same sample number of times in R? In SAS, using the option seed it is possible to reproduce exactly the same sample. Is there any such feature in R which I can use? For further clarity, for (i in

Re: [R] upgrade to 2.5

2007-05-02 Thread Gabor Grothendieck
On 5/2/07, Robert A LaBudde [EMAIL PROTECTED] wrote: At 01:41 PM 5/2/2007, you wrote: On 5/2/07, Sundar Dorai-Raj [EMAIL PROTECTED] wrote: Iasonas Lamprianou said the following on 5/2/2007 8:25 AM: Hi I am using R version 2.4.1. How can I upgrade to version 2.5 without having to

Re: [R] How to install previous packages after upgrading to R 2.5.0?

2007-05-02 Thread Gabor Grothendieck
On 5/2/07, Jeffrey Wood [EMAIL PROTECTED] wrote: Hello, I have just upgraded from R-2.4.1 to R-2.5.0 for Windows. I had installed a large number of add-on packages under 2.4.1. Is there an easy way to install (or load, if that's the easier way) those packages under 2.5.0, without

Re: [R] Is R's fast fourier transform function different from fft2 in Matlab?

2007-05-02 Thread Sundar Dorai-Raj
Li Li said the following on 5/2/2007 4:06 PM: Hi All, I found mvfft in R and fft2 in Matlab give different result and can't figure out why. My example is: In R: matrix(c(1,4,2,20), nrow=2) [,1] [,2] [1,]12 [2,]4 20 mvfft(matrix(c(1,4,2,20), nrow=2)) [,1]

[R] about using read.table

2007-05-02 Thread hassen62
Hi,Dear R users, I have a file text nommed chif which contains 16 lines and 4 columns in the disc dur. I have a difficulty to read this file in R console I have used the following command chif - read.table(c:/chif.txt, header=T, sep= ) I have obtained from R console:undefined file! Can you

[R] Building package: What does this message about rcompgen imply?

2007-05-02 Thread Peter Dunn
Hi all I'm updating a package for submission to CRAN, and I am getting an error message I never have seen before, and can't find out what it implies or means (and hence, what I need to do to fix it). After running R CMD check, I get this: snip * checking foreign function calls ... OK * checking

Re: [R] how to reproduce the same sampled units?

2007-05-02 Thread Kenneth Cabrera
?set.seed On Wed, 02 May 2007 18:25:18 -0500, Santanu Pramanik [EMAIL PROTECTED] wrote: Hi all, Is it possible to generate the same sample number of times in R? In SAS, using the option seed it is possible to reproduce exactly the same sample. Is there any such feature in R which I can

Re: [R] Building package: What does this message about rcompgen imply?

2007-05-02 Thread Prof Brian Ripley
We don't have any of the basic details the posting guide asks for! If this is R 2.5.0, rcompgen is a recommended package and ought to be installed. If it is not, then .check_Rd_refs cannot cross-check references against all the standard and recommended packages, and that looks like the cause

Re: [R] Building package: What does this message about rcompgen imply?

2007-05-02 Thread Dirk Eddelbuettel
Peter, On 3 May 2007 at 10:43, Peter Dunn wrote: | After running R CMD check, I get this: | | snip | * checking foreign function calls ... OK | * checking R code for possible problems ... OK | * checking Rd files ... OK | * checking Rd cross-references ... WARNING | Error in

Re: [R] Building package: What does this message about rcompgen imply?

2007-05-02 Thread Peter Dunn
Dirk On 3 May 2007 at 10:43, Peter Dunn wrote: | After running R CMD check, I get this: | | snip | * checking foreign function calls ... OK | * checking R code for possible problems ... OK | * checking Rd files ... OK | * checking Rd cross-references ... WARNING | Error in

[R] Survival statistics--displaying multiple plots

2007-05-02 Thread Gregory Pierce
Hello all! I am once again analyzing patient survival data with chronic liver disease. The severity of the liver disease is given by a number which is continuously variable. I have referred to this number as meld--model for end stage liver disease--which is the result of a mathematical

Re: [R] Building package: What does this message about rcompgen imply?

2007-05-02 Thread Dirk Eddelbuettel
Peter, On 3 May 2007 at 12:23, Peter Dunn wrote: | Are you by chance on Debian (or Ubuntu)? | | Bingo. | | The short answer seems to be that I do need rcompgen to be | installed, and my Debian system didn't do that. Easy fixed though. Yes, sorry. I realized relatively late that we needed

Re: [R] Is R's fast fourier transform function different from fft2 in Matlab?

2007-05-02 Thread Li Li
Thanks for both replies. Then I found the ifft2 from Matlab gives different result from fft( , inverse=T) from R. An example: in R: temp - matrix(c(1,4,2, 20), nrow=2) fft(temp) [,1] [,2] [1,] 27+0i -17+0i [2,] -21+0i 15+0i fft(temp,inverse=T) [,1] [,2] [1,] 27+0i -17+0i

Re: [R] Survival statistics--displaying multiple plots

2007-05-02 Thread Gregory Pierce
I should clarify. I can generate plots for each category individually but not for all three on the same chart. Greg -Original Message- From: Gregory Pierce [mailto:[EMAIL PROTECTED] Sent: Wednesday, May 02, 2007 10:21 PM To: 'r-help@stat.math.ethz.ch' Subject: Survival

Re: [R] Is R's fast fourier transform function different from fft2 in Matlab?

2007-05-02 Thread Robert A LaBudde
Discrete Fourier transforms can be normalized in different ways. Some apply the whole normalization to the forward transform, some to the reverse transform, some apply the square root to each, and some don't normalize at all (in which case the reverse of the forward transform will need

Re: [R] Survival statistics--displaying multiple plots

2007-05-02 Thread Robert A LaBudde
? layout() ? par() E.g., layout(matrix(c(1,2,3),1,3,byrow=TRUE) #3 plots side-by-side Then use plot() three times to generate each of your graphs. At 11:14 PM 5/2/2007, Greg wrote: I should clarify. I can generate plots for each category individually but not for all three on the same chart.

Re: [R] Is R's fast fourier transform function different from fft2 in Matlab?

2007-05-02 Thread Sundar Dorai-Raj
Li Li said the following on 5/2/2007 7:53 PM: Thanks for both replies. Then I found the ifft2 from Matlab gives different result from fft( , inverse=T) from R. An example: in R: temp - matrix(c(1,4,2, 20), nrow=2) fft(temp) [,1] [,2] [1,] 27+0i -17+0i [2,] -21+0i 15+0i

Re: [R] Survival statistics--displaying multiple plots

2007-05-02 Thread Petr Klasterecky
Have a look at ?lines ?points ?plot - the option add, add=TRUE You will have to specify proper limits on both axes, otherwise you may only see parts of the graphs. Some functions in the survival library also allow stratified analyses in which case the plots account for different strata. Petr

  1   2   >