My goal is to store the DWT coefficients from a number of time series in such a
way to be able to extract subsets satisfying giben conditions.
I have defined the following 11 matrices whose x-axis represent time intervals
and y-axis represent the number of time series:
MS <- 11 d1.mat <- matrix (data=list(), nrow=TotNumCycles,
ncol=2^(MS-1))
d2.mat <- matrix (data=list(), nrow=TotNumCycles, ncol=2^(MS-3))
d3.mat <- matrix (data=list(), nrow=TotNumCycles, ncol=2^(MS-3))
d4.mat <- matrix (data=list(), nrow=TotNumCycles, ncol=2^(MS-4))
d5.mat <- matrix (data=list(), nrow=TotNumCycles, ncol=2^(MS-5))
d6.mat <- matrix (data=list(), nrow=TotNumCycles, ncol=2^(MS-6))
d7.mat <- matrix (data=list(), nrow=TotNumCycles, ncol=2^(MS-7))
d8.mat <- matrix (data=list(), nrow=TotNumCycles, ncol=2^(MS-8))
d9.mat <- matrix (data=list(), nrow=TotNumCycles, ncol=2^(MS-9))
d10.mat <- matrix (data=list(), nrow=TotNumCycles, ncol=2^(MS-10))
d11.mat <- matrix (data=list(), nrow=TotNumCycles, ncol=2^(MS-11))
Each di.mat (i=1,2,...,11) matrix slot is assigned a DWT coefficient modulus
and the signal and cycle number it belongs to:
coef <- abs (as.numeric (aa.dwt$data[[n]])) #GET DWT COEFFICIENTS
..........................................................
for (k in 1:length(coef)) {
d7.mat [[nCycle ,k]] <- list (data=coef
[k],variable=paste(insig,":",cyc,sep=""))
}
As a consequence d7.mat contains the following:
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,] List,2 List,2 List,2 List,2 NULL NULL NULL NULL NULL NULL My
goal is to store the DWT coefficients from a number of time series in such a
way to be able to extract subsets satisfying giben conditions.
I have defined the following 11 matrices whose x-axis represent time intervals
and y-axis represent the number of time series:
MS <- 11 d1.mat <- matrix (data=list(), nrow=TotNumCycles,
ncol=2^(MS-1))
d2.mat <- matrix (data=list(), nrow=TotNumCycles, ncol=2^(MS-3))
d3.mat <- matrix (data=list(), nrow=TotNumCycles, ncol=2^(MS-3))
d4.mat <- matrix (data=list(), nrow=TotNumCycles, ncol=2^(MS-4))
d5.mat <- matrix (data=list(), nrow=TotNumCycles, ncol=2^(MS-5))
d6.mat <- matrix (data=list(), nrow=TotNumCycles, ncol=2^(MS-6))
d7.mat <- matrix (data=list(), nrow=TotNumCycles, ncol=2^(MS-7))
d8.mat <- matrix (data=list(), nrow=TotNumCycles, ncol=2^(MS-8))
d9.mat <- matrix (data=list(), nrow=TotNumCycles, ncol=2^(MS-9))
d10.mat <- matrix (data=list(), nrow=TotNumCycles, ncol=2^(MS-10))
d11.mat <- matrix (data=list(), nrow=TotNumCycles, ncol=2^(MS-11))
Each di.mat (i=1,2,...,11) matrix slot is assigned a DWT coefficient modulus
and the signal and cycle number it belongs to:
coef <- abs (as.numeric (aa.dwt$data[[n]])) #GET DWT COEFFICIENTS
..........................................................
for (k in 1:length(coef)) {
d7.mat [[nCycle ,k]] <- list (data=coef
[k],variable=paste(insig,":",cyc,sep=""))
}
As a consequence d7.mat contains the following:
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,] List,2 List,2 List,2 List,2 NULL NULL NULL NULL NULL NULL
[2,] List,2 List,2 List,2 List,2 NULL NULL NULL NULL NULL NULL
[3,] List,2 List,2 List,2 List,2 NULL NULL NULL NULL NULL NULL
[4,] List,2 List,2 List,2 List,2 NULL NULL NULL NULL NULL NULL
[5,] List,2 List,2 List,2 List,2 NULL NULL NULL NULL NULL NULL
[6,] List,2 List,2 List,2 List,2 NULL NULL NULL NULL NULL NULL
[7,] List,2 List,2 List,2 List,2 NULL NULL NULL NULL NULL NULL
[8,] List,2 List,2 List,2 NULL NULL NULL NULL NULL NULL NULL
[9,] List,2 List,2 List,2 List,2 NULL NULL NULL NULL NULL NULL
................................................................................................................
> d5.mat[[3,4]]
$data
[1] 0.6290917
$variable
[1] "10146:3"
> d5.mat[[6,4]]
$data
[1] 1.187133
$variable
[1] "10146:6"
> d5.mat[[8,4]]
NULL
> d5.mat[[10,4]]
$data
[1] 0.7789786
$variable
[1] "10146:10"
> d5.mat[[50,4]]
$data
[1] 0.8909976
$variable
[1] "10146:50"
How can I extract, for instance, all the data from the i_th column of the
matrix that satisfy a given conditon?
In the case of a matrix A made up of numbers I would write :
which (A [,i] > 0.1) ...?
Thank you in advance for your help.
Maura
tutti i telefonini TIM!
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
