R, если выполнены определенные условия

Question

Hei, Я хочу создать диаграмму рассеяния данных. У меня есть dataframe, который включает (а не только) квартал года, температуру, осадки и gam_indices для вида в виде столбцов. Теперь я хотел бы построить температуру (у) против gam_index (х), но только если это 1. Квартал До сих пор у меня есть основной диаграммы рассеяния-формула

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 

Я вроде застрял здесь, так что я благодарен за любой вклад! Спасибо!

Answer 1

Try:

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