search for: lehmann

Dear r-helpers, I have a small dataset (n<50), and I want to compute the Hodges Lehmann exact confidence interval. So far, I know that "pairwiseCI" has the function "HL.diff". The description is as follows : HL.diff calculates the Hodges-Lehmann confidence interval for the difference of locations by calling wilcox.exact in package exactRankTests ; But when I che...

Does anyone know of an implementation in R of the Hodges-Lehmann nonparametric difference between two groups? I am interested in the estimate of the difference and the CI or significance of that difference. I did some quick searching and didn't see it, but I may not have been looking for the right name, etc. Thanks, Sean

...t linear relationship in our x design-matrix, since all power-values for each subject sum up to 1. How shall we solve this problem: just eliminate one column of x, or introduce a restriction which says exactly that our power data sum up to 1 for each subject? Thanks a lot Christoph -- Christoph Lehmann <lehmann at puk.unibe.ch> University Hospital of Clinical Psychiatry -- Christoph Lehmann <christoph.lehmann at gmx.ch>

...'t like the hard-coded way and would like to use my.group <- array(c("one group", "other group")) but now c(my.group[1] = class.weight[2], my.group[2] = class.weight[1]) gives an error how can I solve this small problem? thanks for a hint christoph -- Christoph Lehmann <christoph.lehmann at gmx.ch>

Dear R-users How can I get the eigenvalues out of an lda analysis? thanks a lot christoph -- Christoph Lehmann <christoph.lehmann at gmx.ch>

...n indexed color, many of its pixels should not cover the basis pixmapGrey pixel, means, for this "in pixmapIndexed not defined pixels" it should be transparent. What would you recommend me to do? Should I go for another solution than pixmap? Many thanks Cheers Christoph -- Christoph Lehmann <lehmann at puk.unibe.ch> University Hospital of Clinical Psychiatry

...ithm does not converge 1: Algorithm did not converge in: (if (is.empty.model(mt)) glm.fit.null else glm .fit)(x = X, y = Y, 2: fitted probabilities numerically 0 or 1 occurred in: (if (is.empty.model(mt)) glm.fit.null else glm.fit)(x = X, y = Y, Many thanks for your help Christoph -- Christoph Lehmann <christoph.lehmann at gmx.ch>

While using the Hodges Lehmann Mean in DescTools (DescTools::HodgesLehmann), I found that it generated incorrect answers (see <https://github.com/AndriSignorell/DescTools/issues/97> https://github.com/AndriSignorell/DescTools/issues/97). The error is driven by the existence of tied values forcing wilcox.test in Base R to s...

...> e <- eigen(t(x) %*% x) better (x must not have a first column of ones): > e <- eigen(cor(x)) > e$val [1] 6.6653e+07 2.0907e+05 1.0536e+05 1.8040e+04 2.4557e+01 2.0151e+00 > sqrt(e$val[1]/e$val) [1] 1.000 17.855 25.153 60.785 1647.478 5751.216 thanks christoph -- Christoph Lehmann Phone: ++41 31 930 93 83 Department of Psychiatric Neurophysiology Mobile: ++41 76 570 28 00 University Hospital of Clinical Psychiatry Fax: ++41 31 930 99 61 Waldau lehmann at puk.unibe.ch CH-3000 Bern 60...

...('Design' library, together with 'Hmisc'). I would like to try therefore a partial proportional odds model Please, could anybody tell me, where to find the code and how to specify such a model ..or any potential alternatives many thanks for your kind help christoph -- Christoph Lehmann <christoph.lehmann at gmx.ch>

...ined treshold <- function(pred) { if (pred < 0.5) pred <- 0 else pred <- 1 return(pred) } and want to use apply it on a vector sapply(mylist[,,3],threshold) but I get: Error in match.fun(FUN) : Object "threshold" not found thanks for help cheers chris -- Christoph Lehmann <christoph.lehmann at gmx.ch>

...the goodness of fit? (2) And: what are the assumptions which must be fulfilled, so that the results of polr() (t-values, etc.) are valid? How can I test these prerequisites most easily: I have a three-level (ordered factor) response and four metric variables. many thanks Christoph -- Christoph Lehmann <christoph.lehmann at gmx.ch>

...ld you recommend to compare classification methods such as LDA, classification trees (rpart), bagging, SVM, etc: 10-fold cv (as in Ripley p. 346f) or leaving-one-out (as e.g. implemented in LDA)? my data-set is not that huge (roughly 200 entries) many thanks for a hint Christoph -- Christoph Lehmann <christoph.lehmann at gmx.ch>

...t;sense" but it is possible, that something like this occurs in my data-set. Since I call gls() for every cell of my array in a loop, I don't want such errors to occur, since this breaks my loop. what is the problem here? What are potential solutions? Many thanks Christoph -- Christoph Lehmann <christoph.lehmann at gmx.ch>

...lem. Quite often it happens, that the graphics are not at the place where I expect them, but (often on a separate page) later on in the pdf. How can I fix this, means how can I define, that I want a graphic exactly here and now in the document? Many thanks and best regards Christoph -- Christoph Lehmann <christoph.lehmann at gmx.ch>

...ync > Before posting, read: http://www.catb.org/~esr/faqs/smart-questions.html > > --------------------------------------- > ------------------------------------------------------------------------ > > Today's Topics: > > 1. error in rsync protocol data stream (Marten Lehmann) > 2. How to exclude directories from source with --relative > (Stuart Anderson) > 3. Re: How to exclude directories from source with --relative > (Wayne Davison) > 4. rsync with proxy squid - authenticated users - doesn't work > (Jean-Jacques ROUGET...

...llow=ulaw allow=gsm allow=alaw ---------------------- in extensions.conf I only added this to lines under section [demo] for testing the calls from gs1 -> gs2 exten => 2209,1,Dial(SIP/2209) exten => 2210,1,Dial(SIP/2210) ------------------------- -- Bye Ernst --------- Ernst Lehmann Email: lehmann@acheron.franken.de

Hello list, i've got some problems getting a E1 line with a E100P up and running (germany). # cat /proc/zaptel/1 Span 1: WCT1/0 "Digium Wildcard E100P E1/PRA Card 0" HDB3/CCS/CRC4 YELLOW RED "YELLOW RED" sounds not so good. When launching asterisk, enabling pri debug on that span, i see outgoing attempts: ;---------- snip -------- > [00 01 7f ] > Unnumbered

...n and then write my own procedure of including/removing variables? or what would be the golden standard for such a case (my "case" is that I have 2 groups (n1=30, n2=15, number of potential variables: 37, no equal variance in the two groups)) many thanks cheers christoph -- Christoph Lehmann <christoph.lehmann at gmx.ch>

The point estimates produced by wilcox.test are perverse (not wrong, just brain damaged). The Hodges-Lehmann estimator that goes with the signed rank test is the median of the Walsh averages. The Hodges-Lehmann estimator that goes with the rank sum test is the median of the pairwise differences. wilcox.test agrees except that it uses the following very peculiar definition of "sample median": i...