similar to: expressing functions

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 All, I am using the allShortestPath function based on Floyd's algorithm in e1071 package. It runs great when I have less than 5000 nodes. But when I tried to work on more than 5000 nodes, I ran into memory problem. The problem I really want to solve has 10000-15000 nodes. Does anybody know how to deal with this problem? Are there any other packages in R that can handle this problem?

Hi, I couldn't figure out how to use the functions from the ctest library. I'm using the r-base package that comes with debian potato. library("ctest") told me that no such package existed. I checked the CRAN, but no such package was availiable, instead I was told that it would be part of the standard installation. But functions from ctest like shapiro-wilk don't work. The

From: aguitatierra@hotmail.com Sent: Friday, April 05, 2013 11:47 PM To: r-help@r-project.org ; R Help Subject: Reversing data transformation Hi everybody, I would be very grateful if you could give me your thoughts on the following issue. I need to perform Box-Cox (bcPower€) transformation on my data. To do this, I calculated lambda using the function '€powerTransform'€.

Hi everybody, I would be very grateful if you could give me your thoughts on the following issue. I need to perform Box-Cox (bcPower€) transformation on my data. To do this, I calculated lambda using the function '€powerTransform'€. powerTransform(data) However, I got an error message when performing this function: Convergence failure: return code = 52 I was told by John Fox

Hello I have used recently the kolmogorov smirnov test, which is a test of normality. This test is named ks.test() in ctest library of R. I wonder if the results of ks.test () are true, because the results are strange, time to time. thank you for help meriema -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hi! I'm doing Monte Carlo analyses of the distribution of the t-statistics of the parameters of models evaluated with the lm( ) function. Is there an easy way to recover the t-statistics (similarly to using coef to recover the coefficients)? Thanks, joseph --------------------------------- [[alternative HTML version deleted]]

I want to overlay 50 denisty plots on a single plot. For each plot there are 10,000 data points and i want the empirical density of the data. I have not been able to find an easy way to achieve this (I have scoured the manula and website so sorry if I missed it!), does anyone have any suggestions? Thankyou.

I simulate a uniform data and then transformed into Burr(1,3,1) data, which is of pdf: f(x)=[3*(x^2)] / [(1+x^3)^2], x>0 How can I perform a goodness-of-fit test (k-s, anderson-darling,chisq,cramer-von mises,...) on it (should highly accept) to get test-statistics & p-values? Thanks! Sincerely, Shelton Jin -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-

I'm pretty new to R and are trying to do some reliable normality testing... but, can't find the Shapiro Wilk test in R Does some experienced user have such a function that will be wanting to share with me? Or there is maybe some other way to get hte Shapiro Wilk test done... I'll appreciate any hint on this, Thanks -- *********************** Horacio Samaniego Dep. Ecologia P.

Hi, Does anyone know if the D'agostino test is available with R ? Alex

Full_Name: Wei Xu Version: 1.5.1 OS: WindowsME Submission from: (NULL) (63.215.238.92) The P-value for a two.sided test is not consistent with the confidence interval. For example, P-value=0.1372, but the 95% CI doesn't include the H0 value(1). > x [1] 0.80 0.83 1.89 1.04 1.45 1.38 1.91 1.64 0.73 1.46 > y [1] 1.15 0.88 0.90 0.74 1.21 >

Hello, I am trying to call a FORTRAN subroutine from R. The Fortran code is @: http://lib.stat.cmu.edu/apstat/206 It performs a bivariate isotonic regression on a rectangular grid (m X n) matrix. I used the g77 compiler and successfully created a dll file and it also loads successfully from R. But somehow the programs fails to run properly. (I do get the correct result when I compile the

> library(Kendall) > Kendall(1:3,1:3) WARNING: Error exit, tauk2. IFAULT = 12 <<<<<< tau = 1, 2-sided pvalue =1 I believe Kendall tau is well-defined for this case and the reported value is correct; isn't it a bug to give a warning? (And if, e.g., the pvalue is not well-defined in this case, wouldn't it be better to return NA or NaN or something?) Also,

Dear McLeod, I am a PhD student working with the Kendall package in R. I used this package to make the so-called TauKr matrix correlation analysis (Hemelrijk, 1990). However, lately I've been getting this warning: Error exit, tauk2. IFAULT = 12. From the forums I've been looking into, I noticed that when this warning pops up, it means that the p-value of the correlation may be wrong

Hello everybody, I observed some strange behaviour of R (0.64.1/Linux 2.0.25) when trying to use FORTRAN-libraries. I first downloaded "ALGORITHM AS 304" from StatLib. All I changed at the code, was to substitute "REAL"-routines by "DOUBLE PRECISION". Compiling went o.k.: > g77 -fpic -O2 -shared -o rtest.o fisher.f ("g77 --version" gives 2.7.2.1)

Actually doing both sounds good to me. Can you explain further about ringing them all at once? Here is how I tried to make mine work and failed... {global} PHONES0=SIP/2000 PHONES1=SIP/2001 [local] exten => 6001,1,Dial(${PHONES0&PHONES1),20,trf) When I dial 6001 I see my debugger tell me that I am using the wrong syntax. Do you know the correct syntax for ringing them all at once? I

I attempted to install an X100P card but it was not correctly recognized by my Redhat 9 install. I had a test install running without any cards which was working great minus the outward dialing since no cards existed. Now that I have a card, I want to add it to the system. Do I have to scratch the whole current install in order to get the X100P running on my system or is there a way to get it

Does anyone have an update patch file to get Spandsp installed? I'm running asterisk CVS-HEAD-11/19/04-21:53:37 on redhat 9.0 I installed spandsp-0.0.2 when runnig the patch I get patching file Makefile Hunk #1 FAILED at 41. Hunk #2 FAILED at 69. 2 out of 2 hunks FAILED -- saving rejects to file Makefile.rej -------------- next part -------------- An HTML attachment was scrubbed...

Hi. I am trying to modify kmeans function. It seems that is failing something obvious with the workspace. I am a newbie and here is my code: myk = function (x, centers, iter.max = 10, nstart = 1, algorithm = c("Hartigan-Wong", + "Lloyd", "Forgy", "MacQueen")) + { + do_one <- function(nmeth) { + Z <- switch(nmeth, { + Z