I am running SKAT Binary Robust and am running into an issue with some of my datasets.
I will run:
outRobustBurden <-SKATBinary_Robust.SSD.All(
SSD.INFO, obj, method = "Burden")
outRobustSKAT <-SKATBinary_Robust.SSD.All(
SSD.INFO, obj, method = "SKAT")
outRobustSKATO <-SKATBinary_Robust.SSD.All(
SSD.INFO, obj, method = "SKATO")
for the same data set but for "SKAT" and "SKAT-O" I will get an error for most of the genes I am testing:
Error to run SKATBinary for (Insert gene here): Error in if (abs(q - m1)/sqrt(var1) < Cutoff) { :
missing value where TRUE/FALSE needed
If I use the older command, "SKATBinary.SSD.All" instead of the newer Robust function, I do not get these errors.
Thoughts?
Bests,
JFK