similar to: R-beta: Questions

Greetings! I want to have the coefficients that R uses in shapiro.test() for the Shapiro-Wilk test for a prticular sample size, i.e. the a[i] in W = Sum(a[i]*x[i])/(Sum(x[i] - mean(x))^2) (where the x[i] are sorted). Two questions: Q1: Is there a readymade R function from which I can extract these? Q2: I was wondering if I might be able to modify the code for the function shapiro.test() so

Hi everyone, ? I am trying to get a point of intersection between a polyline and a straight line ?.. and get the x and y coordinates of this point. For exemplification consider this: ? ? set.seed(123) ? k1 <-rnorm(100, mean=1.77, sd=3.33) ?k1 <- sort(k1) q1 <- rnorm(100, mean=2.37, sd=0.74) q1 <- sort(q1, decreasing = TRUE) plot(k1, q1, xlim <- c((min(k1)-5),

I have two CSV files (exported from Excel), say file1 and file2. The have the same number of rows, and each has several columns, with names on the first line; and some of the columns in file1 are repeated in file2. Using the "foreign" package, I can read these in separately to dataframes say d1 and d2 with > d1<-read.csv("file1") >

The formula attribute of the builtin CO2 dataset seems a bit strange: > formula(CO2) Plant ~ Type + Treatment + conc + uptake What is one supposed to do with that? Certainly its not suitable for input to lm and none of the examples in ?CO2 use the above.

Dear All, I got these results from the example in the function "dbeta": >x <- seq(0, 1, length=21) > dbeta(x, 1, 1) [1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 Any Idea? TIA Giovanni dr. Giovanni Parrinello Section of Medical Statistics Department of Biosciences University of Brescia 25127 Viale Europa, 11 Brescia Italy Tel: +390303717528 Fax: +390303701157

Dear all, I think this question has been raised very recently but I think I didn't see repsonses on the list nor did I find the thread in the archives. (If it is there I'd appreciate directions ...) I am about to install Linux on a Compaq notebook. One important thing I want to do is being able to run R. I have fiddled with different distributions of Linux (LST, SuSE, Debian) and can

I've put up R-0.62.0.tgz up for FTP from Auckland some minutes ago. As usual, don't get it from there unless you are desperate, but wait for it to be mirrored at a CRAN site near you within a day or two. Note that beginning with this version, you need Fortran in some form or other to compile R. (f2c is good enough) Here is the relevant part of the CHANGES file: CHANGES IN R VERSION

I've put up R-0.62.0.tgz up for FTP from Auckland some minutes ago. As usual, don't get it from there unless you are desperate, but wait for it to be mirrored at a CRAN site near you within a day or two. Note that beginning with this version, you need Fortran in some form or other to compile R. (f2c is good enough) Here is the relevant part of the CHANGES file: CHANGES IN R VERSION

I've put up R-0.62.0.tgz up for FTP from Auckland some minutes ago. As usual, don't get it from there unless you are desperate, but wait for it to be mirrored at a CRAN site near you within a day or two. Note that beginning with this version, you need Fortran in some form or other to compile R. (f2c is good enough) Here is the relevant part of the CHANGES file: CHANGES IN R VERSION

Can someone clarify what I cannot make out from the documentation? The function 'clogit' in the 'survival' package is described as performing a "conditional logistic regression". Its return value is stated to be "an object of class clogit which is a wrapper for a coxph object." This suggests that its usefulness is confined to the sort of data which arise in

I noticed that boxplot computes a 95% CI for the median by using median +/- 1.58*IQR./sqrt(n) Where does the 1.58 constant come from? -- Rick White Statistical Consultant U.B.C. Vancouver B.C. Canada rick at stat.ubc.ca -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

Hi Folks, I'm interested in looking (in a beginner's amateurish way) at prosodic/phonetic analysis of recorded speech. In particular I would want to use R to formulate and evaluate specific models. So I would like to ask R people for their recommendations for a program which would a) Take as input a sound file in one of the common formats (".wav", ".au") b)

Hi Folks, I'm seeking advice about how to resume an outer loop following failure of a function within the loop (which issues an error message). Essentially, I'm repeating (simulating) a process which involves random sampling, EM and MCMC. I'm walking on very edge of rather thin ice -- data rather thinly spread over many parameters! Occasionally some component of the loop fails, with

Hi Folks, This is off-topic R-wise, but it may be close to the heart of many R-users, so I think it may be the best place to ask! Users of 'gv' (the "front end" to ghostscript) will be aware of the little window which gives you the x-y coordinates (in points = 1/72 inch) of the position of the "cross-hair" mouse cursor. These coordinates are those of the corresponding

Hi Folks, What is the best way to avoid a function being read in anew (and masking an exiting function) when a definition of it has already been established in R? Reason: Fernando Tusell and I are working up Schafer's 'CAT' for R (basically done now, just needs some cosmetic tidying up). This uses a function 'slice.index', present in S but not in the versions of R we were

Hi Folks, I note that, while the "chisq" functions dchisq(x, df, ncp=0, log = FALSE) pchisq(q, df, ncp=0, lower.tail = TRUE, log.p = FALSE) qchisq(p, df, ncp=0, lower.tail = TRUE, log.p = FALSE) rchisq(n, df, ncp=0) all have a slot for the non-centrality parameter "ncp", of the functions for the t and F distributions: dt(x, df, log = FALSE)

Sorry Folks, I'm sure I could suss out the answer myself but I need it soon ... ! 1. Given a set of 4 variables X,Y,Z,W in a dataframe DF, I make a scatter-plot matrix using splom(DF). 2. I do all regressions of U on V using lm(U~V), where U and V are all 12 different ordered pairs from X,Y,Z,W. 3. Now I would like to superpose the regression lines from (2) onto the corresponding

Greetings Folks, When R code (as entered or read from a courced file) is executed, is it interpreted from the input form every time having once been read in, or do subsequent invocations use an "intermediate" (pre-interpreted) form? Or, putting it another way, is the execution of R code faster second time time round (and later) because the pre-interpretation has already been done once

Hi Folks, I'm currently looking at missing data/imputation methods (including multiple imputation). S-Plus has a "missing data library". What similar resources are available within R? Or does one roll one's own? Best wishes to all, Ted. -------------------------------------------------------------------- E-Mail: (Ted Harding) <Ted.Harding at nessie.mcc.ac.uk>

The following can't be right, first rw2010: > 1 %% 0.001 [1] 0.001 Then rw2001: > 1 %% 0.001 [1] -2.081668e-17 > and the last seems about right. -- Kjetil Halvorsen. Peace is the most effective weapon of mass construction. -- Mahdi Elmandjra -- No virus found in this outgoing message. Checked by AVG Anti-Virus.