用R語言進行相關性分析:畫出散點圖思杯、擬合線及標注相關系數(shù)和P值

一. 讀取數(shù)據(jù)

因為我們一般自己錄的數(shù)據(jù)都是excel表格的贞铣,對于讀取excel闹啦,此處有一個小技巧

讀取excel方法
1.把excel文件另存為CSV文件
1png
2.png
3.png
4.png
2.按照讀取csv文件方法讀取
dat <- read.table("excel1.csv",sep=",",header = TRUE)
此名為excel1文件被R語言讀取

二. 利用ggplot2畫圖

1.加載ggplot2包

library(ggplot2)#加載ggplot2包
注意此處初學者容易出現(xiàn)一個錯誤,ggplot并不是R包辕坝,而是ggplot2包中的函數(shù)窍奋,所以你如果加載或安裝ggplot那肯定會是提示無法加載或安裝的
問題是它還會加載進度條,最后安裝失敗酱畅,所以反反復復會讓你懷疑人生琳袄,我就這樣試過
QQ截圖20190322211336.jpg
QQ截圖20190322211351.jpg

2.畫出散點圖

ggplot(data=dat, aes(x=gene, y=gene2))+geom_point(color="red")
#data=后跟需要畫圖的數(shù)據(jù)的文件名
#X=后跟作為X軸的數(shù)據(jù)的那一列的列名
#Y=后跟作為Y軸的數(shù)據(jù)的那一列的列名
#geom_point函數(shù)是個性化設置散點圖點的形狀,顏色纺酸,大小等窖逗,此處只設置了顏色,有需要可自行加入

Rplot.jpeg

3.畫出擬合曲線

ggplot(data=dat, aes(x=gene, y=gene2))+geom_point(color="red")+stat_smooth(method="lm",se=FALSE)
#stat_smooth是畫擬合曲線的函數(shù)
#se=FALSE意思為不畫出置信區(qū)間
有se=FALSE的運行結果

無se=FALSE的運行結果

4.進行相關性分析

library(ggpubr)
ggplot(data=dat, aes(x=gene, y=gene2))+geom_point(color="red")+stat_smooth(method="lm",se=FALSE)++stat_cor(data=dat, method = "pearson")
#stat_cor(data=dat, method = "pearson")意為用pearson相關進行相關性分析吁峻,可以自行更改方法
具有散點圖滑负、擬合線、相關系數(shù)用含、P值的相關分析圖片就出來了

注意此處有一個易錯點矮慕,就是忘記加載ggpubr,因為stat_cor是ggpubr包中的啄骇,如果未加載ggpubr包痴鳄,則會出現(xiàn)


忘記加載ggpubr
此處有一個小技巧分享,當我們新手發(fā)現(xiàn)某個函數(shù)不可用時缸夹,利用痪寻??stat_cor去查說明書,可以看到此函數(shù)屬于哪個包虽惭,再加載或安裝這個包橡类,函數(shù)不可用這個問題就解決了
利用??stat_cor去查說明書

得知函數(shù)屬于哪個包

最后

感謝jimmy的生信技能樹團隊芽唇!
感謝導師岑洪老師顾画!
感謝健明取劫、孫小潔,慧美等生信技能樹團隊的老師一路以來的指導和鼓勵研侣!
最后編輯于
?著作權歸作者所有,轉載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末谱邪,一起剝皮案震驚了整個濱河市,隨后出現(xiàn)的幾起案子庶诡,更是在濱河造成了極大的恐慌惦银,老刑警劉巖,帶你破解...
    沈念sama閱讀 217,907評論 6 506
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件末誓,死亡現(xiàn)場離奇詭異扯俱,居然都是意外死亡,警方通過查閱死者的電腦和手機基显,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,987評論 3 395
  • 文/潘曉璐 我一進店門蘸吓,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人撩幽,你說我怎么就攤上這事库继。” “怎么了窜醉?”我有些...
    開封第一講書人閱讀 164,298評論 0 354
  • 文/不壞的土叔 我叫張陵宪萄,是天一觀的道長。 經(jīng)常有香客問我榨惰,道長拜英,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 58,586評論 1 293
  • 正文 為了忘掉前任琅催,我火速辦了婚禮居凶,結果婚禮上,老公的妹妹穿的比我還像新娘藤抡。我一直安慰自己侠碧,他們只是感情好,可當我...
    茶點故事閱讀 67,633評論 6 392
  • 文/花漫 我一把揭開白布缠黍。 她就那樣靜靜地躺著弄兜,像睡著了一般。 火紅的嫁衣襯著肌膚如雪瓷式。 梳的紋絲不亂的頭發(fā)上替饿,一...
    開封第一講書人閱讀 51,488評論 1 302
  • 那天,我揣著相機與錄音贸典,去河邊找鬼视卢。 笑死,一個胖子當著我的面吹牛廊驼,可吹牛的內(nèi)容都是我干的据过。 我是一名探鬼主播颊埃,決...
    沈念sama閱讀 40,275評論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼蝶俱!你這毒婦竟也來了?” 一聲冷哼從身側響起饥漫,我...
    開封第一講書人閱讀 39,176評論 0 276
  • 序言:老撾萬榮一對情侶失蹤榨呆,失蹤者是張志新(化名)和其女友劉穎,沒想到半個月后庸队,有當?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體积蜻,經(jīng)...
    沈念sama閱讀 45,619評論 1 314
  • 正文 獨居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 37,819評論 3 336
  • 正文 我和宋清朗相戀三年彻消,在試婚紗的時候發(fā)現(xiàn)自己被綠了竿拆。 大學時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 39,932評論 1 348
  • 序言:一個原本活蹦亂跳的男人離奇死亡宾尚,死狀恐怖丙笋,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情煌贴,我是刑警寧澤御板,帶...
    沈念sama閱讀 35,655評論 5 346
  • 正文 年R本政府宣布,位于F島的核電站牛郑,受9級特大地震影響怠肋,放射性物質發(fā)生泄漏。R本人自食惡果不足惜淹朋,卻給世界環(huán)境...
    茶點故事閱讀 41,265評論 3 329
  • 文/蒙蒙 一笙各、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧础芍,春花似錦杈抢、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,871評論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至虏缸,卻和暖如春鲫懒,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背刽辙。 一陣腳步聲響...
    開封第一講書人閱讀 32,994評論 1 269
  • 我被黑心中介騙來泰國打工窥岩, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留,地道東北人宰缤。 一個月前我還...
    沈念sama閱讀 48,095評論 3 370
  • 正文 我出身青樓颂翼,卻偏偏與公主長得像晃洒,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子朦乏,可洞房花燭夜當晚...
    茶點故事閱讀 44,884評論 2 354

推薦閱讀更多精彩內(nèi)容