[RkWard-devel] wilcoxon & ansari_bradley tests

SJR stefan_roediger at gmx.de
Sun Feb 4 21:33:52 UTC 2007


Am Sonntag, 4. Februar 2007 21:19 schrieb Thomas Friedrichsmeier:
> On Sunday 04 February 2007 20:54, SJR wrote:
> > Thomas please run the test with less or greater. I get rk.temp.x and so
> > on as output in ".GlobalEnv" and don't know how to fix this.
>
> Can't reproduce with any of the tests, and any type of hypothesis. Does
> the "cleanup" stage get shown in the command log?
>
> Regards
> Thomas

Yes it does

I set for Wilcoxon Test
* first is greater
* conf.int=T
a<-rnorm(100), b<-rnorm(100)


see image in appendix

## Vorbereiten
## Berechne
rk.temp.x <- substitute (a)
rk.temp.y <- substitute (b)
rk.temp <- wilcox.test (eval (rk.temp.x), eval (rk.temp.y), alternative = 
c("greater"), mu = 0.00000000, paired = FALSE, exact = TRUE, correct = TRUE, 
conf.int = TRUE )
## Drucke Ergebnisse
rk.header ("Wilcoxon Test",
	parameters=list ("Comparing", paste (rk.get.description (rk.temp.x, 
is.substitute=TRUE), "against", rk.get.description (rk.temp.y, 
is.substitute=TRUE)),
	"H1", if (rk.temp$alternative == "less")
		paste (rk.get.short.name (rk.temp.y), "is greater than", rk.get.short.name 
(rk.temp.x))
	      else if (rk.temp$alternative == "greater")
		paste (rk.get.short.name (rk.temp.x), "is greater than", rk.get.short.name 
(rk.temp.y))
	      else
		paste (rk.get.short.name (rk.temp.x), "and", rk.get.short.name 
(rk.temp.y), "differ"),"Compute Confidence Interval", "TRUE", "Continuity 
correction in normal approximation for p-value", "TRUE", "Compute exact 
p-value", "TRUE","mu", "0.00000000"))

rk.results (list (
	'Variable Names'=rk.get.description (rk.temp.x, rk.temp.y, 
is.substitute=TRUE),
	'statistic'=rk.temp$statistic,
	'Location Shift'=rk.temp$null.value,
	'Hypothesis'=rk.temp$alternative,
	p=rk.temp$p.value,
	'confidence interval percent'=(100 * attr(rk.temp$conf.int, "conf.level")),
	'confidence interval of difference'=rk.temp$conf.int,
	'Difference in Location' = rk.temp$estimate ))
## Aufräumen
rm (list=grep ("^rk.temp", ls (), value=TRUE))


-------------- next part --------------
A non-text attachment was scrubbed...
Name: graphics device5.png
Type: image/png
Size: 103974 bytes
Desc: not available
URL: <http://mail.kde.org/pipermail/rkward-devel/attachments/20070204/6c443fad/attachment.png>


More information about the Rkward-devel mailing list