[問題] 想用R寫一個統計回歸分析的程式

作者: eegame (yu)   2015-11-20 14:24:02
初學者(兩個月) R是第一個程式語言
不知道各位前輩有沒做過類似的程式
!!!假設!!!
想用R統計10條繩子的長度
再根據統計特性的率隨機產生10組數據
程式如下
G <- 10 #10組數據
rowname <- c(1:G) #Y軸名稱
colname <- c("c1","c2","c3","c4") #X軸名稱
F <- matrix(,nrow=G,ncol=4,dimnames=list(rowname,colname))?
#開一個10*4的矩陣並以rowname和colname
F[1,"c1"] <- 19
F[2,"c1"] <- 59
F[3,"c1"] <- 48
F[4,"c1"] <- 23
F[5,"c1"] <- 31
F[6,"c1"] <- 48
F[7,"c1"] <- 56
F[8,"c1"] <- 60
F[9,"c1"] <- 74
F[10,"c1"] <- 83
#在矩陣c1的位置放入這10條繩子的調查數據
#假設1符合常態分布
F[,"c2"] <- rnorm(10,10,90)?
#在矩陣c2的位置放入根據常態分佈的機率隨機產生一組新的數據介於10~90之間的10筆數

#假設2符合均勻分布
F[,"c3"] <- dunif(10,10,90)
#在矩陣c3的位置放入根據均勻分佈的機率隨機產生一組新的數據介於10~90之間的10筆數

#假設3符合Poisson distribution分布
F[,"c4"] <- dpois(10, lambda=1)
#在矩陣c4的位置放入根據Poisson distribution分布的機率隨機產生一組新的數據介於1
0~90之間的10筆數據
#結果如下
F
c1 c2 c3 c4
1? 19 -70.255568 0.0125 1.013777e-07
2? 59 9.346476 0.0125 1.013777e-07
3? 48? -2.139539 0.0125 1.013777e-07
4? 23 -33.231206 0.0125 1.013777e-07
5? 31? 72.630937 0.0125 1.013777e-07
6? 48? -6.032507 0.0125 1.013777e-07
7? 56? -8.420308 0.0125 1.013777e-07
8? 60 -69.542557 0.0125 1.013777e-07
9? 74? 12.044775 0.0125 1.013777e-07
10 83 -13.474063 0.0125 1.013777e-07
問題1.我函數的指令是不是給的不對導致c2 c3 c4的結果錯誤?
問題2.關於隨機分布的應用 如何做到根據統計特性產生結果?
問題3.關於數據統計特性如何檢定比較符合 常態分布 均勻分布 負指數分布等等的分布
特性? 是用卡方檢定 T-test檢定? 還是用數值回歸???
這兩個月沒有看到類似的文章
麻煩各位前輩指點迷津
拜託 拜託
作者: celestialgod (天)   2015-11-20 14:28:00
1. 你知道各分配參數意義嗎?2. 你知道第一個問題就知道怎樣根據你觀察的資料去產生適當的隨機數3. 你的目標是檢定什麼?十條繩子長度符合某分配?4. 你學過適合度檢定嗎?另外,我很好奇這跟回歸(標題)有什麼關係啊?5. 你產生隨機數的目的是什麼?6. 會不會你等一下又要刪文啦?
作者: allen1985 (我要低調 拯救形象)   2015-11-21 17:03:00
所以你的問題跟迴歸沒有關係?最合適的分佈特性? 你是說資料最符合哪種分配嗎?你問題中的自變數 應變數 是?建議你簡化一下你的題目 一步一步問 從頭開始大概知道你的問題是什麼了 你有聽過qq-plot嗎對我來說 你用ggplot2 還是什麼function畫圖 不重要重要的是你想要畫什麼圖 你的圖想表達什麼qq-plot 是用來檢測資料是否為常態分配的一種方法而你現在最欠缺的是 你得先把統計方法學好其次才是用R把結果寫出來
作者: celestialgod (天)   2015-11-22 14:25:00
推樓上,重點不是R,他只是工具,你怎麼呈現你的想法才是重點

Links booklink

Contact Us: admin [ a t ] ucptt.com