alanarnholt / bsda Goto Github PK

Hello,

I am using the z.test function from the BSDA package in an attempt to analyze some data and I keep getting an output of "NA". There are no errors when I run the test, but all of my outputs ( z=, p-value =, etc) simply say NA. Below is the code I have written to run the analysis, hopefully you can provide some insight.

Dat1 <- Dat[Dat$Year == 1,]
Dat2 <- Dat[Dat$Year == 2,]

sd16 <- sd(Dat1$ColumnA, na.rm = T)
sd17 <- sd(Dat2$ColumnA, na.rm = T)

z.test(x = Dat1$ColumnA, sigma.x = sd16, y = Dat2$ColumnA, sigma.y = sd17)

Hello, It seems that z.test function from BSDA package does not work when we have NA in the sample, because of the value of sigma.