[問題類型]: 請把以下不需要的部份刪除 程式諮詢(我想用R 做某件事情，但是我不知道要怎麼用R 寫出來) [軟體熟悉度]: 請把以下不需要的部份刪除 入門(寫過其他程式，只是對語法不熟悉) [問題敘述]: 請簡略描述你所要做的事情，或是這個程式的目的 大家好 我想問問大家如果遇到這樣的內容 會怎麼用R寫？ 我在網路上查到aggregate可以針對不同的groups 算出組平均 例如 group X1 A obs1 obs2 X2 B obs1 obs2 用 aggregate就可以把A和B的分組平均算出 我也想用這方法 但我需要的要再多一個步驟 我現在有6248個觀測值 他們來自175個不同的LotID 每個LotID的樣本個數（就是觀測值個數）在20,32,50和80之間 我要根據不同的樣本數隨機抽樣15,20,25,和35 我的做法是 先用table查看每個LotID有多少觀測值 然後建立一個欄位叫做Osize 就可以告訴我每個對應的LotID的本來的size 然後再建立一個新欄位NSize 根據Osize放入我預計要抽樣的數目 然後建立index=sample(Osize, Nsize, replace=FALSE) 這會是每一個LotID我要取的index value 現在就是我要怎麼把個別的LotID分離出來（不一定要放到不同的欄位 但要分開取平均） 我目前的做法是用loop for (i in 1:nrow(Data)) 這讓我從row1跑到最後一個6248 row (if indentical(Data$LotID[i],Data$LotID[i-1])==TRUE) ##這就告訴我第i row 會開始新的LotID 但我困難點在於如何把i這row之前的那i-1得到數值並根據其數量來取樣呢？ 所以我就設定 for(j in 1:175) 因為總共有175個LotID 假設Observation這個欄位是我要的觀察值 A <- as.vector(Data$Observation) A <- as.numeric(A) ##我把觀測值全部抓出來變成一個數列似乎對我個人比較好懂 Data$mean <- mean(Data$Observation[i:i+j*Osize[i]]) 因為如果以i為21列來看 如果前面i-1是20個數值 且都屬於同一個LotID 那我就是要從這20個值當中 取出15個（因為Osize=20 Nsize=15) 但我的困擾在於j是跟著LotID i是跟著觀測值 兩者其實有關係 但我的寫法似乎沒辦法把兩者牽在一起 aggregate求組平均的算法很容易 aggregate(Data$Observation, Data$LotID,summarize, mean) 大這是這樣 以上是我憑記憶打的 不確定 所以如果可以把觀測值根據隨機取樣的樣本先分好 再用aggregate去取平均 是一種做法 但感覺直接根據要娶隨機的size去直接取平均似乎是更簡易的方式 但我目前只會寫loop 不知道大家有沒有什麼建議？ 感謝大家了 [程式範例]: [環境敘述]: 請提供 sessionInfo() 的輸出結果， 裡面含有所有你使用的作業系統、R 的版本和套件版本資訊， 讓版友更容易找出錯誤 [關鍵字]: 選擇性，也許未來有用 -- --

※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 68.180.87.229 (美國) ※ 文章網址: https://webptt.com/m.aspx?n=bbs/R_Language/M.1562119714.A.431.html