1

我有一个简化如下的数据框： 

head(dendro)
  X            DateTime   ID diameter dendro ring DOY month mday year Rain_mm_Tot Through_Tot temp
1 1 2012-06-21 13:45:00 r1_1     5482      1    1 173     6   22  113          NA          NA   NA
2 2 2012-06-21 13:45:00 r2_3       NA      3    2 173     6   22  113          NA          NA   NA
3 3 2012-06-21 13:45:00 r1_2     5534      2    1 173     6   22  113          NA          NA   NA
4 4 2012-06-21 13:45:00 r2_4       NA      4    2 173     6   22  113          NA          NA   NA
5 5 2012-06-21 13:45:00 r1_3     5606      3    1 173     6   22  113          NA          NA   NA
6 6 2012-06-21 13:45:00 r2_5       NA      5    2 173     6   22  113          NA          NA   NA

数据框首先被“ID”分割，所以它是一个 ID 列表之后我应用了一个包含循环的函数，结果是一个新列“Diameter2”，我希望从函数中得到结果，它可以工作好的：

dendro_sp <- split(dendro, dendro$ID)

library(changepoint)

dendro_sp <- lapply(dendro_sp, function(x){

 x <- subset(x, !is.na(diameter))
 cpfit <- cpt.mean(x$diameter, method="BinSeg")

 x$diameter2 <- x$diameter
 cpts <- cpfit@cpts
 means <- param.est(cpfit)$mean
 meanZero <- means[1]

 for(i in 1:(length(cpts)-1)){
   x$diameter2[(cpts[i]+1):cpts[i+1]] <- x$diameter2[(cpts[i]+1):cpts[i+1]] + (meanZero - means[i+1])
 }

 return(x)
})
dendro2 <- do.call(rbind, dendro_sp)
rownames(dendro2) <- NULL

我的问题是我希望它有条件地应用它，例如应用到 r1_1 和 r1_3，并在新列“diameter2”中获取 r3 的“直径”值，而不是为其余 ID 应用该函数： 

ifelse(diameter$ID==c("r1_1","r1_3"), apply_the_function_to_r11_and_r13_to_calculate_diameter2, otherwise_write_diameter_value_in_diameter2_column)

请记住，数据框“dendro”是按 ID 拆分的，我不知道这对于定义多个 ID 的条件是否重要。谢谢

4

1 回答 1

0

我不确定我是否正确理解了这个问题。我试着回答。我假设您想将函数应用于“直径”数据帧的“直径”字段，以“ID”字段为条件，并在相应的直径 2 字段中返回结果。我不知道该功能是如何工作的，如果这不起作用，请原谅我。

选定字段

直径$diameter2[diameter$ID=="r1_1"|diameter$ID=="r1_3"]<- yourfun(diameter$diameter[diameter$ID=="r1_1"|diameter$ID=="r1_3"]

未选择的字段

直径$diameter2[直径$ID!="r1_1" & 直径$ID=="r1_3"]<- 直径$diameter[直径$ID=="r1_1"|diameter$ID=="r1_3"]

于 2013-05-02T16:37:47.090 回答