问题 在R中将矩阵的对角线设置为零


我在这里尝试了一些提议的解决方案。但它不符合我的情况。 我有一个代码，这里：

a <- read.table("Whirr_127.csv", header=T, sep=",", row.names=1) # task assignment / people vs task
b <- read.table("Files_Whirr_127.csv", header=T, sep=",", row.names=1) #task vs files 
a
b

#calc cr , cr = ta * tf * transpose(ta)
cr <- as.matrix(a) %*% (as.matrix(b) %*% as.matrix(t(b)) %*% as.matrix(t(a)))
cr

#set value to 1, to initialize table
cr[cr>=1]<-1
cr

#identify diagonal matrix, set to zero
cr<-as.matrix(0,ncol=ncol(cr),nrow=nrow(cr))
cr<-diag(cr,x=0)

我想将对角线值设置为零。似乎最后两行中使用的代码对我的情况不起作用。

另外，我想使用文件名 a，并将其保存为AB_Files_Whirr_127.csv 我试着用

write.csv(cr,file = paste("CR_", a,".csv")

但是，我的目录中没有任何内容。

cr的样本输出：

               Adrian Cole Alison Wong Andrei Savu Bruno Dumon Edward J. Yoon Eugene Koontz Jakob Homan Kelvin Kakugawa Kirk True Lars George Soren Macbeth Stu Hood
Adrian Cole               0           0           0           0              0             0           0               0         0           0             0        0
Alison Wong               0           0           0           0              0             0           0               0         0           0             0        0
Andrei Savu               0           0           1           0              0             0           0               0         0           1             1        0
Bruno Dumon               0           0           0           0              0             0           0               0         0           0             0        0
Edward J. Yoon            0           0           0           0              0             0           0               0         0           0             0        0
Eugene Koontz             0           0           0           0              0             0           0               0         0           0             0        0
Jakob Homan               0           0           0           0              0             0           0               0         0           0             0        0
Kelvin Kakugawa           0           0           0           0              0             0           0               0         0           0             0        0
Kirk True                 0           0           0           0              0             0           0               0         0           0             0        0
Lars George               0           0           1           0              0             0           0               0         0           1             1        0
Soren Macbeth             0           0           1           0              0             0           0               0         0           1             1        0
Stu Hood                  0           0           0           0              0             0           0               0         0           0             0        0
Tibor Kiss                0           0           0           0              0             0           0               0         0           0             0        0
Tom White                 0           0           1           0              0             0           0               0         0           1             1        0
Unassigned                0           0           0           0              0             0           0               0         0           0             0        0
                Tibor Kiss Tom White Unassigned
Adrian Cole              0         0          0
Alison Wong              0         0          0
Andrei Savu              0         1          0
Bruno Dumon              0         0          0
Edward J. Yoon           0         0          0
Eugene Koontz            0         0          0
Jakob Homan              0         0          0
Kelvin Kakugawa          0         0          0
Kirk True                0         0          0
Lars George              0         1          0
Soren Macbeth            0         1          0
Stu Hood                 0         0          0
Tibor Kiss               0         0          0
Tom White                0         1          0
Unassigned               0         0          0

11837
2017-09-18 06:52


起源

答案:

a 不能在输出文件的名称中使用，因为它不是字符变量，而是数据框。

infile <- "Whirr_127.csv"
a <- read.table(infile, header=T, sep=",", row.names=1)
....
diag(cr) <- 0
write.csv(cr, file = paste0("CR_", infile, ".csv")

的语法 diag 对于新的R用户来说，这些行看起来很有趣，但它实际上只是调用赋值函数的另一种语法 diag<-，即 diag(x) <- 0 被解释为 diag<-(x, 0)。

更新：多个文件

如果要对多个配对文件重复上述操作，可以执行此操作。

a.files <- grep("^Whirr", dir(), value=TRUE)
b.files <- paste0("Files_", a.files)
for(i in length(a.files)){
    a <- read.table(a.files[i], ...)
    b <- read.table(b.files[i], ...)
    ...
    write.csv(cr, paste0("CR_", a.files[i], ".csv"))
}

12
2017-09-18 06:55


谢谢，对角线解决方案适合我，但文件没有出现在我的工作目录中... - user1676484
我可能在你读它的同时编辑了这篇文章。解决方案 paste0("CR_", infile, ".csv") 应该定义工作。如果不确定的话 paste0 命令返回一个字符标量。 - Backlin
它只是一个简单的诊断解决方案。傻我！无论如何，比方说，我有多个文件，这意味着我必须用filepath替换infile（file = paste0（“CR_”，infile，“。csv”）？它是否正确？ - user1676484
在每次迭代中执行循环并将infile设置为不同的文件。我怀疑这些文件成对出现 "Whirr_###.csv" 和匹配 "Files_Whirr_###.csv" 不过，是吗？ - Backlin
你是对的 - user1676484

答案:

a 不能在输出文件的名称中使用，因为它不是字符变量，而是数据框。

infile <- "Whirr_127.csv"
a <- read.table(infile, header=T, sep=",", row.names=1)
....
diag(cr) <- 0
write.csv(cr, file = paste0("CR_", infile, ".csv")

的语法 diag 对于新的R用户来说，这些行看起来很有趣，但它实际上只是调用赋值函数的另一种语法 diag<-，即 diag(x) <- 0 被解释为 diag<-(x, 0)。

更新：多个文件

如果要对多个配对文件重复上述操作，可以执行此操作。

a.files <- grep("^Whirr", dir(), value=TRUE)
b.files <- paste0("Files_", a.files)
for(i in length(a.files)){
    a <- read.table(a.files[i], ...)
    b <- read.table(b.files[i], ...)
    ...
    write.csv(cr, paste0("CR_", a.files[i], ".csv"))
}

12
2017-09-18 06:55


谢谢，对角线解决方案适合我，但文件没有出现在我的工作目录中... - user1676484
我可能在你读它的同时编辑了这篇文章。解决方案 paste0("CR_", infile, ".csv") 应该定义工作。如果不确定的话 paste0 命令返回一个字符标量。 - Backlin
它只是一个简单的诊断解决方案。傻我！无论如何，比方说，我有多个文件，这意味着我必须用filepath替换infile（file = paste0（“CR_”，infile，“。csv”）？它是否正确？ - user1676484
在每次迭代中执行循环并将infile设置为不同的文件。我怀疑这些文件成对出现 "Whirr_###.csv" 和匹配 "Files_Whirr_###.csv" 不过，是吗？ - Backlin
你是对的 - user1676484

热门问题

不使用eval / new函数的JavaScript模板库 当涉及内部类时，Java继承如何工作 .NET Windows服务的奇怪问题 在.ipa或.app下查找App ID 快速入门XSLT参考[关闭] 如何找出Android应用程序中未使用的资源 Ruby中并发的同步方法[重复] 将std :: chrono :: system_clock :: time_point转换为struct timeval并返回 Google Drive API V3（javascript）更新文件内容 Bootstrap 3.0 - 将元素推到底部 受密码保护的.NET ClickOnce部署？ 如何用postgresql安装wordpress coq Set或Type如何成为命题 硒滚动元素进入（中心）视图 在Spring Transaction JUnit测试中自动装配Hibernate会话的正确方法 Git的Dockerfile策略 如何在FOS_PICKFOLDER中使用IFileDialog，同时仍在对话框中显示文件名 在Firefox扩展中复制Google Chrome浏览器操作弹出效果 CakePHP找到MAX 芹菜 - 完成任务的召唤功能 从使用fmemopen创建的流中读取宽字符 .NET是否为每个程序集创建一个字符串实习池？ DefaultModelBinder不绑定嵌套模型 Navigator.MediaDevices.getUserMedia（）使用了哪些相机通信标准？ 选择命名空间名称时应该知道什么？ cout Swagger Codegen CLI Java客户端 - 如何正确使用它 一个很好的哈希函数用于采访整数，字符串？ Maven 3 ciManagement配置的目的是什么？ 如何通过语言文化获取代码页？
licensed under cc by-sa 3.0 with attribution.
隐私政策 往来