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

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:

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

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

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

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,

Hello all. I wish to run k-means with "manhattan" distance. Since this is not supported by the function "kmeans", I turned to the "pam" function in the "fpc" package. Yet, when I tried to have the algorithm run with different starting points, I found that pam ignores and keep on starting the algorithm from the same starting-points (medoids). For my

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,

Hello, When using PAM (partitioning around medoids), I would like to skip the build-step and give the fonction my own medoids. Do you know if it is possible, and how ? Thank you very much. Isabel

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

Hello Martin Maechler and All, A simple question (I hope): How can I compute the "sum of the dissimilarities" that appears in the pam command (from the cluster package) ? Is it the "manhattan" distance (such as the one implemented by "dist") ? I am asking since I am running clustering on a dataset. I found 7 medoids with the pam command, and from it I have the

I am using the PAM algorithm in the CLUSTER library. When I allow PAM to seed the medoids using the default __build__ algorithm things work well: > pam(stats.table, metric="euclidean", stand=TRUE, k=5) But I have some clusters from a Hierarchical analysis that I would like to use as seeds for the PAM algorithm. I can't figure what the mediod argument wants. When I put in the