2013-08-21 60 views
4

当我执行下面的我的“预测”数据集是正确填充:R hdf5数据集写入不正确?

library(rhdf5) 
library(forecast) 
library(sltl) 
library(tseries) 

fid <- H5Fcreate(output_file) 
## TODO: compute the order p 
p <- 4 
# write predictors 
h5createDataset(output_file, dataset="predictors", c(p, length(tsstl.remainder) - (p - 1)), storage.mode='double') 
predictors <- as.matrix(tsstl.remainder) 
for (i in 1:(p - 1)) {  
    predictors <- as.matrix(cbind(predictors, Lag(as.matrix(tsstl.remainder), i))) 
} 
predictors <- as.matrix(predictors[-1:-(p-1),]) 
head(predictors) 
h5write(predictors, output_file, name="predictors") 
H5Fclose(fid) 

生成的(正确)的输出head(predictors)是:

  [,1]  [,2]  [,3]  [,4] 
[1,] 0.3089645 6.7722063 5.1895389 5.2323261 
[2,] 8.7607228 0.3089645 6.7722063 5.1895389 
[3,] -0.9411553 8.7607228 0.3089645 6.7722063 
[4,] -14.1390243 -0.9411553 8.7607228 0.3089645 
[5,] -26.6605296 -14.1390243 -0.9411553 8.7607228 
[6,] -8.1293076 -26.6605296 -14.1390243 -0.9411553 

然而,当我看到它的结果是不纠正:

tsmatrix <- t(as.matrix(h5read(output_file, "predictors"))) 
head(tsmatrix) 

错误输出:

  [,1]  [,2]  [,3]  [,4] 
[1,] 0.3089645 8.760723 -0.9411553 -14.13902 
[2,] -26.6605296 -8.129308 -9.8687675 31.52086 
[3,] 54.2703126 43.902489 31.8164836 43.87957 
[4,] 22.1260636 36.733055 54.7064107 56.35158 
[5,] 36.3919851 25.193068 48.2244464 57.12196 
[6,] 48.0585673 72.402673 68.3265518 80.18960 

为什么我写的东西与我回来的东西不符?我再次检查并且hdfview HDF5查看器也为“预测者”数据集显示了这个不正确的值。

这里有什么问题?

+0

请提及您正在使用的任何和所有软件包 – Dason

+1

尝试使用't'吗? – blmoore

+0

图书馆提供,与t是完全相同的输出,我在hdfview中是非常奇怪的。 –

回答

4

rhdf5 docs

请注意,阵列显示为转置矩阵具有C-程序(h5dump或HDFView)打开时,它 。这是因为C上最快的变化维数是最后一个,但在R上它是第一个(如在Fortran中)。