Re: [R] Histograms of lots of variables

2007-04-17 Thread Petr Klasterecky
Please, experiment a little bit :-) I have trouble using Win metafile graphics, so I just converted it to jpeg. The code works - there was a copy/paste error. Try this with your 'total' dataframe - it works with my 'total' given there: total -

Re: [R] Please help.... i know its trivial

2007-04-17 Thread Petr Klasterecky
Patrick Callier napsal(a): rownames is the command you want, I think. rownames(A) - seq(1:8) Yes. More generally rownames(A) - 1:length(A[,1]) And keep in mind that nobody is going to help you next time, if you won't give a clear and informative SUBJECT in your message. Here it should

Re: [R] Setting rownames (was Please help.... i know its trivial)

2007-04-17 Thread Prof Brian Ripley
On Tue, 17 Apr 2007, Petr Klasterecky wrote: Patrick Callier napsal(a): rownames is the command you want, I think. rownames(A) - seq(1:8) Yes. More generally rownames(A) - 1:length(A[,1]) Better rownames(A) - seq_len(nrow(A)) which works correctly if the length is 0 or there are zero

Re: [R] Greek symbols in xtable rows

2007-04-17 Thread Andris Jankevics
Thank You for answer. I Tried code that You provided in reply to my question on both PC's with Linux and Widows OS. On linux box output is: xtable(diag(c($\\sigma_1^2$,$\\sigma_2^2$))) % latex table generated in R 2.4.0 by xtable 1.4-3 package % Tue Apr 17 09:18:31 2007 \begin{table}[ht]

Re: [R] It is possible to Expand RAM with HD space?

2007-04-17 Thread Uwe Ligges
Jorge Cornejo-Donoso wrote: I'm using R in a winXp machine. I have problems with the size of my database (xyz), anyone know if is possible to use some kind of swap memory, in order to expand my real RAM and use HD space for the processes? This is done automatically by your operating system.

[R] Extracting approximate Wald test (Chisq) from coxph(..frailty)

2007-04-17 Thread Mohammad Ehsanul Karim
Dear List, How do I extract the approximate Wald test for the frailty (in the following example 17.89 value)? What about the P-values, other Chisq, DF, se(coef) and se2? How can they be extracted? ## kfitm1 Call: coxph(formula = Surv(time,

Re: [R] colSum() in Matrix objects

2007-04-17 Thread Martin Maechler
Hi Jose, Jose I'd like to simply add column-wise using Matrix objects (Csparse). Jose It looks like one can apply mosty any base function to these objects Jose (i.e., apply, colSums), but there is a nasty conversion to traditional Jose matrix objects if one does that. not

[R] background color

2007-04-17 Thread yannig goude
hi, I want to add different colors on the background of a classical plot. Each color is associated to an interval of the x axis. example: the background is red on the interval [1,10], blue on [11,20]. I try the rect function but it isn't appropriate for the background. Can any one can

Re: [R] It is possible to Expand RAM with HD space?

2007-04-17 Thread Prof Brian Ripley
On Tue, 17 Apr 2007, Uwe Ligges wrote: Jorge Cornejo-Donoso wrote: I'm using R in a winXp machine. I have problems with the size of my database (xyz), anyone know if is possible to use some kind of swap memory, in order to expand my real RAM and use HD space for the processes? This is

[R] help comparing two median with R

2007-04-17 Thread Pedro A Reche
Dear R users, I am new to R and I would like to ask your help with the following topic. I have three sets of numeral data, 2 sets are paired and a third is independent of the other two. For each of these sets I have obtained their basic statistics (mean, median, stdv, range ...). Now I

Re: [R] background color

2007-04-17 Thread Petr Klasterecky
What do you mean by background? Maybe this is enough: plot(seq(-3,3,.01), dnorm(seq(-3,3,.01)), type=n, xlab=x, ylab=f(x), main=Normal density) polygon(x=c(-4,0,0,-4), y=c(-1,-1,.5,.5), col=red) polygon(x=c(4,0,0,4), y=c(-1,-1,.5,.5), col=blue) lines(seq(-3,3,.01), dnorm(seq(-3,3,.01)), type=l,

[R] Use of argument '...'

2007-04-17 Thread Juan Lewinger
Dear R list, I've read the function writing sections on both An introduction to R and R language Definition manuals but still don't understand why the following gives an error message: fun - function(x, ...) x + y fun(1, y=2) I get: Error in fun(1, y = 2) : object y not found I'd appreciate

Re: [R] Use of argument '...'

2007-04-17 Thread Soare Marcian-Alin
Dear Juan, You have to define your function like this: fun - function (x,y=0, ...) { x+y } The ... means that u can use the other parameters in par! The variables, which you want to use, these you have to define in your function! In this function i put y=0 on a standard value zero, if

Re: [R] Use of argument '...'

2007-04-17 Thread Peter Dalgaard
Juan Lewinger wrote: Dear R list, I've read the function writing sections on both An introduction to R and R language Definition manuals but still don't understand why the following gives an error message: fun - function(x, ...) x + y fun(1, y=2) I get: Error in fun(1, y = 2) :

Re: [R] help comparing two median with R

2007-04-17 Thread gyadav
hi +++ Now I want to compare if these sets differ. I could compare the mean doing a basic T test . However, I was looking for a test to compare the medians using R. +++ @@@i hope but i am not sure that wilcoxon signed rank test may help :-) cheers Pedro A Reche [EMAIL PROTECTED]

[R] СНИЖЕНИЕ ВЕСА

2007-04-17 Thread Guillermo Berry
ÇÎËÎÒÀß ÑÅÐÜÃÀ - ÈÇÁÀÂËÅÍÈÅ ÎÒ ËÈØÍÈÕ ÊÈËÎÃÐÀÌÌΠ(4 Êà  ÌÅÑ.)!! ÎÃÐÎÌÍÎÅ ÊÎËÈ×ÅÑÒÂÎ ËÞÄÅÉ ÂÎ ÂÑÅÌ ÌÈÐÅ ÏÎÕÓÄÅËÈ ÁËÀÃÎÄÀÐß ÇÎËÎÒÎÉ ÑÅÐÜÃÅ ÁÅÇ ÈÑÒßÇÀÞÙÈÕ ÄÈÅÒ, ÔÈÇÈ×ÅÑÊÈÕ ÓÏÐÀÆÍÅÍÈÉ, ÒÀÁËÅÒÎÊ! ÑÒÎÈÌÎÑÒÜ ÂÑÅÃÎ ËÈØÜ ÇÎÎÎ ÐÓÁËÅÉ!!! Ïðèíöèï ìåòîäà «Çîëîòàÿ Ñåðüãà» ñîçäàí íà îñíîâå

Re: [R] help comparing two median with R

2007-04-17 Thread Jim Lemon
Pedro A Reche wrote: Dear R users, I am new to R and I would like to ask your help with the following topic. I have three sets of numeral data, 2 sets are paired and a third is independent of the other two. For each of these sets I have obtained their basic statistics (mean, median,

[R] About PLR

2007-04-17 Thread Nitish Kumar Mishra
Hello Sir/Mam I am Nitish Kumar Mishra from IMTECH, Chandigarh, India. I want to calculate the Principal Component Analysis(PCA), PLS of the input file and find top 30 PCAs of this file using PLS in R. I am asking regarding Linux(Red Hat 9) operating system. I have downloaded PLS from Crains site

Re: [R] Find zeros of analytic functions

2007-04-17 Thread Alberto Monteiro
Robin Hankin wrote: If iterative methods are appropriate, it's perhaps worth pointing out that Newton-Rapheson works nicely for complex functions. Hmmm... I think there are many cases where Newton-Raphson diverges for complex functions, like those that generate beautiful fractals. Alberto

Re: [R] help comparing two median with R

2007-04-17 Thread Frank E Harrell Jr
Thomas Lumley wrote: On Tue, 17 Apr 2007, Robert McFadden wrote: -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Jim Lemon Sent: Tuesday, April 17, 2007 12:37 PM To: Pedro A Reche Cc: r-help@stat.math.ethz.ch Subject: Re: [R] help comparing two

Re: [R] extracting intercept from ppr fit

2007-04-17 Thread Vadim Ogranovich
Hi, Is there a way, documented or not, to extract the intercept term (the alpha_0 the MASS book) from a ppr() (Projection Persuit Regression) fit? Thanks, Vadim ## Example: n - 1000 data - data.frame(x=rnorm(n), y=rnorm(n)) a - 10 data$z - evalq(a + atan(x + y) + rnorm(n), data)

Re: [R] Modelling Heteroscedastic Multilevel Models

2007-04-17 Thread Doran, Harold
I think there are many who can help, but this question is quite vague. This assumes we have access to the book you note and can make sense of your question w/o sample data. If you cannot find a sample data set please create a sample data file. However, there are so many sample data sets in the

[R] extracting intercept from ppr fit

2007-04-17 Thread Vadim Ogranovich
Hi, Is there a way, documented or not, to extract the intercept term (the alpha_0 the MASS book) from a ppr() (Projection Persuit Regression) fit? Thanks, Vadim [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing

Re: [R] About PLR

2007-04-17 Thread Bjørn-Helge Mevik
From within R, you can give the command install.packages(pls) and R will download and install it for you (as long as you have access to the Internet). To install an already downloaded package, you can use R CMD INSTALL pls_2.0-0.tar.gz in a terminal window. -- Bjørn-Helge Mevik

Re: [R] help comparing two median with R

2007-04-17 Thread Prof Brian Ripley
On Tue, 17 Apr 2007, Robert McFadden wrote: -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Jim Lemon Sent: Tuesday, April 17, 2007 12:37 PM To: Pedro A Reche Cc: r-help@stat.math.ethz.ch Subject: Re: [R] help comparing two median with R Pedro

Re: [R] background color

2007-04-17 Thread Gabor Grothendieck
On 4/17/07, yannig goude [EMAIL PROTECTED] wrote: I want to add different colors on the background of a classical plot. Each color is associated to an interval of the x axis. example: the background is red on the interval [1,10], blue on [11,20]. I try the rect function but it isn't

[R] sample size 2-sample proportion test

2007-04-17 Thread Berta
Hi R-users, I want to calculate the sample size needed to carry out a 2-sample proprotion test (with alfa=0.05, beta=0.8) 1.- probability of success in subpopulation A: 0.8 2.- probability of success in subpopulation B: 0.05 3.- percentage of population in subpopulation A = 5%, 4.- percentage

[R] Special Limited Time Discount Offer from Chapman Hall/CRC Press

2007-04-17 Thread English, Nadja
EXTRA SPECIAL 50% OFF DISCOUNT ON SELECT TITLES including: Robust Statistical Methods with R by Jana Jureckova and Jan Pickek List Price: $84.95 SALE PRICE: $42.48 An Introduction to the Bootstrap by Bradley Efron and R.J. Tibshirani List Price: $99.95 SALE PRICE:

[R] Inverse of one function

2007-04-17 Thread enrico\.foscolo2
Dear Partecipants to the list, I am searching a R function, which can calculate the inverse of one real value function. Does it exists a R code in order to make it? Many thanks for any kind of help and for Your availability. Enrico Foscolo

Re: [R] My First Function: cryptic error message

2007-04-17 Thread Joel J. Adamson
Weiwei Shi writes: I cannot reproduce the error w/o the data called probands I realized this when I posted; sharing the data with anyone is not possible (I hope you understand). I'm hoping a more experienced useR can spot a syntax or procedural error that I made. I have made some progress.

[R] A problem about all possible sequences

2007-04-17 Thread Paul Smith
Dear All Suppose a sequence of length 10 generated by the following rule: the first element is 0 with 50% of probability or 1 with the same probability; the second element likewise; and so on. Is there some R command to obtain all possible different sequences formed by the above rule? I am aware

Re: [R] A problem about all possible sequences

2007-04-17 Thread Paul Smith
On 4/17/07, Robin Hankin [EMAIL PROTECTED] wrote: f - function(n){expand.grid(rep(list(0:1),n))} f(10) [snip] Thanks, Robin, that is it! Paul On 17 Apr 2007, at 15:26, Paul Smith wrote: Dear All Suppose a sequence of length 10 generated by the following rule: the first element is

[R] my.symbols and zoomplot functions now available in TeachingDemos package

2007-04-17 Thread Greg Snow
In February (2/26) Jonathan Lees asked a question about plotting symbols whose shape/aspect ratio is independent of of the user coordinate system. A couple of other posts have also been made asking about plotting symbols of a given fixed size/aspect ratio. The latest version of the

Re: [R] A problem about all possible sequences

2007-04-17 Thread Robin Hankin
Paul f - function(n){expand.grid(rep(list(0:1),n))} f(10) [snip] HTH Robin On 17 Apr 2007, at 15:26, Paul Smith wrote: Dear All Suppose a sequence of length 10 generated by the following rule: the first element is 0 with 50% of probability or 1 with the same probability; the second

Re: [R] PROC DISCRIM vs. lda( ) in MASS

2007-04-17 Thread Prof Brian Ripley
On Tue, 17 Apr 2007, Greg Tarpinian wrote: Hello, I am using WinXP, R version 2.3.1, and SAS for PC version 8.1. I have mostly used SAS over the last 4 years and would like to compare the output of PROC DISCRIM to that of lda( ) with respect to a very specific aspect. My data have k=3

Re: [R] help comparing two median with R

2007-04-17 Thread Thomas Lumley
On Tue, 17 Apr 2007, Robert McFadden wrote: -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Jim Lemon Sent: Tuesday, April 17, 2007 12:37 PM To: Pedro A Reche Cc: r-help@stat.math.ethz.ch Subject: Re: [R] help comparing two median with R Pedro A

Re: [R] Export multiple data files from R

2007-04-17 Thread Vladimir Eremeev
You didn't describe the exact format of the .gpr files. There are 32 heading lines in each file, which are now hidden from R community. You 'skip' 31 of them in read.table, and one more plays header ('header=T'). Since you are using read.table, your files are usual ascii files. You shoud use

Re: [R] read.spss (package foreign) and SPSS 15.0 files

2007-04-17 Thread John Kane
--- Charilaos Skiadas [EMAIL PROTECTED] wrote: On Apr 16, 2007, at 10:41 AM, John Kane wrote: --- Charilaos Skiadas [EMAIL PROTECTED] wrote: It is not an export option, it is a save as option. I don't have a 14 to check, but on a 15 I just go to File - Save As, and change the

Re: [R] A problem about all possible sequences

2007-04-17 Thread Ted Harding
On 17-Apr-07 14:26:15, Paul Smith wrote: Dear All Suppose a sequence of length 10 generated by the following rule: the first element is 0 with 50% of probability or 1 with the same probability; the second element likewise; and so on. You don't say whether the elements of the sequence are

Re: [R] Greek symbols in xtable rows

2007-04-17 Thread James W. MacDonald
Roger Bivand wrote: Andris Jankevics andza at osi.lv writes: Thank You for answer. I Tried code that You provided in reply to my question on both PC's with Linux and Widows OS. On linux box output is: ... System information: sessionInfo () R version 2.4.0 (2006-10-03)

Re: [R] Greek symbols in xtable rows

2007-04-17 Thread Roger Bivand
On Tue, 17 Apr 2007, James W. MacDonald wrote: Roger Bivand wrote: Andris Jankevics andza at osi.lv writes: Thank You for answer. I Tried code that You provided in reply to my question on both PC's with Linux and Widows OS. On linux box output is: ... System information:

Re: [R] read.spss (package foreign) and SPSS 15.0 files

2007-04-17 Thread Frank E Harrell Jr
John Kane wrote: --- Charilaos Skiadas [EMAIL PROTECTED] wrote: On Apr 16, 2007, at 10:41 AM, John Kane wrote: --- Charilaos Skiadas [EMAIL PROTECTED] wrote: It is not an export option, it is a save as option. I don't have a 14 to check, but on a 15 I just go to File - Save As, and

[R] GREP - Choosing values between two borders

2007-04-17 Thread Felix Wave
Hello, I import datas from an file with: readLines But I need only a part of all measurments of this file. These are between two borders START and END. Can you tell me the syntax of grep(), to choose values between two borders? My R Code was not succesful, and I can't finde anything in the help.

Re: [R] installing R on Ubuntu

2007-04-17 Thread Roland Rau
I am not sure what the question is but maybe this provides a starting point? http://cran.r-project.org/bin/linux/ubuntu/README HTH, Roland On 4/16/07, Erin Hodgess [EMAIL PROTECTED] wrote: __ R-help@stat.math.ethz.ch mailing list

Re: [R] Extracting approximate Wald test (Chisq) from coxph(..frailty)

2007-04-17 Thread Paul Artes
Assign the output of coxph to some object, and use the $ extractor function to obtain what you need. ie: rtfm - coxph(formula = Surv(time, status) ~ age + sex + disease + frailty(id, dist = gauss), data = kidney) Age - coef(rtfm)[age] OR Sex - rtfm$coef[sex] Hope this helps. Paul Mohammad

Re: [R] Extracting approximate Wald test (Chisq) from coxph(..frailty)

2007-04-17 Thread Charles C. Berry
On Tue, 17 Apr 2007, Mohammad Ehsanul Karim wrote: Dear List, How do I extract the approximate Wald test for the frailty (in the following example 17.89 value)? The example you give silently invokes print.coxph() to produce that output. You _can_ use tmp - capture.output( print(

Re: [R] A problem about all possible sequences

2007-04-17 Thread Chuck Cleland
Paul Smith wrote: Dear All Suppose a sequence of length 10 generated by the following rule: the first element is 0 with 50% of probability or 1 with the same probability; the second element likewise; and so on. Is there some R command to obtain all possible different sequences formed by

[R] PROC DISCRIM vs. lda( ) in MASS

2007-04-17 Thread Greg Tarpinian
Hello, I am using WinXP, R version 2.3.1, and SAS for PC version 8.1. I have mostly used SAS over the last 4 years and would like to compare the output of PROC DISCRIM to that of lda( ) with respect to a very specific aspect. My data have k=3 populations and there are 3 variates in the feature

Re: [R] Greek symbols in xtable rows

2007-04-17 Thread Roger Bivand
Andris Jankevics andza at osi.lv writes: Thank You for answer. I Tried code that You provided in reply to my question on both PC's with Linux and Widows OS. On linux box output is: ... System information: sessionInfo () R version 2.4.0 (2006-10-03) i686-redhat-linux-gnu

Re: [R] help comparing two median with R

2007-04-17 Thread Robert McFadden
-Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Jim Lemon Sent: Tuesday, April 17, 2007 12:37 PM To: Pedro A Reche Cc: r-help@stat.math.ethz.ch Subject: Re: [R] help comparing two median with R Pedro A Reche wrote: Dear R users, I am new

[R] Tcltk

2007-04-17 Thread Sofia Wikström
I have problems with Tcl/Tk in R 2.4.1, when running it on Windows Vista (see error message below). Regards, Sofia library(tcltk) Loading Tcl/Tk interface ... Error in fun(...) : Can't find a usable init.tcl in the following directories: {C:\Program\R\R-2.4.1/Tcl/lib/tcl8.4}

Re: [R] Modelling Heteroscedastic Multilevel Models

2007-04-17 Thread Douglas Bates
May I suggest the Exam data from the mlmRev package as an example. If you wish to have a random effect for Sex by school you could write the model as lmer(test.result ~ homework + Sex + (Sex|school)) which gives correlated random effects for the overall achievement in schools and the

Re: [R] GREP - Choosing values between two borders

2007-04-17 Thread Gabor Grothendieck
You can adapt this to your situation: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/22195.html On 4/17/07, Felix Wave [EMAIL PROTECTED] wrote: Hello, I import datas from an file with: readLines But I need only a part of all measurments of this file. These are between two borders START and

Re: [R] langage R

2007-04-17 Thread R. Villegas
2007/4/16, elyakhlifi mustapha [EMAIL PROTECTED]: bonjour, je vous écris pour savoir si vous aviez trouver réponse pour l'écriture du test de Newman Keuls sous R car j'en ai aussi besoin merci encore cordialement. - [[alternative HTML

[R] Greatest common divisor of two numbers

2007-04-17 Thread AbouEl-Makarim Aboueissa
Dear Sir/Madam: Could you please let me know which function shall I use to get the Greatest common divisor of two numbers. Thank you so much for your attention to this matter, and i look forward to hear from you soon. Regards; Abou == AbouEl-Makarim Aboueissa, Ph.D.

Re: [R] Greatest common divisor of two numbers

2007-04-17 Thread Patrick Burns
S Poetry has a function for that. Patrick Burns [EMAIL PROTECTED] +44 (0)20 8525 0696 http://www.burns-stat.com (home of S Poetry and A Guide for the Unwilling S User) AbouEl-Makarim Aboueissa wrote: Dear Sir/Madam: Could you please let me know which function shall I use to get the Greatest

[R] how to estimate dose from respond given drc package result

2007-04-17 Thread John McNeil
Dear all, I can use the very nice drc package (multdrc()) to model and plot a dataframe containing dose and response values. I can also use predict.drc() to yield response values given a dose. I need to do the opposite, estimate a dose given the response. The general predict documentation

[R] about contrasts

2007-04-17 Thread Cleber Borges
Hi all R-users Please, I would like to learn about 'contrasts'. I do not know much about importance and several types (of contrasts) over regression and data analysis. I would like references for my study! Thanks in advanced... klebyn

Re: [R] Greatest common divisor of two numbers

2007-04-17 Thread Roland Rau
Hi, wouldn't it be enough for simple purposes (as stated in the original message: Greatest common divisor of two numbers) to use: gcd - function(a,b) ifelse (b==0, a, gcd(b, a %% b)) gcd(12,4) [1] 4 gcd(4,12) [1] 4 gcd(123456789,987654321) [1] 9 gcd(987654321,123456789) [1] 9 Best, Roland

[R] value of complexity parameter in ridge regression

2007-04-17 Thread sikander hayat
Hi, What is the optimum range to look for a value of lambda while doing ridge regression. Can/ should lambda be greater than 1 ? I have conflicting (or what appears conflicting to me) sources that use lambda = 0, without any upper limit, but that makes the search space infinite.. right ??

[R] ICA for less common data

2007-04-17 Thread Monica Pisica
Hi everyone, I am not sure this is the appropriate list I should put this question to, but I hope you will re-direct me to the most appropriate one if necessary. I am doing an independent component analysis on a dataset that represents different metrics for patchreefs such as depth, area,

[R] Is this a bug?

2007-04-17 Thread Luca Braglia
I have found a strange ifelse behaviour (I think) This works: ifelse(T,1+1,1+2) [1] 2 ifelse(F,1+1,1+2) [1] 3 Maybe I missed something about R internals, but why ifelse(T,print(hello),print(goodbye)) [1] hello [1] hello ifelse(F,print(hello),print(goodbye)) [1] goodbye [1] goodbye values

Re: [R] Is this a bug?

2007-04-17 Thread Weiwei Shi
one is returned value, the other one is the result from print t0 - ifelse(T, print(h), print(e)) [1] h t0 [1] h HTH, weiwei On 4/17/07, Luca Braglia [EMAIL PROTECTED] wrote: I have found a strange ifelse behaviour (I think) This works: ifelse(T,1+1,1+2) [1] 2 ifelse(F,1+1,1+2) [1]

Re: [R] Is this a bug?

2007-04-17 Thread Roland Rau
On 4/17/07, Luca Braglia [EMAIL PROTECTED] wrote: I have found a strange ifelse behaviour (I think) Don't you think it is rather consistent behavior? ifelse(T,1+1,1+2) [1] 2 ifelse(F,1+1,1+2) [1] 3 ifelse(T,hello,goodbye) [1] hello ifelse(F,hello,goodbye) [1] goodbye

Re: [R] Greatest common divisor of two numbers

2007-04-17 Thread AbouEl-Makarim Aboueissa
Thank you so much All. It works. Abou == AbouEl-Makarim Aboueissa, Ph.D. Assistant Professor of Statistics Department of Mathematics Statistics University of Southern Maine 96 Falmouth Street P.O. Box 9300 Portland, ME 04104-9300 Tel: (207) 228-8389 Email: [EMAIL

Re: [R] help comparing two median with R

2007-04-17 Thread Peter Dalgaard
Thomas Lumley wrote: There actually is an exact test for the median that does not assume a location shift: dichotomize your data at the pooled median to get a 2x2 table of above/below median by group, and do Fisher's exact test on the table. Fascinating. But can one be sure that the

Re: [R] Is this a bug?

2007-04-17 Thread Luca Braglia
On 17/04/07 - 14:59, Roland Rau wrote: On 4/17/07, Luca Braglia [EMAIL PROTECTED] wrote: ifelse(T,1+1,1+2) [1] 2 ifelse(F,1+1,1+2) [1] 3 ifelse(T,hello,goodbye) [1] hello ifelse(F,hello,goodbye) [1] goodbye ifelse(T,print(hello),print(goodbye)) [1] hello [1] hello

Re: [R] Tcltk

2007-04-17 Thread Prof Brian Ripley
I suspect tcl's own version of 'access', but can you please confirm that this still happens under 'Run as Administrator', assuming 'C:\Program' is a system area in Swedish Windows Vista? I will be able to take a closer look, but not before 2.5.0 (which is in code freeze and I have limited

[R] Question of the impact of the pilot experiment on the overal statistic interpretation of the subsequent work

2007-04-17 Thread Bruce Ling
Hi, I have a question regarding the impact of the pilot experiment on the overall statistic interpretation of the subsequent work. The context is as following: In a lab there are one Professor and three graduate students A, B, and C. They are working on analysis of some disease to discover

Re: [R] Dealing with data frame column names beginning with a numeric

2007-04-17 Thread Patrick Connolly
On Mon, 16-Apr-2007 at 10:12PM +0100, Prof Brian Ripley wrote: | On Mon, 16 Apr 2007, Duncan Murdoch wrote: | The name change happens in the conversion to a dataframe, so why not | change the name afterwards? That is: | | df - data.frame(mat) | names(df)[3] - 5T | boxplot(df, main=blah blah

Re: [R] extracting intercept from ppr fit

2007-04-17 Thread Vadim Ogranovich
Sorry for triple-posting : I seem to have a problem w/ my mail client. Hi, Is there a way, documented or not, to extract the intercept term (the alpha_0 the MASS book) from a ppr() (Projection Persuit Regression) fit? Thanks, Vadim ## Example: n - 1000 data - data.frame(x= rnorm (n),

[R] Runing R in a bash script

2007-04-17 Thread Ulrik Stervbo
Hello! I am having issues trying to plot to a ong (or jpg) when the R-code in a bash script is executed from cron. I can generate a pdf file, but when I try to write to a png, the file is created, but nothing is written. If I execute the bash script from my console, everything works file. Any

Re: [R] Question of the impact of the pilot experiment on the overal statistic interpretation of the subsequent work

2007-04-17 Thread Charles C. Berry
Bruce, Far below you ask for comment on an issue of statistical interpretation of a collection of biological experiments. Your university (judging from your email handle) has one of the best statistics departments in the world and some of the best biostatisticians in the galaxy. You would

Re: [R] Runing R in a bash script

2007-04-17 Thread Jeffrey Horner
Ulrik Stervbo wrote: Hello! I am having issues trying to plot to a ong (or jpg) when the R-code in a bash script is executed from cron. I can generate a pdf file, but when I try to write to a png, the file is created, but nothing is written. If I execute the bash script from my

Re: [R] Extracting approximate Wald test (Chisq) from coxph(..frailty)

2007-04-17 Thread Mohammad Ehsanul Karim
Dear list, I need to extract the approximate Wald test (Chisq) so that I can put it in a loop. str seemed like a great idea, but I cannot seem to find the approximate Wald test for frailty (in the example data below: 17.89 and its p-value 0.12000) there. I cannot seem to find it in capture.output

Re: [R] Extracting approximate Wald test (Chisq) fromcoxph(..frailty)

2007-04-17 Thread Steven McKinney
Hi Mohammed, Here's one way to get the information you need. First I examined the output of your coxph() call: library(survival) kfitm1-coxph(formula = Surv(time, status) ~ age + + sex +disease + frailty(id, dist = gauss), + data = kidney) class(kfitm1) [1] coxph.penal coxph

Re: [R] help comparing two median with R

2007-04-17 Thread Thomas Lumley
On Tue, 17 Apr 2007, Frank E Harrell Jr wrote: The points that Thomas and Brian have made are certainly correct, if one is truly interested in testing for differences in medians or means. But the Wilcoxon test provides a valid test of x y more generally. The test is consonant with the

Re: [R] Runing R in a bash script

2007-04-17 Thread Henrik Bengtsson
Or see png2() in R.utils, which imitates png() but uses bitmap(), which in turn uses postscript-to-png via ghostscript. BTW, personally I think PNGs generated via bitmap() look way better than the ones generated via png(). /Henrik On 4/17/07, Jeffrey Horner [EMAIL PROTECTED] wrote: Ulrik

[R] Data Manipulation using R

2007-04-17 Thread Anup Nandialath
Dear Friends, I have data set with around 220,000 rows and 17 columns. One of the columns is an id variable which is grouped from 1000 through 9000. I need to perform the following operations. 1) Remove all the observations with id's between 6000 and 6999 I tried using this method. remdat1

[R] Einladung - Invito

2007-04-17 Thread Wirtschaftskammer Schweiz-Bulgarien
(Versione italiana più sotto) Sehr geehrte Damen und Herren Wir haben bemerkt, dass viele Firmenemailadressen keine HTML-Emails empfangen dürfen. Aus diesem Grund senden wir Ihnen unsere Einladung auch in Textformat zu, damit Sie diese rechtzeitig erhalten, denn die Anmeldefrist läuft bald ab.

Re: [R] how to estimate dose from respond given drc package result

2007-04-17 Thread Gabor Grothendieck
Try this: library(drc) # get dose corresponding to response (y = 2) model1 - multdrc(rootl ~ conc, data = ryegrass) f - function(x, y) y - predict(model1, data.frame(conc = x))[1] uniroot(f, c(0, 10), y = 2) On 4/17/07, John McNeil [EMAIL PROTECTED] wrote: Dear all, I can use the very nice

Re: [R] Data Manipulation using R

2007-04-17 Thread Charilaos Skiadas
On Apr 17, 2007, at 8:03 PM, Anup Nandialath wrote: Dear Friends, I have data set with around 220,000 rows and 17 columns. One of the columns is an id variable which is grouped from 1000 through 9000. I need to perform the following operations. 1) Remove all the observations with id's

Re: [R] how to estimate dose from respond given drc package result

2007-04-17 Thread John McNeil
Wow, thanks, that does it perfectly John On Apr 17, 2007, at 5:51 PM, Gabor Grothendieck wrote: Try this: library(drc) # get dose corresponding to response (y = 2) model1 - multdrc(rootl ~ conc, data = ryegrass) f - function(x, y) y - predict(model1, data.frame(conc = x))[1] uniroot(f,

[R] R-OSX error: 'memory not mapped'

2007-04-17 Thread Atte Tenkanen
Hi, I often get a following error with R *** caught segfault *** address 0x78807e00, cause 'memory not mapped' Possible actions: 1: abort (with core dump) 2: normal R exit 3: exit R without saving workspace 4: exit R saving workspace Selection: The system is OSX 4.9 and R-version 2.4.1. Is

Re: [R] Runing R in a bash script

2007-04-17 Thread Ulrik Stervbo
As I har problems installing the Cairo package, I went for Henriks solution - and it works almost perfect. I would like to have been able to generate transparent png. Thanks for the help Ulrik On 18/04/07, Henrik Bengtsson [EMAIL PROTECTED] wrote: Or see png2() in R.utils, which imitates png()