similar to: Anyone familiar with Cumulants or good reference books?

Good day everyone, I want to use the Gram-Charlier series expansion to model some data. To do that, I need functions to: 1) Calculate 'n' moments from given data 2) Transform 'n' moments to 'n' central moments, or 3) Transform 'n' moments to 'n' cumulants 4) Calculate a number of Hermite polynomials Are there R-functions to do any of the above?

Hey, R-listers I was recommended to try using smooth.spline function for estimating 2-Dimensinal curve given a data set. So will you please tell me where to get this R function? Or which package provides this function? Thanks for your point. Fred

Hi, R-listers I tried to plot several graphs in a sigle x-y coordinate settings, like the following: |(y) s | ****** s | ***** s | sssssssssssssssssss |_______________________________(x) where "*" and "s" denote two diffrent plots. However, when I used plot(data1); % data1 is the data points of "*"

Hey, I want to draw several plots sequently, but have to make them dispaly in one figure. So how to achieve this? Thanks. Fred

When I try to install the R RPM under Mandrake Linux 7.1 (intel), I get the following error: [root at edgeworth programs]# rpm -ivh R-base-1.1.1-1.i386.rpm error: failed dependencies: libreadline.so.3 is needed by R-base-1.1.1-1 [root at edgeworth programs]# I do have realine installed: [root at edgeworth programs]# rpm -qa | grep readline readline-4.1-4mdk and

Dear R-listers, I have a dxr matrix Z, where d > r. And the product Z*Z' is a singular square matrix. The problem is how to get the left inverse U of this singular matrix Z*Z', such that U*(Z*Z') = I? Is there any to figure it out using matrix decomposition method? Thanks a lot for your help. Fred

I am sorry to ask such a simple question, but my knowledge of programming language is extremely limited, and I couldn't find the answer in any of the normal resources. I want to do a for loop in which m takes on the values 1, 5, 10, 15, 20. What is the syntax for doing that? I had been doing a loop for m in 1:20, but I only want those values above. Thanks, Brian

hi all can anyone tell me what is head() function in r. what is it's utility? thanks in advance --------------------------------- [[alternative HTML version deleted]]

dear R group: I have the probably fairly common problem that I would like to have one code.R file do different things if it is invoked from a symbolic link, which should be easy to uncover. $ ln -s code.R code-0.R $ ln -s code.R code-1.R $ R CMD BATCH code-1.R what needs to be in code-1.R to put code-1.r into a character vector? help appreciated. regards, /ivo welch PS : I read

I have just upgraded from Red Hat Linux 5.2 to version 6.0 (with GNOME) on an intel box. I have also upgraded to the latest R version (0.64.1-2 rpm) and have also installed the corresponding GNOME rpm. However, I am not able to start R: [cribari at edgeworth install]$ R /usr/lib/R/bin/R.binary: error in loading shared libraries: /usr/lib/R/bin/R.binary: undefined symbol: __setfpucw [cribari at

Full_Name: David Simcha Version: 2.10 OS: Windows XP Home Submission from: (NULL) (173.3.208.5) > a <- c(1:10) > b <- c(1:10) > cor.test(a, b, method = "spearman", alternative = "greater", exact = TRUE) Spearman's rank correlation rho data: a and b S = 0, p-value < 2.2e-16 alternative hypothesis: true rho is greater than 0 sample estimates:

Hi All, I am trying to estimate a functional form in terms of cumulants / moments from the results of a MC simulation. My first idea was to use an Edgeworth expansion, but this looks pretty difficult to implement. Before I go ahead and write some code to do this, could anyone tell me if this has already been done, and which package I could find it in if it has? If it hasn't, does anyone know

Do you know of a reference that discusses alternative choices for plotting positions for a normal probability plot? The documentation for qqnorm says it calls ppoints, which returns qnorm((1:m-a)/(m+1-2*a)) with "a" = ifelse(n<=10, 3/8, 1/2)? The help pages for qqnorm and ppoints just refer to Becker, Chambers and Wilks (1988) The New S Language (Wadsworth & Brooks/Cole),

I am trying to install R under Red Hat Linux 5.2, but the following error happens: [root at edgeworth cribari]# rpm -ivh R-base-0.63.2-1.i386.rpm failed dependencies: libncurses.so.3.0 is needed by R-base-0.63.2-1 I must be missing something very simple. I would appreciate any guidance. Thanks. FC. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help

> 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,

Hello, I thougth that there is a function which does the kendall test in R, I writed on the console apropos("kendall") and I didn't found anything can you tell me how could I do to use the kendall test? Thanks. _____________________________________________________________________________ l [[alternative HTML version deleted]]

Dear R users I want to compute Kendall's Tau between two vectors x and y. But x and y may have zeros in the same position(s) and I wrote the following function to be sure to drop out those "double zeros" "cor.kendall" <- function(x,y) { nox <- c() noy <- c() # for (i in 1:length(x)) if (x[i]!= 0 | y[i] != 0) nox[length(nox)+1]<- x[i] for (i in

Hi every body: I need some help with kendall.global. The example in the manual seems not working well, and cannot used with my data, always the same error. data(mite) > mite.hel <- decostand(mite, "hel") > > # Reproduce the results shown in Table 2 of Legendre (2005), a single group > mite.small <- mite.hel[c(4,9,14,22,31,34,45,53,61,69),c(13:15,23)] >

I know this is a dumb question but I'm pulling out what little hair I have left. I download the package and got the following: trying URL ' http://ftp.ussg.iu.edu/CRAN/bin/windows/contrib/2.14/Kendall_2.2.zip' Content type 'application/zip' length 52732 bytes (51 Kb) opened URL downloaded 51 Kb package ‘Kendall’ successfully unpacked and MD5 sums checked The downloaded

Hey, All In principal component analysis (PCA), we want to know how many percentage the first principal component explain the total variances among the data. Assume the data matrix X is zero-meaned, and I used the following procedures: C = covriance(X) %% calculate the covariance matrix; [EVector,EValues]=eig(C) %% L = diag(EValues) %%L is a column vector with eigenvalues as the elements percent