Question

我有一个如下所示的数据框：

structure(list(A = c(70, 70, 70, 70, 70, 70), T = c(0.1, 0.2, 
0.3, 0.4, 0.5, 0.6), X = c(434.01, 434.01, 434.75, 434.75, 434.75, 
434.01), Y = c(454.92, 454.92, 454.92, 454.92, 454.18, 454.92
), V = c(0, 0, 21.128, 0, 14.94, 14.94), thetarad = c(0.151841552716899, 
0.151841552716899, 0.150990672182432, 0.150990672182432, 0.150177486839524, 
0.151841552716899), thetadeg = c(8.69988012340509, 8.69988012340509, 
8.6511282599214, 8.6511282599214, 8.6045361718215, 8.69988012340509
)), .Names = c("A", "T", "X", "Y", "V", "thetarad", "thetadeg"
), row.names = 1423:1428, class = "data.frame")

我希望以30秒的间隔对R中的特定时间点进行子集化。我可以通过手动对每个我想要的时间点进行子集来实现：

a1=subset(binA, T==0.1)
a2=subset(binA, T==30)
a3=subset(binA, T==60)
a4=subset(binA, T==90)
a5=subset(binA, T==120)
a6=subset(binA, T==150)
a7=subset(binA, T==180)
a8=subset(binA, T==210)
a9=subset(binA, T==240)
a10=subset(binA, T==270)
a11=subset(binA, T==300)
a12=subset(binA, T==330)
a13=subset(binA, T==360)
a14=subset(binA, T==390)
a15=subset(binA, T==420)
a16=subset(binA, T==450)
a17=subset(binA, T==480)
a18=subset(binA, T==510)
a19=subset(binA, T==540)
a20=subset(binA, T==570)
a21=subset(binA, T==599.5)

我尝试使用sapply和seq功能进行子集化，但结果令人困惑。我还想计算每个数据子集中的唯一A.我也知道我可以使用count包中的plyr函数来执行此操作。

a1=count(unique(subset(binA, T==0.1)))

但是计数将适用于一个数据框而不是多个数据框（如果我错了，请纠正我）。我还想为每个子集采用thetadeg的方法（这应该只在一个数据帧中很容易获得）。所以我需要有关如何编写具有特定seq点的函数的帮助。

我知道这个问题很简单，但我们将不胜感激。

由于

Answer 1

假设数据在df数据框中，那么试试这个：

sapply(c(0.1,seq(30,599,30),599.5),
       function(x)
         length(unique(df[ df$T==x, "A"])))

Answer 2

您应该能够使用以下代码来获得所需内容。这不是0.1和599.5，但应该易于操作。

timeintervals <- seq(0,600, 30)
for(i in 1:length(timeintervals)
{
  # create the subsets for each time interval
  assign(
    paste0("a",i),
    df[df$T == timeintervals[i],]
    )

  # get all unique As
  assign(
    paste0("b",i),
    unique(df[df$T == timeintervals[i],"A"])
  )

}

Answer 3

我认为您想要的功能是split：

 subsetted.by.T <- split(dfrm, dfrm$T)
lapply(subsetted.by.T, nrow)

$`0.1`
[1] 1

$`0.2`
[1] 1

$`0.3`
[1] 1

$`0.4`
[1] 1

$`0.5`
[1] 1

$`0.6`
[1] 1

> subsetted.by.T[[1]]
      A   T      X      Y V  thetarad thetadeg
1423 70 0.1 434.01 454.92 0 0.1518416  8.69988

如果您想为这些单独的项目命名，那么names<-功能将是合适的：

names(subsetted.by.T) <- paste0("a", seq(length(subsetted.by.T) ) )

如果“T”列的值有些不规则，那么使用cut在常规中断处创建类别可能对分割有用。如果“T”实际上是时间值，则可以澄清问题。目前它是一个“数字”值，但是有日期时间类的剪切方法。

Answer 4

如果只是为了获得平均值，唯一计数等，你不需要子集。还有一件事，id T因素是连续的你需要制作箱子吗？在这里我假设因素

这是plyr

的一种方法

ddply(df,~T,summarise,l=length(unique((A))))
ddply(df,~T,summarise,m=mean(thetadeg))

在R中的特定时间点序列中进行子集化，我可以使用seq吗？

4 个答案: