###R物件.rds R物件有檔案小與讀取快速的優點,如果在R程式處理資料後必須儲存一份以供後續分析 ... 利用上述日期向量資料,搭配迴圈,依序設定since 和until參數. ... <看更多>
r 迴圈 存檔 在 [程式] 矩陣輸出及輸出存檔問題for R - Statistics | PTT Web 的推薦與評價
[程式]矩陣輸出及輸出存檔問題forR@statistics,共有11則留言,3人參與討論 ... 4 F 推gsuper:在回圈中加入一行print(class(aaa(data[i,],data2[i,]))) ... ... <看更多>
r 迴圈 存檔 在 [求助]迴圈中使用fopen、textscan - 看板MATLAB - PTT數位 ... 的推薦與評價
... 此時filename為cell之型式F=strvcat(filename); %轉成char是為了方便後續存檔時製造檔案名稱file=F(n,:); %令file=第n個檔案名稱fileID = fopen((file),'r'); ... ... <看更多>
r 迴圈 存檔 在 [程式] R迴圈結果輸出的問題- 看板Statistics - 批踢踢實業坊 的推薦與評價
------------------------------------------------------------------------
[軟體程式類別]:
R
[程式問題]:
應該算是簡單的資料分析與迴圈。 t-test結果輸出
[軟體熟悉度]:
新手(不到1個月)
[問題敘述]:
我的data大概是長下面這個樣子↓
A10 A11 A12.....A19 A20 A21...........
[1] 1 1 2 1 4 5
[2]
~~~
~~~
[30]
我要做的是將每一列針對綠色字(A1X)和黃色字(A2X)為對比分別做t-test
我的想法是把最大的矩陣分割為兩個小矩陣(normal & tumor)
然後要求他以第一列對應第一列的方式來做t-test
換言之最後要做30個t-test
雖然也可以用最笨的方法慢慢寫:
n1=nomal[1,]
t1=tumor[1,]
t.test(n1,t1)
↑然後重複30遍
不過這樣的寫法實在太冗長了,
所以想試試看用迴圈的方式,但按照現在的寫法只能跑出最後一筆數據
想問問看要怎麼寫才能夠讓他每一筆資料都出現結果
我曾經試過用sink,append=TRUE的方式來寫,但是失敗了Orz
下面把我寫的程式附上,因為我只是個新手,可能有些看起來很奇怪的程式語法
如果有的話希望大家幫忙指正
或是有什麼更好的寫法也請大家不吝賜教<(_ _)>
謝謝大家。
[程式範例]:(有問題的部分是黃色的字)
rm(list=ls())
library(RODBC)
A <- odbcConnectExcel2007("D:/data/Anne.xlsx")
all <- sqlQuery(A,"select*from[Sheet1$]")
nomal <- all[c(2:133)]
tumor <- all[c(134:244)]
n <- NULL
n <- as.vector(n)
a <- NULL
a <- as.vector(a)
i=1
for(i in 1:30)
{
n = nomal[i,]
a = tumor[i,]
e <- t.test(n,a)
i=i+1
e
}
-----------------------------------------------------------------------------
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 118.165.130.45
所以很對不起,是不是能講得再清楚一些呢?
我的問題在於:
t-test的結果並非一個vector、matrix或是dataframe
這三者的相互轉換方法我知道,但t-test的結果要如何擺到矩陣裡頭呢?
需不需要做些資料類型上的轉換?
剛接觸R,還有很多搞不清楚的地方,有些問題可能很笨.....
還是非常謝謝您<(_ _)>
後來加上 nomal <- as.matrix(nomal[1:132])
tumor <- as.matrix(tumor[1:111]) 這兩行指令以後就轉為矩陣了
感謝指正<(_ _)>
※ 編輯: lovelySky 來自: 118.165.142.164 (02/17 14:27)
... <看更多>
相關內容