我总是得到一个20多岁左右的矩阵。它是对矩阵M的数和索引进行迭代,但是当我每次寻找一个只有唯一值的10x2的矩阵时,它就会写过它。
n = 20;
M = matrix(NA, ncol = 2, nrow = 10);
a = 1
b = 1
for (i in 1:n){
for (r in 1:nrow(M))
for (c in 1:ncol(M))
i -> M[r,c]
print(M)
}
M
在R中拟合重复测量,收敛性问题。我有以下拟合,它是许多数据集中的一个,它不收敛。我做了其他的套装。此数据集和模型在SAS中工作...我能得到一些关于如何在R中进行这项工作的指导吗?要查看的内容(矩阵、选项设置、r/splus的此主题参考...)。
fit.gls <- gls(resp~grpnum+grpnum/day,data=long, corr=cormat,na.action=na.omit)
Error in glsEstimate(object, control = control) :
computed "gls" fit is singular,
我正在编写R码来创建一个方阵。所以我的方法是:
分配一个大小正确的矩阵
循环遍历矩阵中的每个元素,并用适当的值填充它。
我的问题很简单:预先分配矩阵的最好方法是什么?到目前为止,我有两种方法:
> x <- matrix(data=NA,nrow=3,ncol=3)
> x
[,1] [,2] [,3]
[1,] NA NA NA
[2,] NA NA NA
[3,] NA NA NA
或
> x <- list()
> length(x) <- 3^2
> dim(x) <-
我想知道在Julia中是否有一个命令或包允许我们直接提取矩阵的下三角部分,不包括对角线。很明显,我可以为此调用R命令(比如gdata包的lowerTriangle ),但是我想知道朱莉娅是否有类似的东西。例如,假设我有矩阵
1.0 0.751 0.734
0.751 1.0 0.948
0.734 0.948 1.0
我不想创建一个下三角矩阵
NA NA NA
0.751 NA NA
0.734 0.948 NA
但将矩阵的较低部分提取为数组:0.751 0.734 0.948
如何选择R中不包含任何NA值的列?只要一列至少包含一个NA,我就想排除它。做这件事最好的方法是什么?我正试图使用sum(is.na(x))来实现这一目标,但并没有成功。
还有另一个R问题。是否可以使用命令排除包含所有相同值的列?例如,
column1 column2
row1 a b
row2 a c
row3 a c
我的目的是从我的矩阵中排除column1,所以最终的结果是:
column2
row1 b
row2 c
row3 c
我是新的R,下面有一个矩阵:
赚得
0 Name usd/day update_year usd/day update_year
1 John 52.0 2011 NA NA
2 Mary 44.0 2012 NA NA
3 Nicole 44.5 2000 est. 49.2 2010 est.
4 Cynthia 38.1 2001 est. 44.0 2011
我想清理R中的数据,只有3列-名称、美元/日和update_year,
我在R中执行分析,我想用表的内容填充矩阵的第一行。我遇到的问题是,表的内容根据数据的不同而变化,所以有时矩阵中出现的某些标识符不会出现在表中。
> random.evaluate
DNA LINE LTR/ERV1 LTR/ERVK LTR/ERVL LTR/ERVL-MaLR other SINE
[1,] NA NA NA NA NA NA NA NA
> y
DNA LINE LTR/ERVK LTR/ERVL LTR/ERVL-MaLR S
我尝试在R中按列绑定两个矩阵。矩阵具有不同的行大小。下面是一个例子:这是矩阵1:
M <- matrix(0, nr = 3, nc = 2)
和矩阵2:
N <- matrix(2, nr = 1, nc = 3)
输出需要如下所示:
0 0 2 2 2
0 0 NA NA NA
0 0 NA NA NA
尽管这个想法很简单。我没有设法实现它,也没有在互联网上找到解决方案。
我对r还是很陌生的(用来在Matlab中编程),我正在尝试使用并行软件包来加速一些计算。下面是一个例子,我试图计算一个矩阵的滚动标准差(按栏)使用动物园包,与并行和不并行的代码。然而,产出的形式却有所不同。
# load library
library('zoo')
library('parallel')
library('snow')
# Data
z <- matrix(runif(1000000,0,1),100,1000)
#This is what I want to calculate with timing
system.
我想为一个矩阵指定一个名称,而不是像R那样给矩阵的每个组件指定一个名称,例如:
a <- matrix(1:10, nrow=5)
names(a)
NULL
names(a) <- "this is a"
a
[,1] [,2]
[1,] 1 6
[2,] 2 7
[3,] 3 8
[4,] 4 9
[5,] 5 10
attr(,"names")
[1] "this is a" NA NA NA NA
因此,我正在研究R的矩阵,如下所示:
diff_0
SubPop0-1, SubPop1-1, SubPop2-1, SubPop3-1, SubPop4-1,
SubPop0-1, NA NA NA NA NA
SubPop1-1, 0.003403100 NA NA NA NA
SubPop2-1, 0.005481177 -0.002070277
我有两个矩阵,我想从其中提取对应单元格中所有唯一的值对。也就是说,我想将第一个矩阵中单元格1,1的值与第二个矩阵中的单元格1,1配对,然后将这个对列表缩小为只包含唯一对。我该怎么做?
下面是矩阵的示例版本
V1 V2 V3 V4 V5 V6
[1,] "w" "o" "r" "d" NA NA
[2,] "p" "o" "r" "c" NA NA
[3,] "v" "i" "n" NA N
我想用巧妙的方式制作一个三维的曲面图。这很好,但是x和y轴上的值是没有意义的。他们比他们应该的要高得多。我在图中使用了一个矩阵,x和y值是行和列名,z值(称为高程),当然是矩阵本身。实际上,似乎没有使用x和y数据的值。X和y值在矩阵中从0到6不等。这是我矩阵的一个小样本:
head(m.dune)
> 1.90 1.95 2 2.05 2.01 2.15
> 0 NA NA NA NA NA NA
> 0.05 NA NA NA NA NA NA
> 0.10 1.14 1.14 NA
下面的代码询问用户要创建的字数:
String sInt = JOptionPane.showInputDialog("Cuántas palabras vamos a capturar?");
int palabras = Integer.parseInt(sInt);
char[][] charMatrix = new char[palabras][];
其目的是将单词保存在一个矩阵中并打印出来,将每个字符保存在不同的矩阵维数中。
for (int i = 0; i<palabras; i++) {
String str
从R开始,我被用来通过使用索引将一个简单的值赋值给一个矩阵,例如:
> mat <- matrix(NA,2,3)
> mat
[,1] [,2] [,3]
[1,] NA NA NA
[2,] NA NA NA
> mat[1,2] <- 12
> mat
[,1] [,2] [,3]
[1,] NA 12 NA
[2,] NA NA NA
>
现在,我不得不用Python运行一些脚本,并希望收集2D矩阵中的值。矩阵应该有列和行名,我也想通过索引访问它们。到目前为止,我想到了两种方
在我的代码中,我从同一矩阵的每一列中减去一列矩阵。然后,我计算出有多少新列只有小于r的元素。我对矩阵的每一列都这样做。你可以在下面看到我的代码。我忽略了将值放入矩阵中的部分。
有什么方法可以提高这个代码的性能吗?我似乎想不出更快的办法
B = matrix(NA,(m),(window_step))
B_m_r = c(1:(window_step))
for (i in 1:(window_step)){
B_m_r[i] = sum(apply(abs(B[,-i]-B[,i]), 2,function(x) max(x) < r))
}
解决方案
B = matrix(N
我在R中有一个数据帧,如下所示:
a b c out1 out2
0 1 0 NA NA
1 0 1 NA NA
1 1 0 NA NA
DF <- data.frame(a = c(0,1,1), b= c(1,0,1), c = c(0,1,0), out1=NA, out2=NA)
我使用a、b和c的值作为创建两个矩阵的函数的输入。目前,如果我试图将一个矩阵赋给数据帧中的一个单元格(例如dataframe$ou
这是我的桌子
c1 c2 c3 c4 c5
r1 1 NA NA NA NA
r2 1 1 NA NA NA
r3 1 1 1 NA NA
r4 1 1 1 1 NA
r5 1 1 1 1 1
我想在每一列的底部移动NA,然后向上移动非NA数,而不需要在矩阵中做任何事情或安排非NA值(假设非NA值是唯一的),我只希望它们向上移动,使NA处于下面这样的位置:
c1 c2 c3 c4 c5
r1 1 1 1 1
我有一堆小矩阵,它们基本上是一个较大矩阵的子集,但值不同。我想从这些子矩阵中获取值,并在较大的矩阵中覆盖相应的值。例如,假设这是我的更大的矩阵:
AB-2000 AB-2600 AB-3500 AC-0100 AD-0100 AF-0200
AB-2000 6.5 NA -1.8 3.65 -17.96 -26.5
AB-2600 NA 7.18 NA NA NA NA
AB-3500 -1.79 NA 5.4 NA -4.63 NA
AC-0100