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

看板R_Language作者allen1985 (我要低調拯救形象)時間10年前 (2015/11/22 10:04)推噓2(2推 0噓 21→)

留言23則, 3人參與討論串2/2 (看更多)

首先這標題可能有誤我的理解這問題跟迴歸沒有關係如果我的理解沒有錯題目是: 現在有一組資料是10條繩子的長度那麼想了解繩子的長度應該是什麼統計分配呢是不是常態分配 or 指數分配 or Poisson分配我猜原po想做的事情是隨機產生10組數據從上述的分配然後看看是否跟繩子長度的分配一樣如果到這邊我的理解還是對的原po可以往下看 1. 請先了解生成隨機變數的指令請google一下例如 rnorm(10,10,90) 跟你的敘述差距很大很大只有第一個10是對的產生10筆資料還有 dunif 跟 runif 的差別等等總之三個指令都誤用了 2. 回到問題如果沒有學過很多檢定方式例如 Kolmogorov–Smirnov test 建議的做法是畫直方圖如果某一組隨機生成的數據長相跟繩子長度類似那麼我想你就接近答案了 3. 跟第一點相關產生隨機變數的指令有些參數需要輸入一個簡單的做法利用資料去算出這些參數例如 sample mean and sample variance 其實在算出這些值之後也可以利用這些值刪除某些不太可能的分配舉例來說 Poisson distribution 的特性是什麼 mean = variance? 4. 最後想要做任何分析之前請把統計課本再拿出來翻一下瞭解一下其意義再來實際做分析就我看來 T-test 迴歸在這裡毫無用處不知道為什麼你會提到 ===========================分隔線======================== Sorry 這篇回文跟R沒什麼關係小弟我實在不敢在這版上回R相關指令問題版上神人太多了如果版主覺得不適合可以刪文謝謝 ※ 引述《eegame (yu)》之銘言： : 初學者（兩個月） 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檢定? 還是用數值回歸??? : 這兩個月沒有看到類似的文章 : 麻煩各位前輩指點迷津 : 拜託拜託 -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 76.181.80.117 ※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1448157872.A.10B.html ※ 編輯: allen1985 (76.181.80.117), 11/22/2015 10:07:34

推

celestialgod

11/22 13:44, , 1^F

11/22 13:44, 1^F

推

eegame

11/22 14:06, , 2^F

11/22 14:06, 2^F

→

eegame

11/22 14:06, , 3^F

11/22 14:06, 3^F

Hello 看到你的推文我想你完全沒有了解我的意思今天你想要看你的資料是趨近於什麼分配你也學了不少分配每個分配都有其特性簡單講你就是要想辦法把資料的特性找出來以及了解你想比較之分配的特性看看有沒有相似之處舉例來說常態分配是鐘形曲線你的資料有嗎如果你的資料是雙峰的那麼基本上就不太可能是常態而你提到一堆檢定下次你打出某個檢定的名字請你在後面打出這個檢定是做什麼用的例如 T-test 是檢定兩組連續行資料的平均數是否相同不然大家很難幫你因為很難了解你到底要做什麼

→

celestialgod

11/22 14:10, , 4^F

11/22 14:10, 4^F

→

celestialgod

11/22 14:10, , 5^F

11/22 14:10, 5^F

→

celestialgod

11/22 14:11, , 6^F

11/22 14:11, 6^F