从Matlab导入数组到R

Question

我相信这是一个简单的问题，但我一直没有找到明显的解决方案。我有一系列由Matlab生成的模型输出数组文件（维度为180、360、12），我需要在R中打开。

我尝试使用R.matlab包，只使用readMat命令，结果得到一个列表对象。将该列表写入矩阵会导致内存分配错误。

我已经尝试过使用unlisting，但没有任何帮助。

如何打开这些Matlab矩阵文件并将其作为矩阵写入R？有什么想法吗？

---

到目前为止，用于读取Matlab的代码非常简单：

> data<-readMat("filename")
> typeof(data)
[1] "list"
> str(data)
List of 1
 $ pco2: num [1:180, 1:360, 1:12] NaN NaN NaN NaN NaN ...
 - attr(*, "header")=List of 3
 ..$ description: chr "MATLAB 5.0 MAT-file, Platform: GLNXA64, Created on: Thu Jul 26 10:36:42 2012 "
 ..$ version : chr "5"
 ..$ endian : chr "little

- Clive Trueman

---

请提供一个小的、可重现的示例，以便我们可以查看您的数据和错误。例如，在Matlab中，创建一个维度为3、4、5的矩阵。导入到R中并运行查询，如typeof(my_matrix)。 - Carl Witthoft

抱歉 - 我不是Matlab专家（这也许就是我遇到问题的原因！） - Clive Trueman

目前为止，Matlab的代码只是简单地读取数据<-readMat("filename")。typeof(data)返回"list"，str(data)返回1个列表 $pco2: num [1:180, 1:360, 1:12] NaN NaN NaN NaN NaN ...

• attr(*,"header")=List of 3 ..$ description: chr "MATLAB 5.0 MAT-file, Platform: GLNXA64, Created on: Thu Jul 26 10:36:42 2012 " ..$ version : chr "5" ..$ endian : chr "little"

- Clive Trueman

如果您编辑问题并将MATLAB代码作为代码块添加进去（通过换行后跟随四个空格来与文本区分开），那么这篇文章就会更易读和解释。 - Drew Steen

看起来所有的值都是“NaN”，这意味着原始数据可能是垃圾，或者存在翻译问题（可能是大小端问题）。 - Carl Witthoft

Answer 1

请仔细检查MATLAB中的文件，确保并非所有值都是NaN。只需在MATLAB中加载文件并检查其内容即可：

load file.mat

这里是我刚做的一个简单测试：

MATLAB

>> x = rand(3,4,2)

x(:,:,1) =
      0.75127      0.69908      0.54722      0.25751
       0.2551       0.8909      0.13862      0.84072
      0.50596      0.95929      0.14929      0.25428
x(:,:,2) =
      0.81428      0.34998      0.61604      0.83083
      0.24352       0.1966      0.47329      0.58526
      0.92926      0.25108      0.35166      0.54972

>> save file.mat x

R

R> library(R.matlab)
R> data <- readMat('file.mat')
R> str(data)
List of 1
 $ x: num [1:3, 1:4, 1:2] 0.751 0.255 0.506 0.699 0.891 ...
 - attr(*, "header")=List of 3
  ..$ description: chr "MATLAB 5.0 MAT-file, Platform: PCWIN, Created on: ..."
  ..$ version    : chr "5"
  ..$ endian     : chr "little"
R> data$x
, , 1

       [,1]   [,2]   [,3]   [,4]
[1,] 0.7513 0.6991 0.5472 0.2575
[2,] 0.2551 0.8909 0.1386 0.8407
[3,] 0.5060 0.9593 0.1493 0.2543

, , 2

       [,1]   [,2]   [,3]   [,4]
[1,] 0.8143 0.3500 0.6160 0.8308
[2,] 0.2435 0.1966 0.4733 0.5853
[3,] 0.9293 0.2511 0.3517 0.5497

- Amro