I am facing a similar error with the waller. test in agricolae. I tried to trace it but I've totally failed to rectify it Here is a glimpse of my data frame: print(obj1) ID code bp stage dose rep day numberinitiated numberhealth numberdead 1 1 BAA n a A 1 0 10 10 0 2 2 BAA n a A 2 0 10 10 0 3 3 BAA n a A 3 0 10 10 0 4 4 BPA n p A 1 0 10 10 0 5 5 BPA n p A 2 0 10 10 0 6 6 BPA n p A 3 0 10 10 0 7 7 BLA n l A 1 0 10 10 0 8 8 BLA n l A 2 0 10 10 0 9 9 BLA n l A 3 0 10 10 0 10 10 BAF b a F 1 0 10 10 0 11 11 BAF b a F 2 0 10 9 1 12 12 BAF b a F 3 0 10 10 0 13 13 BPF b p F 1 0 10 9 1 14 14 BPF b p F 2 0 10 10 0 15 15 BPF b p F 3 0 10 9 1 16 16 BLF b l F 1 0 10 9 1 17 17 BLF b l F 2 0 10 9 1 18 18 BLF b l F 3 0 10 7 3 19 19 BAE b a E 1 0 10 6 4 20 20 BAE b a E 2 0 10 6 4 21 21 BAE b a E 3 0 10 6 4 22 22 BPE b p E 1 0 10 1 9 23 23 BPE b p E 2 0 10 3 7 24 24 BPE b p E 3 0 10 2 8 25 25 BLE b l E 1 0 10 8 2 26 26 BLE b l E 2 0 10 6 4 27 27 BLE b l E 3 0 10 7 3 28 28 BAD b a D 1 0 10 0 10 29 29 BAD b a D 2 0 10 1 9 30 30 BAD b a D 3 0 10 2 8 31 31 BPD b p D 1 0 10 6 4 32 32 BPD b p D 2 0 10 0 10 33 33 BPD b p D 3 0 10 0 10 34 34 BLD b l D 1 0 10 9 1 35 35 BLD b l D 2 0 10 5 5 36 36 BLD b l D 3 0 10 7 3 37 37 BAC b a C 1 0 10 9 1 38 38 BAC b a C 2 0 10 9 1 39 39 BAC b a C 3 0 10 8 2 40 40 BPC b p C 1 0 10 6 4 41 41 BPC b p C 2 0 10 9 1 42 42 BPC b p C 3 0 10 8 2 43 43 BLC b l C 1 0 10 7 3 44 44 BLC b l C 2 0 10 8 2 45 45 BLC b l C 3 0 10 8 2 46 46 BAB b a B 1 0 10 0 10 47 47 BAB b a B 2 0 10 0 10 48 48 BAB b a B 3 0 10 1 9 49 49 BPB b p B 1 0 10 5 5 50 50 BPB b p B 2 0 10 0 10 51 51 BPB b p B 3 0 10 0 10 52 52 BLB b l B 1 0 10 0 10 53 53 BLB b l B 2 0 10 10 10 54 54 BLB b l B 3 0 10 0 10 55 55 MAF m a F 1 0 10 10 0 56 56 MAF m a F 2 0 10 10 0 57 57 MAF m a F 3 0 10 10 0 58 58 MPF m p F 1 0 10 10 0 59 59 MPF m p F 2 0 10 10 0 60 60 MPF m p F 3 0 10 10 0 61 61 MLF m l F 1 0 10 10 0 62 62 MLF m l F 2 0 10 10 0 63 63 MLF m l F 3 0 10 10 0 64 64 MAE m a E 1 0 10 10 0 65 65 MAE m a E 2 0 10 10 0 66 66 MAE m a E 3 0 10 10 0 67 67 MPE m p E 1 0 10 10 0 68 68 MPE m p E 2 0 10 10 0 69 69 MPE m p E 3 0 10 10 0 70 70 MLE m l E 1 0 10 10 0 71 71 MLE m l E 2 0 10 10 0 72 72 MLE m l E 3 0 10 10 0 73 73 MAD m a D 1 0 10 10 0 74 74 MAD m a D 2 0 10 10 0 75 75 MAD m a D 3 0 10 10 0 76 76 MPD m p D 1 0 10 10 0 numberinfected event PctMrt 1 0 0 0 2 0 0 0 3 0 0 0 4 0 0 0 5 0 0 0 6 0 0 0 7 0 0 0 8 0 0 0 9 0 0 0 10 0 0 0 11 0 0 10 12 0 0 0 13 0 0 10 14 0 0 0 15 0 0 10 16 0 0 10 17 0 0 10 18 0 0 30 19 0 0 40 20 0 0 40 21 0 0 40 22 0 0 90 23 0 0 70 24 0 0 80 25 0 0 20 26 0 0 40 27 0 0 30 28 0 0 100 29 0 0 90 30 0 0 80 31 0 0 40 32 0 0 100 33 0 0 100 34 0 0 10 35 0 0 50 36 0 0 30 37 0 0 10 38 0 0 10 39 0 0 20 40 0 0 40 41 0 0 10 42 0 0 20 43 0 0 30 44 0 0 20 45 0 0 20 46 0 0 100 47 0 0 100 48 0 0 90 49 0 0 50 50 0 0 100 51 0 0 100 52 0 0 100 53 0 0 0 54 0 0 100 55 0 0 0 56 0 0 0 57 0 0 0 58 0 0 0 59 0 0 0 60 0 0 0 61 0 0 0 62 0 0 0 63 0 0 0 64 0 0 0 65 0 0 0 66 0 0 0 67 0 0 0 68 0 0 0 69 0 0 0 70 0 0 0 71 0 0 0 72 0 0 0 73 0 0 0 74 0 0 0 75 0 0 0 76 0 0 0

Here is the code: Effmodel <- lm(event ~ dose, data = obj1) summary(Effmodel) anovaEff <- anova(Effmodel) waller.test (Effmodel, "dose", group = TRUE)

Here is the error: waller.test (Effmodel, "dose", group = TRUE) Error in if ((K - IN0/ID0) * (K - IN1/ID1) <= 0) b0 <- t : missing value where TRUE/FALSE needed

Please help me as I#ve been stuck for 4 hours and googling the error is not helping