similar to: Is there an R or S implementation of PAMSIL or PAMMEDSIL

I am managing two Xen containers which are exhibiting a spooky behavior. There are intermittent network failures that affect only the VM''s, and not dom0. This exhibits itself as problems with nfs locks, snmp, and even icmp ping. The arp tables look normal on both the dom0 and the guest. However, it is definitely arp related. Causing arp traffic of any sort results the problem curing

Hello, I was using R (v.2.11.1, 32 bits) and I did the upgrade to R (v.2.12.0, 64 bits). I followed the instructions in R´s FAQ (What´s the best way to upgrade, question 2.8) and updated my packages. However, now, I can´t use the library "vars". When I call it, there is an error message concerning the package "MASS" which couldn´t be updated because it seems to be no more

Dear list, I am looking for an implementation of random rotation matrix generation in R to do a rotation test: I want to use the matrices to create random multivariate normal matrices with common covariance structure and mean based on an observed data matrix. The rRotationMatrix-function in the mixAK-package is an option, but as far as I can tell I need to draw rotation matrices with determinant

Hello, I am looking for some help with this question: how could I test structural breaks in a instrumental variables´s model? For example, I was trying to do something with my model with three time series. tax_ivreg <- ivreg(l_y ~ l_x2 + l_x1+ dl_y | lag(l_x2, -1)+lag(l_x2, -2)+ lag(l_x1, -1)+lag(l_x1, -2)+lag(l_y, -1)+lag(l_y, -2), data=tax1) summary(tax_ivreg) ## after estimating it,

The package `snow' (Simple Network of Workstations), originally announced in early form about a year ago on this list, is now available from CRAN. Snow implements a simple mechanism for using a workstation cluster for ``embarrassingly parallel'' computations in R. The interface, which is based in part on the Python CoW (Cluster of Workstations; http://www.scipy.org/) package, is

Hello, I am trying to replicate the "missing example" of a TSLS estimation in Epple & McCallum (link below) http://wpweb2.tepper.cmu.edu/facultyadmin/upload/ppaper_32774807225408_Epple-McCallum93.pdf According to them, the commands are in: http://www.tepper.cmu.edu/faculty-research/faculty-pages/dennis-epple/simultaneous-equation-econometrics/index.aspx They use the Stata's

[hoping to redeem myself for my last spurious bug report] From ?corRatio: Letting d denote the range and n denote the nugget effect, the correlation between two observations a distance r apart is (r/d)^2/(1+(r/d)^2) when no nugget effect is present and (1-n)*(r/d)^2/(1+(r/d)^2) when a nugget effect is assumed. This disagrees with the C code (corStruct.c) /* Rational class

Hello, I was studying several packages related to time series analysis (urca, vars, tseries). I understand that we can estimate a VECM and also test restrictions on alphas and betas. However, I couldn't find a function that allows me to specify the five cases of VECM (restricted constant, unrestricted constant, restricted and unrestricted trend and no constant). Is there any function that

I have asked this question on Stackoverflow and was told it does not relate to the sites' mission as it is statistical question, thus I brought it here. I am fitting a gam mode in the mgcv package to study associations of environmental pollutants and mortality. The aim is to choose a model with lowest mgcv and also to reduce the PACF to less than < |0.1|. library(gamair) library(mgcv)

I have seen several queries about parameterizing the negative binomial scale parameter. This is called the heterogeneous negative binomial. I have written a function called "nbinomial" which is in the msme package on CRAN. Type ?nbinomial to see the help file. The default model is a negative binomial for which the dispersion parameter is directly related to mu, which is how Stata,

Making data, especially data that have been the subject of published papers, widely available, can be a useful spinoff from the R project, another gift to the scientific community beyond the provision of computing and analytic tools. Nowadays, in a complete publication of a scientific result, there is every reason for the data to be part of that publication. The Gentleman and Lang 2004

I am probably just looking in the wrong place. I am sure there are a number of ways to do this. If anyone could point me in the right direction it would be very much appreciated. Thanks _________________________________________________ Tom Mulholland Senior Policy Officer WA Country Health Service 189 Royal St, East Perth, WA, 6004 Tel: (08) 9222 4062 e-mail:

Yes and no. Given your response it appears that "Proper Case" is not a term that everyone uses. In Excel there is a function "Proper" which in essence changes "this line into something like this" into "This Line Into Something Like This." My look at casefold seesm to be that is is a wrapper of two functions to change text into either Lower or Upper case.So

I am trying to create a new variable which uses the suburb names if HR and HRRES are the same but which uses HRRES if they are different. Any assistance would be appreciated as my brain has just packed up. I'm not sure I can teach myself anymore new tricks this afternoon. HR HRRES SUBURB What I am trying to get 954 Wheatbelt

I am probably missing something quite obvious, but any help would be appreciated. I am continually getting people misreading the lattice plots because they are expecting the strip (with the factor names in them) to be below the graph. Is there anyway of achieving this. Secondly, from a more personal note I find the grid formed by the axes to be a bit overpowering and would like to make it a

Hello R-community: Much of the time I want to use loops to look at graphs, etc. For example, I have 25 plots, for which the names are m.1$medoids, m.2$medoids, ..., m.25$medoids. I want to index the object number (1:25) as below (just to show concept). for (i in 1:25){ plot(m.i$medoids) } I've tried the following, with negative results for ...

I am using some functions from package clusterSim to evaluate the best clusters layout. Here is the features vector I am using to cluater 12 signals: > alpha.vec [1] 0.8540039 0.8558350 0.8006592 0.8066406 0.8322754 0.8991699 0.8212891 [8] 0.8815918 0.9050293 0.9174194 0.8613281 0.8425293 In the following I pasted an excerpt of my program:

Talking about Excel, you can produce excellent graphs in Excel. Yes you have to work at it, but you can get there. The problem is that they are not the default. My gut feeling is that R will make more of an impact in the presentation of graphics than any implementation in Excel. So even if a patent were granted and it made itself into a mainstream package, would it change the world or would the

The question puzzled me at first, because of your use of "library". It looks as if the hdf5 r "package" utilises the "windows hdf5 library" binary. My reading is that you will have to compile the package yourself after you have downloaded the hdf windows dll from hdf.ncsa.uiuc.edu The instructions are in win.readme.txt of the package source which you can download at

Dear all, Can I do a permutation test if the number of individuals in one group is much bigger than in the other group? I searched the literature but I didin´t find any assumption that refers to this subject for permutation tests. Best regards João Fadista Ph.d. student UNIVERSITY OF AARHUS Faculty of Agricultural Sciences Dept. of Genetics and Biotechnology Blichers Allé 20, P.O.