R_Language 板


LINE

你在推文中的補充說明還是不清不楚。 不過我猜你要的只是 eval 硬幹。 a <- rnorm(10) b <- rnorm(10) c <- rnorm(10) d <- data.frame(a, b, c) txt <- paste0(names(d)[1], "~", names(d)[2], "+", names(d)[3]) f <- eval(parse(text = txt)) summary(lm(f, data = d)) 這個例子中,只要被讀進去的 data frame 叫作 d,且其中欄位順序不變, 則欄位的名字(在此例就是 a、b 和 c)是無所謂的。 重點在倒數第二列和第三列,看看是不是你要的吧。 如果我猜錯你要的,請你把問題說明得更具體一點。 這裡熱心的人很多, 但請不要讓熱心的人花大部份的時間在猜發問者的問題到底是什麼啦。 ※ 引述《w2a3n4g5 (MORRIS.W)》之銘言: : - 問題: 當你想要問問題時,請使用這個類別 : [問題類型]: : 程式諮詢(我想用R 做某件事情,但是我不知道要怎麼用R 寫出來) : [軟體熟悉度]: : 入門(寫過其他程式,只是對語法不熟悉) : [問題敘述]: : 大家好我又來問問題了QQ : 目前遇到的問題如下 : 先讀入資料檔 , 再透過commandline給定參數去決定自變數 ,應變數的位置 : 程式碼大略如下 : read.table(data,"路徑") : y <- data[[k1]] : x1 <- data[[k2]] : x2 <- data[[k3]] : y:應變數 , x1 , x2 皆為連續型自變數 , 不須轉換成分類變數 : 其中k1,k2,k3行的名稱分別為 weigh , A ,B : output <- lm(y ~ x1 + x2,data) : summary(output) : 則會出現下列結果 : Estimate Std. Error t value Pr(>|t|) : (Intercept) 26.23730 5.93517 4.421 1.08e-05 *** : data[, k2] 0.01520 0.01138 1.336 0.1819 : data[, k3] -0.90883 8.18376 -0.111 0.9116 : 但我希望summary後的結果像這樣 : Estimate Std. Error t value Pr(>|t|) : (Intercept) 26.23730 5.93517 4.421 1.08e-05 *** : A 0.01520 0.01138 1.336 0.1819 : B -0.90883 8.18376 -0.111 0.9116 : 就是如同 lm (weigh ~ A+B , data=data)的效果 : 請問我再 lm() 中 該怎麼寫呢? : 已經試過很多方式例如 names(data)[k1] 之類的方式lm當中 , 但還是失敗了 : 也看過help(lm)了 ,不太了解內部放的y~x屬於何種型態 : 因為檔案是從外不讀入所以行的名稱不固定不能寫死QQ : 第二個問題是 , 假設我將檔案輸出了 , 程式碼大略如下 : write.table(summary(output),"路徑") : 假設輸出的是txt檔 , 前面都會有數字1,2...... 就像下面那樣 : "x" : "1" "" : "2" "Call:" : "3" "lm(formula = data[, f] ~ data[, c1] + data[, c2] + data[, c3], " : "4" " data = data : 請問是否有辦法再輸出時去除掉那些數字呢? : 謝謝各位QQ --



※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.128.117.75
※ 文章網址: http://webptt.com/m.aspx?n=bbs/R_Language/M.1411538660.A.C46.html ※ 編輯: andrew43 (140.128.117.75), 09/24/2014 14:06:54
1F:推 koai: 推! ^^ 09/24 14:09
2F:→ obarisk: 其實程式不熟的人,不是很好瞭解expression和character 09/24 15:33
3F:→ obarisk: 的差別,特別是人門書又不會多加著墨 09/24 15:33
4F:→ andrew43: 以前學matlab時,大部份的教材對eval還蠻著重的。 09/24 16:36
5F:→ andrew43: 否則我也不會想到用eval硬幹。大多數R的教本沒強調這個. 09/24 16:39
6F:→ celestialgod: 可以直接as.formula 不用用eval 09/24 18:13
7F:推 w2a3n4g5: 抱歉我真的努力把問題說清楚了QQ ,因為在學校沒學過R 09/25 09:48
8F:→ w2a3n4g5: 只學過C,java 但專題老師是統計專長的希望我們寫R 09/25 09:49
9F:→ w2a3n4g5: 所以只能慢慢摸索 , 不好意思~~~~~ 09/25 09:52
10F:→ w2a3n4g5: 另外成功了!! 謝謝大家 09/25 09:53
11F:推 gsuper: 邪惡的貼上 09/25 13:07
12F:→ obarisk: matlab應該沒有formula這個型態XD 09/25 14:58







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