25
4
I have calculated temperature difference between two model experiments. The difference is a 2D matrix (lon×lat). I want to find the grid points where the temperature differences are statistically significant at 95% significance level based on the student's t-test using MATLAB.
我已经计算了两个模型实验之间的温差。不同之处在于2D矩阵(长×长)。我想找出温差在95%显着水平上具有统计学意义的网格点，这是基于学生使用MatLab进行的t检验。
Now, I used ttest function in MATLAB like in the following:
现在，我使用了MATLAB中的Ttest函数，如下所示：
h = ttest(temp_diff,[],'Alpha',0.01)
But, h is coming as a vector of 1×231 (note that temp_diff has dimensions of 152×231, so 231 corresponds to latitude dimension). But, it should come as a 2D matrix with dimension same as temp_diff, i.e, 152×231 so that I can give markers on each significant grid points in the contourf spatial map (latitude vs longitude). So, how to do that in MATLAB? Can anyone please guide me in this regard? That will be very helpful for me. Thank you for your time and consideration.
With regards,
Ankan
H=Ttest(temp_diff，[]，‘Alpha’，0.01)，但是，h是一个1×231的向量(请注意，temp_diff的维度是152×231，因此231对应于纬度维度)。但是，它应该是一个维度与temp_diff相同的2D矩阵，即152×231，这样我就可以在等高线空间地图(纬度与经度)中的每个重要网格点上给出标记。那么，如何在matlab中做到这一点呢？有没有人能在这方面给我指点一下？这对我很有帮助。感谢您抽出时间考虑我的问题。说到这里，安坎
You cannot apply a t-test to a single number. You need multiple measurements to decide if their mean is statistically different from 0 or not. This is why ttest works on columns, it does the statistics on the set of numbers in each column.
您不能对单个数字应用t检验。您需要多次测量才能确定它们的平均值是否在统计上不同于0。这就是Ttest在列上工作的原因，它对每列中的一组数字进行统计。
Thank you @CrisLuengo for your response. So, if I want to find the grid points where the temperature differences are statistically significant at 95% significance level based on the student's t-test using MATLAB, then should I take two temperature mean value from the two experiments (which are used to compute temperature differences) like in the following: [h,p]=ttest2(temperature1_mean,temperature2_mean)?
感谢您@CrisLuengo的回复。因此，如果我想找出网格点，其中温差在95%显着水平的学生使用MATLAB的t检验，那么我是否应该从两个实验的两个温度平均值(用于计算温差)如下：[h，p]=test2(温度1_平均值，温度2_平均值)？
You cannot do what you want to do. The t-test uses multiple measurements to estimate the error in the mean. You have only one measurement per grid point, so you can't do those statistics. What you can do is estimate the error in your measurements in another way, maybe by knowing the device with which you've done the measurements. Or you can take multiple measurements.
你不能做你想做的事。T检验使用多个测量值来估计平均值中的误差。每个网格点只有一个测量值，因此不能进行这些统计。你能做的是用另一种方式估计你的测量误差，也许是通过知道你用来进行测量的设备。或者你也可以进行多次测量。
Thank you @CrisLuengo for your response.
感谢您发送编修。
25
4
0
问题来自 Julia Discourse 我正在使用 Julia 1.2。这是我的测试: a = rand(1000, 1000) b = adjoint(a) c = copy(b) @btime
给定两个 m x n 矩阵 A 和 B，其元素属于集合 S。 问题:A 的行和列可以置换为 B 吗？ 解决这个问题的算法的复杂度是多少？ 行列式有部分帮助(当 m=n 时):必要条件是 det(A)
这是共同计算平移和旋转的正确方法，还是有更好的方法？目前我的代码先翻译然后旋转，这会造成问题吗？ 代码 from math import cos, sin, radians def trig(angl
如何在 core.matrix 中将矩阵和向量元素相乘？ 我正在尝试找到与以下 Octave 音程代码等效的 core.matrix: A = [1 2 3 4; 5 6 7 8] B = [2; 3
我是 Scilab 的新用户(另见 here)。 我定义了一个简单的分段函数，并在使用该函数 ( "Warning adding a matrix with the empty matrix will
我有一个像这样的 RDD: (A,AA,1) (A,BB,0) (A,CC,0) (B,AA,2) (B,BB,1) (B,CC,4) 我想将其转换为以下 RRD: ([1,0,0],[2,1,4])
我使用的矩阵如下 (require '[clojure.core.matrix :as ccm]) (def M (ccm/matrix [[1 2] [3 4]])) (ccm/mset! M 0
我有一个矩阵类，它有一组函数，其中一个是矩阵运算符++(); 构造函数: Matrix(int num_rows,int num_col,int initialization,double initi
我有一个矩阵如下； 1 2 3 4 5 1 0 1 1 0 0 2 0 0 1 1 0 3 1 0 0 0 1 4 0 0 1 0 0
我已经部署了为家庭服务器(synapse)运行的单个实例，并附加了多个域作为 example.com 和 example1.com。我想创建像 [email protected] 这样的用户和 [em
我有 200 个向量；每一个的长度都是 10000。 我想填充一个矩阵，使每一行代表一个向量。 最佳答案 如果你的向量已经存储在一个数组中，那么你可以在这里使用 vcat( ): A = [rand(
如何向现有矩阵添加行或列？我正在尝试添加一个偏差项(一列)作为矩阵的第一行。在 Octave 中我可以这样做: M = [ones(size(M, 1), 1), M]; 最佳答案 您可以使用 joi
我正在使用 GNU GSL 进行一些矩阵计算。我正在尝试将矩阵 B 与矩阵 A 的逆矩阵相乘。 现在我注意到 GSL 的 BLAS 部分有一个函数可以做到这一点，但前提是 A 是三角形。这有什么具体原
我想计算如下:Matrix * Matrix Matrix有大约 6M*3 个元素，如何转换 Matrix至 Matrix这样我就能得到 Matrix结果。 最佳答案 您可以使用 Map 函数将 do
我只是没有看到我的错误。关于此错误消息的问题太多了，答案要么不适用，要么我只是看不到它们适用。也许应该改进错误消息？ Matrix a = Matrix(3, 4); // fill a with v
在android.opengl.Matrix类中有两种旋转矩阵的方法，它们是: static void rotateM (float[] m, int mOffset, float a, float
我正在使用 C++ 进行编码，并且使用的是 FEniCS fenics/2016.1.0。我的部分代码是 Matrix A; Vector f; std::vector> dirichlet_matr
JAMA(用于矩阵计算的 java 库)中的 JAMA:Matrix.times() 与 Matrix.arrayTimes() 有什么区别 如果我有一个d维度 vector x和一个k维度 vect
我试图做的是简单地将 cublasDgemm(矩阵-矩阵乘法)应用于多个具有“双”(8 字节)类型元素的矩阵，所有这些元素都具有一个非常大的维度。在我的例子中，矩阵的大小是 12755046 x 46
我正在尝试使用 android Matrix 对象旋转给定的位图。 我想将它发送到我的服务器，我正在使用 Android API8。 我应该使用 Matrix.setRotate 还是 Matrix.
个人中心
文章发布
退出
我是一名优秀的程序员,十分优秀！