如何输出存入矩阵

Question

列表

数据：

x <- seq(0, 1, len = 1024) 
pos <- c(0.1, 0.13, 0.15, 0.23, 0.25, 0.40, 0.44, 0.65, 0.76, 0.78, 0.81) 
hgt <- c(4, 5, 3, 4, 5, 4.2, 2.1, 4.3, 3.1, 5.1, 4.2) 
wdt <- c(0.005, 0.005, 0.006, 0.01, 0.01, 0.03, 0.01, 0.01, 0.005, 0.008, 0.005) 
pSignal <- numeric(length(x)) 
for (i in seq(along=pos)) { 
    pSignal <- pSignal + hgt[i]/(1 + abs((x - pos[i])/wdt[i]))^4 
} 
df = as.data.frame(rbind(pSignal,pSignal,pSignal)) 
dflist=list(df,df,df) 

我试图运行此pracma包的findpeaks（）函数查找每行的局部最大值在每个data.frame列表，dflist。输出是N×4阵列。 N =峰值数量。所以在第一个data.frame的第一行，如果它找到4个峰值，它将是一个4x4矩阵。我的目标是将此函数循环遍历每个data.frame中的每一行，并将输出的矩阵存储在列表中。

我的代码：

## Find Peaks 
pks=list() 
for (i in 1:length(dflist)){ 
    for (j in 1:length(dflist[[i]])){ 
    row = dflist[[i]][j,] 
    temppks = findpeaks(as.vector(row,mode='numeric') 
         ,minpeakheight = 1.1,nups=2) 
    pks[i][[j]]=rbind(pks,temppks) 
    } 
} 

这似乎并没有做完全是我想要的太多。有任何想法吗？

Answer 1

应用的组合（）和sapply（）可以做的工作：

my.f.row <- function(row) findpeaks(as.vector(row,mode='numeric'), minpeakheight = 1.1, nups=2) 
sapply(dflist, function(df.i) apply(df.i, 1, my.f.row)) 

最终你要重组的结果。