similar to: What's the Edgeworth expansion of a random vector

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

Dear R users, I want to know, given a curve f in d-dimensional space, It is possible to keep the curve’s shape and size unchanged by an arbitrary dxd orthnormal matrix A? That is, the new curve g = A*f is still the same shape and size as f? Thanks for your advices and answers. Fred [[alternative HTML version deleted]]

Hi Everyone, Following the previous discussion on optimizing *wilcox functions, Andreas Loeffler brought to my attention a few other bugs in `wilcox` family functions. It seems like these issues have been discussed online in the past few months, but I haven?t seen discussion on R-devel...unless I missed an email, it seems like discussion never made it to the mailing list. I haven?t seen any bug

I have a dataframe with many hundreds of survey records containing: tripid_nu, lineon, MuchOtherData, X1, X2, X3, X4 tripid_nu = identifier for each "Trip" {E.g. for EB MTA-901 leaving the first stop at 07:24} lineon = identifier for each of the stop locations Where X1 is the ratio (for each trip/stop combination) of some target number to the number of records.

Dear All Maybe the following is a stupid question. Assume I have 3 coordinate points (not limited to be in 2D or 3D space) a, b, c. It is known that these 3 points will define a plane. The problem is how to get the normal direction that is orthogonal to this plane. Is there an easy way to calculate it using the values of a, b, and c? Thanks for any point or help on this. 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

I have the following array: 3 dimensional object, one dimension being year. Object is 3*3*3 library(plyr, reshape) a1<-rep(c(2007,2008,2009),9) a2<-c(rep("a",9),rep("b",9),rep("c",9)) a3<-c(rep(c(rep(1,3),rep(2,3),rep(3,3)),3)) a4 <- rnorm(27) A<-data.frame(cbind(comp=a2,val=a3, year=a1, a4)) A1<-melt(A, id=c("year", "comp",

I believe eigen(), svd() and qr() can all do it. Andy > From: Jonathan Baron > > On 06/17/04 19:04, Fred wrote: > >Dear R-listers: > > > >I am trying to test an algorithm on a set of linearly > independent vectors > >{x1,x2,...,xn}. > > Well, here's an idea, for 10 vectors of length 10, > as columns of a matrix m1. The 11th seems to be needed.

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

> Then the question would be if this test could be replaced with a new > argument to sample, e.g. expandSingle, which has TRUE as default for > backward compatibility, but FALSE if you dont want population to be > expanded to 1:population. It could certainly be useful in some cases, > but you still need to know about the expansion to use it. I think most > of these bugs

Dear Statistics-Experts, Assume you have given a new and untested pseudo-random number generator (prng) and you want to test if it "works". The distribution function (cdf) from which the prng is supposed to sample is known. Further, you are given some finite (large) sample from the prng. If the sample is one-dimensional, we can apply the cdf to it and test the sample for being

Dear R-users, I'd like to announce the release of the new package JM (JM_0.1-0 available from CRAN) for the joint modelling of longitudinal and time-to-event data. The package has a single model-fitting function called jointModel(), which accepts as main arguments a linear mixed effects object fit returned by function lme() of package nlme, and a survival object fit returned by either

Dear R-users, I'd like to announce the release of the new package JM (JM_0.1-0 available from CRAN) for the joint modelling of longitudinal and time-to-event data. The package has a single model-fitting function called jointModel(), which accepts as main arguments a linear mixed effects object fit returned by function lme() of package nlme, and a survival object fit returned by either

tree: https://kernel.googlesource.com/pub/scm/linux/kernel/git/mst/vhost.git linux-next head: f2c4b499aecc0c5a1ec67f3a2a7310cb7168a8ab commit: 4c145987a955269da79312a79ec26638712644bb [8/9] vhost: block speculation of translated descriptors config: mips-malta_kvm_defconfig (attached as .config) compiler: mipsel-linux-gcc (GCC) 7.4.0 reproduce: wget

Hi, the section "Global Parameters" of the current http://rsync.samba.org/ftp/rsync/rsyncd.conf.html says: "You may use references to environment variables in the values of parameters. String parameters will have %VAR% references expanded as late as possible (when the string is used in the program), allowing for the use of variables that rsync sets at connection time, such as

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:

First the R question. I have the results of a rather large survey (thousands of forms, each with dozens of questions) with some existing weights and expansion factors. I wish to add additional weighting factors, based on new information that elements of certain variables should appear in certain proportions. Where should I look in R to develop what is essentially an N-dimensional balancing

Hi R users, I am using the following model to analyze data from a factorial experiment (randomized complete block design with no replication within blocks): model <- glm(survival ~ density * vegetation + (1|block), data=sal2005, family="binomial") Does R use a binomial distribution in this formulation to model random effects or a normal distribution (in which case the analysis is

Dear all, I have a problem formulating random effects when using the nlme package. ?nlme says that the random effects can be a formula or a list of formulae. But when trying to do so, I run into trouble. (R-1.2.2, nlme-3.1-10, Linux) library(nlme) example(nlme) # works nice ... fm1 # the fitted model ### Trying to update with different random effects

Currently on PowerPC, calls to memcmp are not expanded and are left as library calls. In certain conditions, expansion can improve performance rather than calling the library function as done for functions like memcpy, memmove, etc. This patch (https://reviews.llvm.org/D28163) is an initial implementation for PowerPC to expand memcmp when the size is an 8 byte multiple. The approach currently