R_Language 板


LINE

[問題類型]: 程式諮詢(我想用R 做某件事情,但是我不知道要怎麼用R 寫出來) 不曉得標題這樣下對不對,我想要計算多比不等長的資料(存在list) 舉例:共四筆資料,我想要運算1 vs 2,3,4 ; 2 vs 3,4 ; 3 vs 4 (vs:比較或是運算,並且重覆的不算) 效能諮詢(我想讓R 跑更快) 我已經有使用多層for loop將上述描述實作完畢,但是迴圈多層,運算效能非常差... [軟體熟悉度]: 使用者(已經有用R 做過不少作品) 算是有學了一陣子,但感覺程式邏輯還是很差 [問題敘述]: 如上述,我要運算多筆資料,且作兩兩比較,目的是將計算的結果使用矩陣表示 舉例:共四筆資料,1,2,3,4 想要兩兩比較,並以矩陣方式表示 (如下圖) 1 2 3 2 1vs2 3 1vs3 2vs3 4 1vs4 2vs4 3vs4 [程式範例]: 我直接貼程式碼,並以註解的方式解釋 實際上我要做的是修改DTW這個套件,並已經將裡面的內容改成我要的運算方式 sample_list<- list(1,2,3,4) # 有幾筆資料就設定list放幾筆的資料 tmp1 <-list() i <-1 for(i in 1:3){ tmp <- lapply(sample_list, function(x) x+i) tmp1[[i]] <- unlist(tmp) tmp1[[i]][which(tmp1[[i]]>4)] <- NA } # for 1 到 3(4筆資料-1自己),逐步隨著i值相加,最後將大於4的設為NA tmp1 <-lapply(tmp1, function(x) x[!is.na(x)]) # 把NA 刪掉 tmp_vector<-vector() for (i in 1:3){ for(k in 1:length(unlist(tmp1[i]))){ for(j in unlist(tmp1[i])[k]){ kai<-kai_dtw(x=unlist(data.list[i]) ,y=unlist(data.list[j]), g=0, step.pattern = symmetric1) tmp<- kai$distance tmp_vector <- c(tmp_vector,tmp) } } } # 多層迴圈,達到我要兩兩比較的意境,當中的kai_dtw是我修改演算法的函式 它對應的x與y一定要unlist,後面就是將其取distance存在vector中 tmp_vector1 <- matrix(nrow = 4,ncol = 4) tmp_vector1 <- as.dist(tmp_vector1) tmp_vector1[1:6] <- tmp_vector # 創一個空的矩陣(tmp_vector1)m,並轉算距離矩陣的形式,將前述得到的vector 填入,得到我要的結果 [關鍵字]: dtw , 距離矩陣 ,兩兩比較 以上 就是我的問題敘述,不曉得有沒有表達清楚 基本上就是多層迴圈太慢了,我也正在嘗試使用lapply家族進行運算 謝謝~~ --



※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.116.86.113
※ 文章網址: https://webptt.com/m.aspx?n=bbs/R_Language/M.1504073346.A.4C9.html ※ 編輯: pk790127 (140.116.86.113), 08/30/2017 14:20:04 ※ 編輯: pk790127 (140.116.86.113), 08/30/2017 14:21:59
1F:→ andrew43: 我覺得瓶頸應該是你的演算法不夠快,倒不是填值的loop 08/30 14:27
2F:→ celestialgod: preallocation 08/30 20:10
3F:→ celestialgod: 還有 你有必要先搞動[跟[[.... 08/30 20:10
4F:→ celestialgod: 懂 08/30 20:10
5F:→ celestialgod: kai_dtw這個函數怎麼實做講出來比較快 08/30 20:11
6F:→ celestialgod: 這個函數可能是最大的瓶頸 08/30 20:11







like.gif 您可能會有興趣的文章
icon.png[問題/行為] 貓晚上進房間會不會有憋尿問題
icon.pngRe: [閒聊] 選了錯誤的女孩成為魔法少女 XDDDDDDDDDD
icon.png[正妹] 瑞典 一張
icon.png[心得] EMS高領長版毛衣.墨小樓MC1002
icon.png[分享] 丹龍隔熱紙GE55+33+22
icon.png[問題] 清洗洗衣機
icon.png[尋物] 窗台下的空間
icon.png[閒聊] 双極の女神1 木魔爵
icon.png[售車] 新竹 1997 march 1297cc 白色 四門
icon.png[討論] 能從照片感受到攝影者心情嗎
icon.png[狂賀] 賀賀賀賀 賀!島村卯月!總選舉NO.1
icon.png[難過] 羨慕白皮膚的女生
icon.png閱讀文章
icon.png[黑特]
icon.png[問題] SBK S1安裝於安全帽位置
icon.png[分享] 舊woo100絕版開箱!!
icon.pngRe: [無言] 關於小包衛生紙
icon.png[開箱] E5-2683V3 RX480Strix 快睿C1 簡單測試
icon.png[心得] 蒼の海賊龍 地獄 執行者16PT
icon.png[售車] 1999年Virage iO 1.8EXi
icon.png[心得] 挑戰33 LV10 獅子座pt solo
icon.png[閒聊] 手把手教你不被桶之新手主購教學
icon.png[分享] Civic Type R 量產版官方照無預警流出
icon.png[售車] Golf 4 2.0 銀色 自排
icon.png[出售] Graco提籃汽座(有底座)2000元誠可議
icon.png[問題] 請問補牙材質掉了還能再補嗎?(台中半年內
icon.png[問題] 44th 單曲 生寫竟然都給重複的啊啊!
icon.png[心得] 華南紅卡/icash 核卡
icon.png[問題] 拔牙矯正這樣正常嗎
icon.png[贈送] 老莫高業 初業 102年版
icon.png[情報] 三大行動支付 本季掀戰火
icon.png[寶寶] 博客來Amos水蠟筆5/1特價五折
icon.pngRe: [心得] 新鮮人一些面試分享
icon.png[心得] 蒼の海賊龍 地獄 麒麟25PT
icon.pngRe: [閒聊] (君の名は。雷慎入) 君名二創漫畫翻譯
icon.pngRe: [閒聊] OGN中場影片:失蹤人口局 (英文字幕)
icon.png[問題] 台灣大哥大4G訊號差
icon.png[出售] [全國]全新千尋侘草LED燈, 水草

請輸入看板名稱,例如:Boy-Girl站內搜尋

TOP