2015-09-04 85 views
1

Hei, 我想創建我的數據的散點圖。我有一個數據框,其中包括(不僅)一年的四分之一,溫度,降水和一個物種的gam_indices列。 現在我想暗算的gam_index(X)溫度(Y),但前提是quartal爲1 到目前爲止,我有基本的散點圖配方如果滿足某些條件,則R圖

plot(Archon_noduplicate$regional_gam, Archon_noduplicate$prev.temp, main="Archon_noduplicate", 
xlab="logAUC", ylab="prev.temp", pch=19) 

,但沒有真正的想法有關一個簡短的方法來包含條件。當然,有一種方法可以首先創建一個只包含四分之一溫度值的新列,然後繪製該列,但由於我有很多種類和很多變量,所以很容易找到一些不那麼單調乏味的東西。 我的數據,如果你需要一些例子,看起來是這樣的:

transect_id year day month  LST precipitation Quarter  SumPre average.temp MinTemp MaxTemp 
    1  TR001 2011 1  1 15.70134  0.000000  1 161.9494290  14.69150 10.28747 18.95994 
    2  TR001 2011 91  4 19.53845  1.858437  2 106.3809708  23.56554 16.58574 29.81955 
    3  TR001 2011 182  7 28.96053  0.000000  3 0.9290609  29.62564 26.22075 31.78875 
    4  TR001 2011 274 10 26.46966  0.000000  4 66.1280496  18.84587 12.18172 29.72368 
    5  TR001 2012 1  1 13.17917  2.534553  1 154.5095844  13.38382 10.64767 17.80582 
    6  TR001 2012 92  4 17.18347  0.000000  2 9.4365772  24.77507 17.18347 30.88236 
     prev.temp prev.Precip   species regional_gam prop_pheno_sampled 
    1 16.22194   0 Archon apollinus  1.780317   0.4154008 
    2 24.81180   0 Archon apollinus  1.780317   0.4154008 
    3 30.78058   0 Archon apollinus  1.780317   0.4154008 
    4 22.72036   0 Archon apollinus  1.780317   0.4154008 
    5 14.69150   0 Archon apollinus  1.103804   0.4413598 
    6 23.56554   0 Archon apollinus  1.103804   0.4413598 

我有種陷在這裏,所以我很感激任何輸入! 謝謝!

回答

3

嘗試:

with(Archon_noduplicate[Archon_noduplicate$Quarter==1,], plot(regional_gam, prev.temp, main="Archon_noduplicate", xlab="logAUC", ylab="prev.temp", pch=19)) 
+0

酷,的作品!非常感謝! –

相關問題