天天看點

r語言 col_R語言代碼日志之讀取寫入資料

我已經盡最大的可能不談政治,但是,被趕出學校,又沒去醫院,都是因為疫情。現在做一些R語言心得筆記,感謝頭條給的這個免費的網頁,如果疫情持續到3月底,那估計我都可以寫出一個網站了。

R語言,和其他的C、B、++還有什麼PH什麼的到底有什麼關系,我不清楚,但是我最清楚的是,這R語言,是一個統計學軟體,是以code為基礎的,和SPASS(這個軟體衆所周知,是收費的,是圖形界面,傻瓜式操作,容易上手,能夠輸出各種圖像和資料,性能不是一般強大,但是就有一個問題,收費的)不一樣的操作方式。

R語言最厲害的,最吸引人的,不是因為它功能強大,而是因為它免費。配套RGUI軟體使用的是Rstudio界面分割軟體。

r語言 col_R語言代碼日志之讀取寫入資料

蒼天饒過誰,即便他們老态龍鐘換成狼心狗肺(移植)的,也難逃一死。然而那些罪大惡極的他們死後的一百年不被鞭屍,他們的基因還在耀武揚威。這就是我們用R語言研究基因表達的必要性,我們需要鑒别那些特有的,具有差別的基因。

給資料框加上列名test01就是your data.frame

colnames(test01)

儲存資料成TXT檔案活着CSV檔案(txt檔案,是個電腦或者手機都能打開,CSV檔案用辦公軟體excel等表格軟體打開。)

write.table(test01, file="GOALL.txt",sep ="", row.names =TRUE, col.names =TRUE, quote =TRUE)

write.table(test01, file="GOALL.csv",sep ="", row.names =TRUE, col.names =TRUE, quote =TRUE)

其實那個解放黑奴的林肯是一個很會演講的人,林肯講了一個故事:當時南北雙方都面臨嚴峻的問題,奴隸放假時間太長,還天天要吃糧,要取暖,甚至還要吵吵嚷嚷,長時間不幹活,奴隸主們覺得白養了奴隸,奴隸太多多的休息,奴隸主太痛心了,休假以後要讓奴隸多幹活、不休假補償奴隸主。以保證抗擊北方政權的物資儲備。沒想到,到了交戰的後期,林肯被逼的沒辦法,欲取得奴隸的倒戈,以解決膠着的戰争,林肯才成了那個解放黑奴的人。其實沒有人願意解放奴隸。解放奴隸隻不過是曆史發展的副産物。奴隸是曆史上形式不同,概念相同的元素。沒有奴隸這個元素,可能會少了一點曆史,但是不會改變曆史。

r語言 col_R語言代碼日志之讀取寫入資料

對資料進行去重:

##############################################################################################

#####################################################以上是去除@後但是未去重的GLP和表達譜

exprSet9 %

#合并探針的資訊

#inner_join(probe2symbol_df,by="probe_id") %>%

#去掉多餘列

# select(-probe_id) %>%

#重新排列

select(GENE_SYMBOL,everything()) %>%

#求出平均數(這邊的.代表上一步的資料,實際上是exprSet3)

mutate(rowMean =rowMeans(.[grep("GSM", names(.))])) %>%

#把表達量的平均值按從大到小排序

arrange(desc(rowMean))# %>%

# symbol留下第一個

distinct(SYMBOL,.keep_all = T) %>%

#反向選擇去除rowMean這一列

select(-rowMean)

事實證明,在R語言裡面,如果隻是改變行名、列名,還遠遠不夠,因為這樣的資料實際上上是沒有被系統儲存的。系統是不認的。R語言的計算法則,就根本不知道存在着你改名的那列的資料,或者你修改過行名的那列資料,你隻是作為人類,自己認為修改了,但是機器沒感覺到。着就是人類大腦思維方式和機器思考方式的異同。哪天如果機器思維方式不這麼幹了,大約,人類就被機器奴役了。顯而易見的,将來的世界,必定是智能機器的世界,我們所有的工作、努力,其實就是為了制造控制人類的終極大boss。當然,伺候小boss和伺候大boss,對于普通人沒什麼兩樣。因為這些boss,壓根不會想你。

r語言 col_R語言代碼日志之讀取寫入資料