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

Hi all .. there is no doubt a simple answer to this, but it eludes me. In the first session below ( with jarque.bera.test) you will see that p.value prints with a name of X-squared . This is easily fixed by changing the source to assign a more appropriate name - no name is assigned in the source listing below (the original source code of jarque.bera.test() from tseries).. but what I

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,

I may have a deficiency in understanding the procedure for ADS authentication with samba, but most of the server setup works so far. I have bound a Red Hat Enterprise 5 server to our windows domain, it shows up in DNS and ADS, I can ping it, but I can't get samba shares to be accessible to users, or even get the smbclient to return shares properly. wbinfo -g returns the domain groups

I am trying to use the value of an ID variable in an if statement and not getting the results I expected. # ID values for two school districts > with(rf, tapply(DistrictID, DistrictName, min) ) Aberdeen School Dist. # 58 Buhl Joint School District 59340 53409 This creates DNAME as I expected ...

We've just been hit rather badly by a very nasty bug that can cause rsync to silently discard files or fill them with zeroes. It happens when e.g. opendir() succeeds but readdir() returns an error; rsync does not check for an error from readdir() and so simply ignores the error (along with any remaining files in the directory). No error is reported to the user, who will then happily

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)