0

我想将我的一个建模结果的输出分成某种 ID 列。由于我有不同的粒子轨迹（n = 15000）和不同的nrow间距，我不能假设它们之间的间隔相等（n_all = 1767333）。我唯一的分隔指标是以“0”开头的“时间”列。--> 因此，如果出现另一个“0.000”，则 ParticleID 应增加一...

我尝试了不同的 for 循环和 if-else 函数，但我无法这样做......我非常感谢你的帮助。下面的代码示例和数据示例。提前致谢。

其中一个尝试不起作用：

#initial seperation process
n_all <- nrow(all_trac)                                   #number of rows - 1767333 rows
n <- nrow(subset(all_trac, all_trac[,"IntTimeloc"] == 0)) #number of particles - 15806 particles
all_trac$Particle_ID <- 1:n_all

#start loop
for (j in 1:n_all)

  for (i in 1:n){

      if (all_trac$IntTimeloc[j] != 0) {

        all_trac$Particle_ID[j] = i

      }

      else

      all_trac$Particle_ID[j] = i+1

  } 

}
#end loop

到底应该怎么输出...

IntTimeloc  VXloc           VYloc           VZloc       Particle_ID
0.000   0.0009997907    9.038499e-04    -8.904433e-05   1
1.000   0.0009920447    9.680137e-04    -7.579094e-05   1
3.000   0.0009914209    9.718620e-04    -7.473217e-05   1
...
0.000   0.0008281278    -2.395466e-05   -1.537789e-04   2
9.000   0.0008310459    -2.691492e-05   -1.533394e-04   2
...
4

1 回答 1

0

试试这个例如：

transform(dat,id = cumsum(c(1,diff(IntTimeloc)<0)))

 IntTimeloc        VXloc         VYloc         VZloc Particle_ID id
1          0 0.0009997907  9.038499e-04 -8.904433e-05           1  1
2          1 0.0009920447  9.680137e-04 -7.579094e-05           1  1
3          3 0.0009914209  9.718620e-04 -7.473217e-05           1  1
4          0 0.0008281278 -2.395466e-05 -1.537789e-04           2  2
5          9 0.0008310459 -2.691492e-05 -1.533394e-04           2  2
于 2013-03-19T12:35:22.607 回答