使用R做ARL的蒙地卡羅模擬
這是上學期上品質管理的時侯 老師出的作業之一 用任何一種程式語言模擬製程偏移之後 ARL的表現 其實我本來也不會寫 後來在網路上查到有人用C++寫好了 但我不好意思照抄 所以就把C++轉成R 因為品管的老師很喜歡出這個題目 我就放上來解救大家 但是每一個老師的假設都不一樣，所以如果有人要抄的話還是要注意一下 這是假設製程在第100點以後開始偏移，所以前100點是服從N~(0,1)，100點以後就看老師要求的偏移量是什麼，去改k的值 a<-0 # 在陣列中 ' 列 ' 的位移 n<-20 # 樣本數 m1<-0 # 原始 mean sigma<-1 # 原始 sigma L<-3 # 管制界線的寛度 K<-1.5 # 偏移程度 m2<-m1+K*sigma ...