多圖功能R

Question

我是R新手，一直試圖通過谷歌和論壇在過去一週內找到解​​決方案。 我的問題：我有一個數據集，我需要繪製年齡。在40個不同的條件下，有超過1000個變量進行不同的測量。 是這樣的：

Age Variables1 Variable2 (....) Variable1000 > 
|  | 
|  | 
v  v 

我需要做的是劇情的條件（年齡）對各變量和輸出不同地塊的列（這一切僅僅是散點圖）。更重要的是，我希望輸出僅限於具有正趨勢線係數的變量。

所以目前我有這個非常醜陋的代碼，基本上是我真正需要的草稿。

plotest <- function(lung){ 
    # need to add the condition of abline function coefficient > 0 before plotting  
    plot(lung$Age, lung$hsa.let.7a.1, xlab = "Age", ylab = "miRNA") 
    abline(lm(lung$hsa.let.7a.1 ~ lung$Age), col= "red") 
    return(plot) 
} 
par(mfrow=c(2,2)) 
for (i in lung{plotest(i)}) 

我知道這是大多數錯誤。所以對於可怕的一切都很抱歉。

任何人都可以引導我到任何來源，我可能忽略瞭如何在如此大的數據集中指定範圍？和函數語法？我已經做了一些Python，但發現R設定爲多在這方面更加混亂...

感謝所有， 保羅

Answer 1

這應該是非常接近你問什麼，雖然你在做什麼去做1000張圖是超出我的。

# make up some data 
x <- seq(1,10,len=100) 
set.seed(1) # for reproducible example 
df <- data.frame(x,y1=1+2*x+rnorm(100), 
        y2=3-4*x+rnorm(100), 
        y3=2+0.001*x+rnorm(100)) 

# this does the work... 
lapply(colnames(df)[-1],function(col){ 
    form <- formula(paste(col,"x",sep="~")) 
    fit <- lm(form,df) 
    if (coef(fit)[2] >0) { 
    plot(form,df) 
    abline(fit) 
    } 
}) 

你的代碼是不是已經不遠了。此示例除第一個列名（colnames(df)[1]）之外的所有列名都將其傳遞給該函數。該函數使用列名稱和第一列的名稱創建一個公式變量，調用lm(...)，檢查x的係數是否大於0，如果是，則繪製數據和最佳擬合線。

查閱關於formula(...)，lm(...)和coef(...)的文檔。請注意，此示例有一個變量y3，斜率爲正數，但與0不顯着不同。您應該考慮如何處理這種情況。