rank - kruskal-wallis contrast test in r -


all,

i perform equivalent of tukeyhsd on rank ordering median shift test such kruskal wallis

x=matrix(c(1,1,1,1,2,2,2,4,4,4,4,4,1,3,6,9,4,6,8,10,1,2,1,3),ncol=2) anova=aov(x[,2]~factor(x[,1])) tukeyhsd(anova) ## tukey multiple comparisons of means ## 95% family-wise confidence level ## ## fit: aov(formula = x[, 2] ~ factor(x[, 1])) ## ## $`factor(x[, 1])` ## diff lwr upr p adj ## 2-1 1.25 -5.927068 8.427068 0.8794664 ## 4-1 -1.35 -7.653691 4.953691 0.8246844 ## 4-2 -2.60 -9.462589 4.262589 0.5617125 kruskal.test(x[,2]~factor(x[,1])) ## ## kruskal-wallis rank sum test ## ## data: x[, 2] factor(x[, 1]) ## kruskal-wallis chi-squared = 1.7325, df = 2, p-value = 0.4205

i analyze contrasts. please help. thanks. rik

if want multiple comparisons after kruskal-wallis test, need kruskalmc function pgirmess package. before can implement function, need transform matrix dataframe. in example:

# convert matrix dataframe dfx <- as.data.frame(x)  # kruskal-wallis test & output kruskal.test(dfx$v2~factor(dfx$v1))       kruskal-wallis rank sum test  data:  dfx$v2 factor(dfx$v1) kruskal-wallis chi-squared = 1.7325, df = 2, p-value = 0.4205  # post-hoc tests & output kruskalmc(v2~factor(v1), data = dfx)  multiple comparison test after kruskal-wallis  p.value: 0.05  comparisons     obs.dif critical.dif difference 1-2    1.75     6.592506      false 1-4    1.65     5.790265      false 2-4    3.40     6.303642      false

Comments

Post a Comment

Popular posts from this blog

jquery - How do you format the date used in the popover widget title of FullCalendar? -

Image
in fullcalendar, how format date used in popover widget title when clicking '+n more'? code: <div class="fc-popover fc-more-popover"> <div class="fc-header fc-widget-header"> <span class="fc-close fc-icon fc-icon-x"></span> <span class="fc-title">tuesday, april 7</span> ... </div> ... </div> you use daypopoverformat that. determines date format of title of popover created eventlimitclick option. must date formatting string. default value "dddd, mmmm d" english , "ll" other languages.
Read more

asp.net mvc - SSO between MVCForum and Umbraco7 -

i need integrate sso between mvcforum (forum of company) , umbraco7 (manages cms company). however struggling shared authentication work. tried several ways suggested on other threads still not running expected. here did: 1. in mvcforum (forum.mywebsite.com) 1.1. machine key in web.xml: <machinekey validationkey="the same key umbraco7" decryptionkey="the same key umbraco7" validation="sha1" decryption="aes" /> 1.2. authenticate forms in web.xml <forms name="my-sso-auth" protection="all" path="/" timeout="43200" domain="forum.mywebsite.com" loginurl="/members/logon" enablecrossappredirects="true" /> 2. in umbraco7 (news.mywebsite.com) 2.1. machine key in web.xml: <machinekey validationkey="the same key mvcforum" decryptionkey="the same key mvcforum" validation="sha1" decryption="aes" /> 2.2. a...
Read more

Python Tkinter keyboard using bind -

i have problem using bind option can click button in keyboard, , call function once key pressed. tried taking other codes have similar purpose , saw i'm doing pretty same, though still have problem must've missed. thankful if me this. here's code: # -*- coding: utf-8 -*- #imports tkinter import * pil import imagetk, image import time import tkmessagebox #===========================================================================================================================================# #tkinter class class app(frame): def __init__(self, master=none): frame.__init__(self, master) self.pack() #===========================================================================================================================================# #game pieces classes class board: def __init__(self, rows = 7, columns = 6, picture_height = none, picture_width = none): self.rows = rows self.columns = columns self.picture ...
Read more