Mann-Whitney U test after zero-inflation correction
0
0
Entering edit mode
2.6 years ago
exin ▴ 60

I'm trying to perform Mann-whitney u test on my data, but realised the excessive zeros cause an issue. I read about the pscl package that uses the zeroinfl function to correct for this, but is unable to find any example online that uses both together. Can I simply perform the test on zeroinfl-corrected data? How do I do that?

Data structure:
time    treat   count
t0  ct  10
t5  ct  1
t10 ct  0
t15 ct  0
t20 ct  0
t25 ct  0
t30 ct  0
t0  ct  7
t5  ct  4
t10 ct  2
t15 ct  2
t20 ct  2
t25 ct  2
t30 ct  2
t0  dm  9
t5  dm  4
t10 dm  0
t15 dm  0
t20 dm  0
t25 dm  0
t30 dm  0
t0  dm  11
t5  dm  10
t10 dm  3
t15 dm  1
t20 dm  1
t25 dm  1
t30 dm  1

wilcox.test(count~treat, data=data, paired=FALSE)
datanew<- zeroinfl(count ~ time+treat, data = data)
R • 712 views
ADD COMMENT
3
Entering edit mode

Let me see if I got it right. There are two replicates from each condition and seven measurement across 30 minutes. IMHO one can't and shouldn't do any statistical test, just plot the four lines and make whatever claim. Any p<0.1 next to it will be worthless or misleading. If there are more replicates, I don't think there is a problem of zero-inflation. Zeros in the data are allowed! It's a ranked test anyway so don't worry about it too much.

ADD REPLY
3
Entering edit mode

The statistical test you're using is rank based, it won't care about excess zeros.

ADD REPLY
0
Entering edit mode

Thanks guys!!

Yeh there're more replicates.

ADD REPLY

Login before adding your answer.

Traffic: 2336 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6