# Usage: HT_Data_For_Mue_No_Sigma(x,mue0,H1,alpha) HT_Data_For_Mue_No_Sigma = function(x,mue0,H1,alpha) { library(TeachingDemos) library(car) qqPlot(x) n = length(x) df = n-1 if (H1 == "less") tcritical = qt(alpha,df) else if (H1 == "greater") tcritical = qt(alpha,df,lower.tail = FALSE) else { t1 = qt(alpha/2,df) t2 = qt(alpha/2,df,lower.tail = FALSE) tcritical = c(t1,t2) } print('Critical value:') print(round(tcritical,3)) t.test(x, alternative = H1, mu = mue0) }