在MATLAB中创建火山图时出错

我是MATLAB的新手，我的第一个任务是创建火山图。我一直在使用它来了解它并开始

我试着用虚拟值运行它-

a=[1 2 3]
b=[4.6 2.7 4.5]
c=[0.05 0.33 0.45]

然后我跑了——

SigStructure = mavolcanoplot(a, b, c)

我的理解是，

a

表示条件1的基因表达值，

b

表示条件2，

c

是

a

和

b

中3个数据点的p值列表

但是运行此代码时会出现错误-

Index exceeds matrix dimensions.

Error in mavolcanoplot (line 127)
appdata.effect = X(paramStruct.goodVals) - Y(paramStruct.goodVals);

Error in volc (line 4)
SigStructure = mavolcanoplot(a, b, c)

---

有人能解释一下我的错误吗？

您遇到了一个问题，因为您使用的是行向量

在

mavolcanoplot

函数中（您可以通过在命令窗口中键入

edit mavolcanoplot

查看文件），有一个用于检查输入的本地函数，称为

check\u inputdata

您的数据通过了所有验证检查，然后遇到以下部分：

% Here, 'X' and 'Y' are the local names for your inputs 'a' and 'b'
% Below code is directly from mavolcanoplot.m:

% Handle the matrix input. Use its mean values per row
if size(X, 2) > 1
    X = mean(X,2);
end    
if size(Y, 2) > 1
    Y = mean(Y,2);
end

这会将您的输入减少到平均值。稍后在主函数（第127行）中，您会遇到如上所述的错误，其中

paramStruct.goodVals

是一个3元素数组，它现在正在尝试索引一个1元素数组，但失败了

---

这基本上是调试和阅读文档的一节课，其中

DataX，DataY：如果一个[…]矩阵，每一行是一个基因，每一列是一个样本，为每个基因计算一个平均表达值

您应该使用这些等效方法之一创建列向量输入

a=[1 2 3].';       % Using transpose (.') to create a column vector from a row vector
b=[4.6; 2.7; 4.5]; % Creating a column vector using the semi-colon operator to end each row
c=[0.05
   0.33
   0.45];          % Using actual code layout to create a column vector 

---

我不确定它是否有用，但文档中说c应该是一个列向量。