search for: dondrup

Dear all, is there an easy way to sort or get the order of a matrix or data.frame by multiple rows (like ODER BY in SQL), that means if some values in the first column are equal, then sort them (or 'order') by the second etc.? Thanks in advance! Michael -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

...tion(x = numeric()){mean(x)}') # or more evil things > eval(expr) I know I could do something like > system <- function(...) stop ('This is not allowed!') but it's rather likely to miss one of the 'bad' functions. Any ideas would be appreciated. Regards Michael Dondrup

Dear all, I would like to use the Levenberg-Marquardt algorithm for non-linear least-squares regression using function nls.lm. Can anybody help me to find a a way to compute confidence intervals on the fitted parameters as it is possible for nls (using confint.nls, which does not work for nls.lm)? Thank you for your help Michael

I have been using the package vioplot to make boxplot/densityplots. Now I am looking for a variation of this, and I was wondering if someone could give me any tips on how to do what I want. I have nine groups of values, each group containing two sets of numbers which show the difference between two values, such as this: Start dev Stop dev -1 10 5 -2 0

Dear all, the default plotting method for hclust trees looks just fine for few objects like in the example dataset. But when it comes to many objects (some 1000 - I'm trying to visualize clustered microarray data) it renders a tree, that one cannot inspect, because of overlapping text and lines. My question is, is there a way or a plotting parameter for plotting a tree which is wide enough

Hi, Im currently trying to debug a 'error in unprotect: stack imbalance' problem and I am curious about two basic questions on the use of PROTECT and UNPROTECT, which I could not figure out: - which objects have to be protected, namely, if the code is something like: SEXP fun, e; /* get the expression e ... */ fun = eval(e, R_GlobalEnv); /* or like this?: PROTECT(fun = eval(e,