[R] p-Value

2007-08-27 Thread amna khan
Hi Sir When we use Kendall Package to obtain Kendall's Tau statistic. Then we also get two-sided p value. What does two-sided p-value mean? The word two-sided is confusing to understand. Kindly provide help in this regard. -- AMINA SHAHZADI Department of Statistics GC University Lahore,

Re: [R] (coxph, se) Obtaining standard errors of coefficients from coxph to store

2007-08-27 Thread joris . dewolf
David, It would be helpful to give an example of what you would like to extract. I guess you know how to extract elements from vectors and lists. However, sometimes the objects returned by functions can be rather complex (output of coxph() is...) A general method to capture printed output is

[R] Problem with save or/and if (I think but maybe not ...)

2007-08-27 Thread Ptit_Bleu
Hi, I recently discovered the R program and I thought it could be useful to me. I have to analyse data saved as .Px file (x between 0 and 8 - .P0 files have 18 lines at the beginning that I have to skip). New files are generated everyday. This is my strategy : In order to analyse the data, I

[R] I again with shorter message and script

2007-08-27 Thread Ptit_Bleu
Hi, I realized that my first message and the script were (maybe) too long and difficult to read. So I tested this shorter one : - setwd(D:/RWork) chemin-d:/Mydata/ relrfichiers-dir(chemin, pattern=.P) rfichiers-paste(chemin,relrfichiers, sep='')

Re: [R] Problem with save or/and if (I think but maybe not ...)

2007-08-27 Thread Prof Brian Ripley
On Mon, 27 Aug 2007, Ptit_Bleu wrote: Hi, I recently discovered the R program and I thought it could be useful to me. I have to analyse data saved as .Px file (x between 0 and 8 - .P0 files have 18 lines at the beginning that I have to skip). New files are generated everyday.

Re: [R] Calculating diameters of cirkels in a picture.

2007-08-27 Thread Bartjoosen
Hi All, I really like to thank you for the answers, while I was searching for some edge detection and clustering algorithms, Moshe came with a simple but effective solution: use the area to find the diameter! But I tried Moshe's solution, but I couldn't figure out what you mean with

Re: [R] Problem with save or/and if (I think but maybe not ...)

2007-08-27 Thread Ptit_Bleu
Dear Prof Ripley, I thank you for your fast answer. In order to follow your advices : I deleted all the objects and the tfichiers.r already created. I changed all the tfichiers.t of the script into tfichiers.rda Then I launched the script twice. The first time, as tfichiers.rda didnt' exist,

[R] Column naming mystery

2007-08-27 Thread Werner Wernersen
Hi, I hope somebody could help me explain what seems mysterious to me? I use this line on a dataframe ae: summaryBy(total_inflated+total~gr1, data=ae, FUN=sum, na.rm=T) and it returns 3 columns as expected and columns gr1 and total_inflated.sumare correct but the total.sum column consists of

Re: [R] Problem with save or/and if (I think but maybe not ...)

2007-08-27 Thread Ptit_Bleu
Dera Prof Ripley, You wrote : What did you intend there? It is not a test of no difference, but a test that each element of the difference is not 0, and furthermore if() expects a test of length one, not the length of nfichiers. I suspect you intended to test length(nfichiers) 0. And of

[R] [SOLVED] save/load - I finally found (to be honnest : jholtman found)

2007-08-27 Thread Ptit_Bleu
The post of jholtman gave me the solution : http://www.nabble.com/problems-saving-and-loading-%28PLMset%29-objects-tf4179541.html#a11885136 Like Quin Wills, I was trying to assign tfichiers.rda to tfichiers. I've just write load(tfichiers.rda) instead of tfichiers-(tfichiers.rda) And now it

Re: [R] Monmonier algorithm

2007-08-27 Thread Thibaut Jombart
Hello, Here is a late answer, but an answer nonetheless to the question I asked almost one year ago on this list: On Wed, 29 Mar 2006, Thibaut Jombart wrote: Hello list, http://tolstoy.newcastle.edu.au/R/help/06/03/24318.html#24322qlink1 / / / does anyone know if Monmonier algorithm is

[R] Confidence intervals for ccf()

2007-08-27 Thread Gustaf Rydevik
Hello, This is not a purely R-question, but perhaps someone can help me anyway. I am trying to estimate the correlation between two time series (which are both basically different types of measurements of the same phenomena), using both cor.test() (with pearson as method) and ccf(). Now,

[R] [R-pkgs] proftools package now available from CRAN

2007-08-27 Thread Luke Tierney
PROFILE OUTPUT PROCESSING TOOLS FOR R = This package provides some simple tools for examining Rprof output and, in particular, extracting and viewing call graph information. Call graph information, including which direct calls

Re: [R] FAQ 7.x when 7 does not exist. Useability question

2007-08-27 Thread John Kane
--- Duncan Murdoch [EMAIL PROTECTED] wrote: Deepayan Sarkar wrote: On 8/23/07, Duncan Murdoch [EMAIL PROTECTED] wrote: On 8/23/2007 11:28 AM, Prof Brian Ripley wrote: On Thu, 23 Aug 2007, John Kane wrote: The FAQ Section 7 is a very useful place for new users

Re: [R] Column naming mystery

2007-08-27 Thread Werner Wernersen
Sorry that the problem description was not sufficient. Here is a self-contained code replicating the problem: require(doBy) x - as.data.frame(matrix(ncol=3,seq(1,12),dimnames=list(c(),c(hh,total,total.inf summaryBy(total+total.inf~hh,x,FUN=sum) What surprises me are the zeros in the

[R] how to write nicely a condition on a loop for (that is, not like I did)

2007-08-27 Thread Ptit_Bleu
Hi again, This is the follow of my post Problem with save or/and if (I think but maybe not ...). In this post, I wrote that I solved my main problem. And it is true. I also wrote that there was still another problem, which I managed to solve. But I think there must be another way to solve it

[R] How can I interpret this test hypothesis test

2007-08-27 Thread Tom La Bone
Since no reply has been posted yet I will give it a shot. runs.test uses the normal approximation and in your case it returned a z score of -1.8732. This z score has a cumulative probability of pnorm(-1.8732,0,1) [1] 0.03052039 If you are concerned about having too many runs and too few

Re: [R] FAQ 7.x when 7 does not exist. Useability question

2007-08-27 Thread Duncan Murdoch
On 8/27/2007 8:52 AM, John Kane wrote: --- Duncan Murdoch [EMAIL PROTECTED] wrote: I like the first, simple suggestion best; I'll put it into R-devel. (With the slight change to use ul.menu instead of just ul, because FAQ 2.7 includes a plain ul list.) Duncan Murdoch Thanks Deepayan

Re: [R] subset using noncontiguous variables by name (not index)

2007-08-27 Thread Muenchen, Robert A (Bob)
Gabor, That works great! I think this would be a very helpful addition to the main R distribution. Perhaps with a single colon representing numerical order (exactly as you have written it) and two colons representing the order of the variables as they appear in the data frame (your first

Re: [R] R 2.5.1 - Rscript through tee

2007-08-27 Thread Dirk Eddelbuettel
On 26 August 2007 at 22:47, François Pinard wrote: | I met a little problem for which someone might have a solution. Let's | say I have an executable file (named pp.R) with this contents: | |#!/usr/bin/Rscript |options(echo=TRUE) |a - 1 |Sys.sleep(3) |a - 2 | | If I

[R] validate (package Design): error message subscript out of bounds

2007-08-27 Thread Wentzel-Larsen, Tore
Dear R users I use Windows XP, R2.5.1 (I have read the posting guide, I have contacted the package maintainer first, it is not homework). In a research project on renal cell carcinoma we want to compute Harrell's c index, with optimism correction, for a multivariate Cox regression and also

Re: [R] subset using noncontiguous variables by name (not index)

2007-08-27 Thread Thomas Lumley
On Mon, 27 Aug 2007, Muenchen, Robert A (Bob) wrote: Gabor, That works great! I think this would be a very helpful addition to the main R distribution. Perhaps with a single colon representing numerical order (exactly as you have written it) and two colons representing the order of the

Re: [R] R-2.5.1 RedHat EL5 compilation failed

2007-08-27 Thread Stefan Grosse
Original Message Subject: [R] R-2.5.1 RedHat EL5 compilation failed From: Wang Chengbin [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Date: 26.08.2007 15:22 I can't get R-2.5.1 compiled under RedHat EL5 with gcc 4.1.1. Configure failed at the following: You don't need to

Re: [R] p-Value

2007-08-27 Thread Daniel Lakeland
On Mon, Aug 27, 2007 at 11:49:19AM +0500, amna khan wrote: Hi Sir When we use Kendall Package to obtain Kendall's Tau statistic. Then we also get two-sided p value. What does two-sided p-value mean? The word two-sided is confusing to understand. Two-sided is sometimes also called

[R] Robust Standar Errors in Zero-Truncated Poisson

2007-08-27 Thread Pedro Mota Veiga
Hi. I would like to know if is it possible to estimate zero-truncated count models with robust standard errors in R. In Stata that is possible. I already made some searches and attempts but not obtained it. In R I made the estimation of the truncated poisson by the vglm command of VGAM package .

[R] Sequential Rank Test

2007-08-27 Thread Bernardo Rangel Tura
Hi R-Masters I need use a sequential approach in serie of cases, but may data is not normal. If data is normal distribution is very easy create analysis using likelihood ratio like of Wald test. But in my case I need use a non-parametric test (Mann-Whitney). I was use:

Re: [R] Sequential Rank Test

2007-08-27 Thread Henrique Dallazuanna
Hi Bernardo, I think that ?wilcox.test will help you. -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40 S 49° 16' 22 O On 27/08/07, Bernardo Rangel Tura [EMAIL PROTECTED] wrote: Hi R-Masters I need use a sequential approach in serie of cases, but may data is not normal. If

Re: [R] Sequential Rank Test

2007-08-27 Thread Birgit Lemcke
I looked for the same topic today and found ?wilcox.test in the stats package. B Am 27.08.2007 um 17:33 schrieb Bernardo Rangel Tura: Hi R-Masters I need use a sequential approach in serie of cases, but may data is not normal. If data is normal distribution is very easy create

[R] FW: subset using noncontiguous variables by name (not index)

2007-08-27 Thread Muenchen, Robert A (Bob)
Thomas, that's a good point. I was thinking of anscombe[x1::y1] making it clear which one, but you would then want just x1::y1 to have unambiguous meaning on its own, which is impossible. As for x1:xN, it's unambiguous on its own. I thought one of the great advantages of R was that it could use

Re: [R] Formatting Sweave in R-News

2007-08-27 Thread Arjun Narayan
Thank you Paul for your response. Unfortunately that did not work. A figure environment frames it neatly, but still contained in only one column. I have tried various methods, but they all seem to not work, or if the solutions involve manually setting the size, the grey column separator still

Re: [R] Formatting Sweave in R-News

2007-08-27 Thread Arjun Narayan
Dear Paul, I stand corrected. Your solution was the right way. The following code now works: (Apparently I still need to specify the width command as my pdf is incorrectly sized by default) \begin{figure*}[b] \begin{center} \includegraphics[width=8in]{generatedPDF.pdf} \end{center} \end{figure*}

[R] Max vs summary inconsistency

2007-08-27 Thread Adam D. I. Kramer
Hello, I'm having the following questionable behavior: summary(m) Min. 1st Qu. MedianMean 3rd Qu.Max. 1 13000 26280 25890 38550 50910 max(m) [1] 50912 typeof(m) [1] integer class(m) [1] integer ...it seems to me like max() and summary(m)[6] ought to return

[R] subset question

2007-08-27 Thread Kirsten Beyer
I would like to code records in a dataset with a 1 if any of the columns 9-67 contain a particular code, and zero if they don't. I've been working with subset and it seems that something like subset(data, data[9:67]--12345) would work, but I have been unsuccessful so far. It seems like a simple

Re: [R] Max vs summary inconsistency

2007-08-27 Thread Thomas Lumley
On Mon, 27 Aug 2007, Adam D. I. Kramer wrote: Hello, I'm having the following questionable behavior: summary(m) Min. 1st Qu. MedianMean 3rd Qu.Max. 1 13000 26280 25890 38550 50910 max(m) [1] 50912 typeof(m) [1] integer class(m) [1] integer ...it seems

Re: [R] subset using noncontiguous variables by name (not index)

2007-08-27 Thread Muenchen, Robert A (Bob)
Thanks for helping me see why R doesn't have the obvious! -Bob -Original Message- From: Thomas Lumley [mailto:[EMAIL PROTECTED] Sent: Monday, August 27, 2007 2:12 PM To: Muenchen, Robert A (Bob) Subject: RE: [R] subset using noncontiguous variables by name (not index) On Mon, 27

Re: [R] Max vs summary inconsistency

2007-08-27 Thread François Pinard
[Adam D. I. Kramer] I'm having the following questionable behavior: summary(m) Min. 1st Qu. MedianMean 3rd Qu.Max. 1 13000 26280 25890 38550 50910 max(m) [1] 50912 ...it seems to me like max() and summary(m)[6] ought to return the same number. Am I doing

Re: [R] Max vs summary inconsistency

2007-08-27 Thread Adam D. I. Kramer
On Mon, 27 Aug 2007, François Pinard wrote: summary(m) Min. 1st Qu. MedianMean 3rd Qu.Max. 1 13000 26280 25890 38550 50910 max(m) [1] 50912 ...it seems to me like max() and summary(m)[6] ought to return the same number. Am I doing something wrong? Some may

Re: [R] how to include bar values in a barplot?

2007-08-27 Thread Frank E Harrell Jr
Donatas G. wrote: On Tuesday 07 August 2007 22:09:52 Donatas G. wrote: How do I include bar values in a barplot (or other R graphics, where this could be applicable)? To make sure I am clear I am attaching a barplot created with OpenOffice.org which has barplot values written on top of each

[R] use apply function with which

2007-08-27 Thread schuurmans
Dear R-users, For a data frame (say in this example X) I want to look up the corresponding value in a 'look-up data frame' (in this example Y). The for-loop works but is very time-consuming because 'X' in reality is very big. Therefore I would like to have a solution with apply. However, I do not

[R] grouping scat1d/rug and plotting to 2 axes

2007-08-27 Thread Mike
Hi, I'm wondering if anybody can offer a bit of guidance on how to add a couple of features to a plot. I'm using Frank Harrell's Design library to model some survival data in R (2.3.1, windows platform). I'm fairly comfortable with the survival modeling in Design, but am still at a

Re: [R] use apply function with which

2007-08-27 Thread Charles C. Berry
On Mon, 27 Aug 2007, [EMAIL PROTECTED] wrote: Dear R-users, For a data frame (say in this example X) I want to look up the corresponding value in a 'look-up data frame' (in this example Y). The for-loop works but is very time-consuming because 'X' in reality is very big. Therefore I would

[R] Rmpi and x86

2007-08-27 Thread Edna Bell
Dear R Gurus: Is there a problem with Rmpi on x86 with SUSE 10.1, please? I've tried everything and it still won't load. Has anyone else dealt with this please? Thanks, Edna Bell __ R-help@stat.math.ethz.ch mailing list

Re: [R] grouping scat1d/rug and plotting to 2 axes

2007-08-27 Thread Frank E Harrell Jr
Mike wrote: Hi, I'm wondering if anybody can offer a bit of guidance on how to add a couple of features to a plot. I'm using Frank Harrell's Design library to model some survival data in R (2.3.1, windows platform). I'm fairly comfortable with the survival modeling in Design, but

[R] oddity with method definition

2007-08-27 Thread Faheem Mitha
Just wondered about this curious behaviour. I'm trying to learn about classes. Basically setMethod works the first time, but does not seem to work the second time. Faheem.

Re: [R] oddity with method definition

2007-08-27 Thread Duncan Murdoch
On 27/08/2007 5:47 PM, Faheem Mitha wrote: Just wondered about this curious behaviour. I'm trying to learn about classes. Basically setMethod works the first time, but does not seem to work the second time. Faheem.

Re: [R] oddity with method definition

2007-08-27 Thread Thomas Lumley
On Mon, 27 Aug 2007, Faheem Mitha wrote: Just wondered about this curious behaviour. I'm trying to learn about classes. Basically setMethod works the first time, but does not seem to work the second time. Faheem.

Re: [R] Calculating diameters of cirkels in a picture.

2007-08-27 Thread Moshe Olshansky
Hi Bart, Let's assume that you situation was simpler - you have a BW (Black and White) image containing circles (in white) and you need to find the diameter of each circle (and of course to know how many circles you have). This can be done with labeling of connected components. You say that two

Re: [R] subset question

2007-08-27 Thread jim holtman
Here is one way of checking to see if a row contains a particular value and setting the contents of a new column: n - 20 # create test data x - data.frame(sample(letters,n),sample(letters,n),sample(letters,n),sample(letters,n)) # add a column indicating if the row contains 'a', 'b' or 'c' x$a -

Re: [R] validate (package Design): error message subscript out of bounds

2007-08-27 Thread Frank E Harrell Jr
Wentzel-Larsen, Tore wrote: Dear R users I use Windows XP, R2.5.1 (I have read the posting guide, I have contacted the package maintainer first, it is not homework). In a research project on renal cell carcinoma we want to compute Harrell's c index, with optimism correction, for a

Re: [R] How to provide argument when opening RGui from an external application

2007-08-27 Thread Sébastien
Thanks everyone. I actually thought about ?Rscript.exe but, having used only Rgui, I thought it was a instruction specific to this interface. I will look into it. Sebastien Gabor Grothendieck a écrit : There are also some batch files that can be used with Rscript on XP and info in the README

[R] Excel

2007-08-27 Thread David Scott
A common process when data is obtained in an Excel spreadsheet is to save the spreadsheet as a .csv file then read it into R. Experienced users might have learned to be wary of dates (as I have) but possibly have not experienced what just happened to me. I thought I might just share it with

[R] Problem with lme using glht for multiple comparisons

2007-08-27 Thread Christian Kost
Hi everyone, I am new to R and have a question that relates to unplanned post-hoc comparisons using the multcomp package after a mixed effects model. I couldn't find the answer to it in the archive or in any manual. I have a dataset in which several plants have been treated in a particular

Re: [R] Excel

2007-08-27 Thread Robert A LaBudde
If you format the column as Text, you won't have this problem. By leaving the cells as General, you leave it up to Excel to guess at the correct interpretation. You will note that the conversion to a date occurs immediately in Excel when you enter the value. There are many formats to enter

Re: [R] Excel

2007-08-27 Thread David Scott
On Tue, 28 Aug 2007, Robert A LaBudde wrote: If you format the column as Text, you won't have this problem. By leaving the cells as General, you leave it up to Excel to guess at the correct interpretation. Not true actually. I had converted the column to Text because I saw the

Re: [R] Excel

2007-08-27 Thread Moshe Olshansky
As far as I understand, changing the format changes the way data is displayed by Excel but this does not change the data itself - if while reading the data Excel decided that it was a date, it is being converted to an integer (the number of days since January 1, 1900 - and they mistakenly think