use 1-pchisq instead of lower.tail=FALSE to avoid warnings when ncp > 80