Search Archives

Tried to look for package GLLMGibbs-(package for fitting mixed models by Gibbs samppling contributed by Myles' and Clayton) in CRAN but not found it, Does anybody know if it is still available? thanks Oarabile

I'm trying to produce empirical bayes estimates based on the lognormal model in disease mapping Is there a way this can be done in R? thanks Oarabile

I am trying to determine whether two samples are identical or not. I'm aware that somebody can use the Kolmogorov-Smirnoff test to compare empirical distributions, but since my samples have ties I'm not sure if I'm...

I am trying to write a function that will simulate observed counts X and Y as below, I want the function to be able give many replicates, can somebody advise? thanks Oarabile #alpha,n sdx and sdy are constant, and...

I'm trying to write a function that takes a vector of length n and then takes the first value of the vector i.e j=1 and forms a new vector of length n (i.e replicate the first...

Thanks it worked. Oarabile Gabor Csardi wrote: >apply(a, 2, function(x) min(x[x!=0]) ) > >should do it. Might need some improvement if all numbers in a column >can be zero, try it. > >Gabor > >On Tue, May 08, 2007...

Thanks it worked Oarabile Dimitris Rizopoulos wrote: > try this: > > apply(a, 2, function(x) min(x[x > 0])) > > > I hope it helps. > > Best, > Dimitris > > ---- > Dimitris Rizopoulos > Ph.D. Student > Biostatistical Centre > School of Public Health > Catholic University of Leuven > > Address...

Here are replies from Ted and Jasjeet . Thank you both for your help. Oarabile Jasjeet Singh Sekhon wrote: >A bootstrap Kolmogorov-Smirnoff test will have the correct test level >even if there are ties---i.e., even if non-continuous...

Thanks for the help Petr it works, and all who contributed managed with the suggested function below as I needed to remove the zeros and the upper triangular portion of the matrix.: IED <- function(risk) { n <- length(risk) mrisk <- matrix...

I want to do kernel smoothing on the standardised mortality rates at specified location(x-y coordinates) and be able to produce contour map. I've tried to use the package called fields in R (function smooth.2d) but the...