similar to: cross-validation for count data

Dear Germano, Thank you for your fast reply, In the above code, *MYData *is the actual data set. Do not we need to convert *MYData to *the dissimilarity matrix using *pam(as.dist(**MYData**), k = 10, diss = TRUE*)* code line?* *Regards.* On Thu, Aug 17, 2017 at 2:58 PM, Germano Rossi <germano.rossi at gmail.com> wrote: > try this > > MYdata <-

Sorry, I never use pam. In the help, you can see that pam require a dataframe OR a dissimilarity matrix. If diss=FALSE then "euclidean" was use.So, I interpret that a matrix of dissimilarity is generated automatically. Problems may be in your data. Indeed pam(ruspini, 4)$diss write a dissimilaty matrix while pam(MYdata,10)$diss wite NULL 2017-08-17 16:03 GMT+02:00 Sema Atasever

[This email is either empty or too large to be displayed at this time]

try this MYdata <- read.csv2("data.txt",dec='.') library(cluster) cluster.pam = pam(MYdata,10) table(cluster.pam$clustering) filenameclu = paste("clusters", ".txt") write.table(cluster.pam$clustering, file=filenameclu,sep=",") 2017-08-17 10:28 GMT+02:00 Sema Atasever <s.atasever at gmail.com>: > Dear Authorized Sir / Madam, > >

Someone knows the existence of a function to sum the elements of the same place A[i,j] B[i,j] of a matrix?thank you -- Dr. Daniele Riggi, PhD student University of Milano-Bicocca Department of Statistics Building U7, Via Bicocca degli Arcimboldi, 8 20126 Milano, Italy cell. +39 328 3380690 mailto: daniele.riggi@gmail.com [[alternative HTML version deleted]]

Dear R-users, I'm working with the Rpart package and trying to understand how the procedure select the best split in the case the method "class" and the splitting index "Gini" are used. In particular I'd like to have look to the source code that works out the best split for un unordered predictor. Does anyone can suggest me which functions in the sources I should

Hi, I would like to know if it is possible to have a "R code" to estimate the parameters of a mixture of bivariate (or multivariate) normals via EM Algorithm. I tried to write it, but in the estimation of the matrix of variance and covariance, i have some problems. I generate two bidimensional vectors both from different distribution with their own vector means and variance and

Hi R-users,I'm trying to solve a non linear equation, to find the degrees of freedom of a mixture of t student. I'm sure i wrote the minimization equation in the right way, but when i try to run the EM algorithm to estimate the parameters of the mixture, the following error will appear: Error in uniroot(function(z) log(z/2) - digamma(z/2) + 1 - log((z + d)/2) + : f() values at end

Dear Gregory R. Warnes, I'm a phd student in statistics, at the University of Milano Bicocca. I'm interested to the methods of estimate the parameters of the Dirichlet distribution. Do you have implemented an algorithm in R? If so, can you give me the script? or, in general, some helps? Best Regards Gianna --------------------------------------------------------------- Gianna Monti

Dear R-listers, I would like to carry out a very basic descriptive analysis of homicides rates in Italy, taking into account both the spatial dimension (103 provinces) and the temporal dimension (10 years), but no covariates. In practice, what I would like to do is to describe spatio-temporal variation of homicide rates, identifying those combinations of province-year where the homicide rate

Dear all, I'm trying to manage with user defined split function in rpart (file rpart\tests\usersplits.R in http://cran.r-project.org/src/contrib/rpart_3.1-34.tar.gz - see bottom of the email). Suppose to have the following data.frame (note that x's values are already sorted) > D y x 1 7 0.428 2 3 0.876 3 1 1.467 4 6 1.492 5 3 1.703 6 4 2.406 7 8 2.628 8 6 2.879 9 5 3.025 10 3 3.494

Dear R users, I am in trouble using the function solve(matrix) in a loop. When it happens that I have a singular matrix, I get an error message saying so and the loop stops. My question is if there in any way to keep on going in the loop. Thank you in advance, Laura *---------------------------------------------------------* Laura Antolini, PhD Department of

http://bugzilla.netfilter.org/show_bug.cgi?id=706 Jan Engelhardt <jengelh at medozas.de> changed: What |Removed |Added ---------------------------------------------------------------------------- CC| |jengelh at medozas.de --- Comment #1 from Jan Engelhardt <jengelh at medozas.de> 2011-03-03 14:20:30

Great! Please let the mailing list know if you run into any further problems. -- Peter Paolo Pedroni wrote: > > Alle 16:38, venerd=EC 25 agosto 2006, hai scritto: > > Paolo, > > > > I recently fixed some bugs regarding the report descriptor retrieval. > > Please try it with the newest development version from SVN. You can > > follow the instructions at >

hello R-Wizards!!!!! again i'm invoking your presence!! since all the fuss about the paradoxes on a computer algorithm generating "caos" or un-determinancy, i recently grew quite curious about the mechanism underlying the procedure of NUMBER RADOMIZATION. could anyone of you, masters, attach me the algorithm (or source code? is it right?) behind the *sample()* function, so i can

After a crash, in my zpool tree, some dataset report this we i do a ls -la: brwxrwxrwx 2 777 root 0, 0 Oct 18 2009 mail-cts also if i set zfs set mountpoint=legacy dataset and then i mount the dataset to other location before the directory tree was only : dataset - vdisk.raw The file was a backing device of a Xen VM, but i cannot access the directory structure of this dataset. However i

Hi we use samba 3.0.25a on FreeBSD 6.2 and we would like to have in the logs only the successfull and unsuccessfull login attempts. We would like to read in the log file: Authentication for user [%s] -> [%s] -> [%s] succeeded or Authentication for user [%s] -> [%s] FAILED with error We have tried with log level 1 and we get only the successfull logins. We have tried with log level 2

Hello there, i am trying to create a VM using xcp on my system, with os Ubuntu 11.10. I’ve got 3 server on my testbed and in the near future i am going to try a VM live migration from a server to another, but actually I am having few problems.  The first one : i installed the xcp-xapi on 2 of my servers, and when i give the command #xe template-list One of the server gives me the list, while the

Dear R-users, after seeking for help in R-search I did not find any hint on my particular problem. Countless help on "true" overlay, but nothing on this. Please consider the following: par(mfrow = c(2, 1)) T <- seq(0, 20, by = 0.01) ## PLOT 1 plot(T, 30*exp(-0.65*T), type = "l" ) points(T, 30*exp(-0.26 * T), type = "l", lty = "F8") points(T,

Dear R users, today I've got the following problem. Here you are a dataframe as example. There are some SAMPLES for which a CONCentration was recorded through TIME. The time during which the concentration was recorded is not always the same, 10 points for Sample A, 7 points for Sample B and 11 for sample C Also the initial concentration was not the same for the three samples. I would like