[R] individual likelihoods

2004-01-16 Thread Bruno Giordano
Dear all, is there a way to extract individual likelihoods from a glm/lrm object? By individual likelihoods, I mean the likelihoods whose product give the overall likelihood of the model. I guess the code in the base package, used to compute the Akaike Information Criterion may help me. However, I

[R] non-nested models testing

2004-01-14 Thread Bruno Giordano
Hi, I'm trying to figure out some stastistics for non-nested models testing. Among them: [1] Vuong statistic (Vuong, 1989) [2] Distribution-free test by Clarke (2003). Did someone wrote a piece of [R] code to compute any of them? Many thanks, Bruno Bruno L.

[R] non-uniqueness in cluster analysis

2003-12-03 Thread Bruno Giordano
Hi, I'm clustering objects defined by categorical variables with a hierarchical algorithm - average linkage. My distance matrix (general dissimilarity coefficient) includes several distances with exactly the same values. As I see, a standard agglomerative procedure ignores this problems, simply

Re: [R] non-uniqueness in cluster analysis

2003-12-03 Thread Bruno Giordano
What I did was, in presence of equal values distances, to randomize the selection of them, and compute the distortion of the solution using cophenetic correlation. I computed 1 random trees for each of three methods: average, single and complete linkage. Among the randomly selected solutions,

Re: [R] cclust - cindex - binary data

2003-12-01 Thread Bruno Giordano
If anybody had the same problem I had in calculating the cindex with the clustIndex function (binary data), replace this line of code: if (length(unique(x))==2) with this one: if (length(unique(as.vector(x)))==2) located in the call to the procedure for calculating the cindex, in the cclust

[R] cclust - cindex - binary data

2003-11-28 Thread bruno . giordano
Hi, I'm trying to debug a function I wrote to calculate the cindex for a hierarchical tree. For this it is useful to compare my calculations with those in output from the clustindex function, in the cclust library. There's no way, however, to have the cindex value for a given output of the cclust

[R] cclust - cindex - binary data

2003-11-27 Thread Bruno Giordano
Hi, I'm trying to debug a function I wrote to calculate the cindex for a hierarchical tree. For this it is useful to compare my calculations with those in output from the clustindex function, in the cclust library. There's no way, however, to have the cindex value for a given output of the cclust