2015-09-04 126 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

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

相关问题