R语言矩阵计算 - 经管之家

8567

收藏 2013-11-25

大家好，我手里有一个矩阵，

想计算两个人是否毕业于同一个学校。
最后的结果应该是

用户1	用户2	是否毕业于同一学校
A	B	否
A	C	是
B	C	否

用R该怎么写？多谢~

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

2013-11-25 17:21:44

这是mysql的专长，用R做不太合适，速度慢，内存消耗高

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2013-11-25 18:08:04

nuomin 发表于 2013-11-25 17:21
这是mysql的专长，用R做不太合适，速度慢，内存消耗高

哦，多谢。我写的循环
n <- length(y <- matrix(x[,1]))
y <- matrix(0,nrow=3,ncol=3)
for (i in (1:(n-1))) {
for (j in (1:(n-i))){
relation <- length(intersect(x[i,2],x[i+j,2]))/length(union(x[i,2],x[i+j,2]))
#print (x[n,2])
print(x[i,1])
print(x[i+j,1])
print(relation)
}
}

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

qoiqpwqr

2013-11-25 18:48:34

复制代码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2013-11-26 09:39:55

qoiqpwqr 发表于 2013-11-25 18:48

谢谢，写得太专业了，感觉我的思维还没有转到R上面来

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2013-11-26 23:22:15

感觉可以不用循环
尝试了下，过程比较狗血
其想法是把 school 数字化factor（且对于ID)，A和T(A),相乘变矩阵，矩阵的列可以理解为school或者ID，行也是school或者ID，如果同校数值相同，则sqrt(此矩阵)为school的数值（某整数）。此时取矩阵的对角线下部分整数即可，并得到坐标，从而找到对应的ID。
而在sqrt()时，2*2,3*3,4*4，，，同校，但要避免1*4,2*8这类的数sqrt()也为整数，但并不同校，可采取数值-0.1的变通下，sqrt(..-0.1)+0.1...

复制代码

命名、细节有些糙有补充的吗？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

点击查看更多内容…

2013-11-27 21:01:44

有R味道的代码要没有循环，以向量运算为主。我也来贴一段代码，比以上的简洁一些，简单一些。

复制代码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2013-11-27 22:50:38

nuomin 发表于 2013-11-27 21:01
有R味道的代码要没有循环，以向量运算为主。我也来贴一段代码，比以上的简洁一些，简单一些。

原来有 combn函数

那就更味儿一点的

[titter]

复制代码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2013-11-28 08:45:27

jmpamao 发表于 2013-11-27 22:50
原来有 combn函数

那就更味儿一点的

简洁、味儿浓，哈哈

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2013-11-28 10:13:10

nuomin 发表于 2013-11-28 08:45
简洁、味儿浓，哈哈

越来越看不懂了。。。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2013-11-28 10:13:14

nuomin 发表于 2013-11-28 08:45
简洁、味儿浓，哈哈

越来越看不懂了。。。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2013-11-28 10:33:31

nuomin 发表于 2013-11-27 21:01
有R味道的代码要没有循环，以向量运算为主。我也来贴一段代码，比以上的简洁一些，简单一些。

与matlab一样，R也是以向量为主。。。我之前用过FORTRAN，也是主要以写循环为主，我都有点思维定势了

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2013-11-28 21:20:31

sosmall1181 发表于 2013-11-28 10:33
与matlab一样，R也是以向量为主。。。我之前用过FORTRAN，也是主要以写循环为主，我都有点思维定势了

在我看来能用像fortran、lisp这样的古老语言的都是牛人

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2013-11-29 09:30:38

nuomin 发表于 2013-11-28 21:20
在我看来能用像fortran、lisp这样的古老语言的都是牛人

我都用简单的加加减减，复杂的也是不会，都不精啊。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

hayes

2013-11-29 10:46:17

jmpamao 发表于 2013-11-27 22:50
原来有 combn函数

那就更味儿一点的

高人啊~~~~~膜拜

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2013-11-29 14:03:35

jmpamao 发表于 2013-11-27 22:50
原来有 combn函数

那就更味儿一点的

如果一个人对应两个学校，用intersect也不可以，该用什么函数判断是否有交集呢？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2013-11-29 14:57:34

sosmall1181 发表于 2013-11-29 14:03
如果一个人对应两个学校，用intersect也不可以，该用什么函数判断是否有交集呢？

复制代码

是这种吗？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2013-11-29 15:12:40

jmpamao 发表于 2013-11-29 14:57
是这种吗？

膜拜。。。作为菜鸟，需要更加用功，谢谢！

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2013-11-29 15:22:32

sosmall1181 发表于 2013-11-29 15:12
膜拜。。。作为菜鸟，需要更加用功，谢谢！

一起努力吧我是跟那些高手后面跑跑龙套的
只知道些前期的数据清洗。后面的要向在座的学习

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2013-11-29 15:56:46

jmpamao 发表于 2013-11-29 15:22
一起努力吧我是跟那些高手后面跑跑龙套的
只知道些前期的数据清洗。后面的要向在座的学习

一起努力。但是如果矩阵比较大的话，用combn函数可能效率就比较低

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2013-11-29 16:56:37

qoiqpwqr 发表于 2013-11-25 18:48

请教一下，能用apply函数替代那个循环么，不用combn函数。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2013-11-29 21:07:54

sosmall1181 发表于 2013-11-29 16:56
请教一下，能用apply函数替代那个循环么，不用combn函数。

apply不是用那种格式的，最少我没见过。我把之前的运用 matrix的方法改了下。觉得ok

复制代码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2013-11-30 11:34:39

不用combn函数的，纯纯的矩阵

复制代码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2013-11-30 13:48:15

nuomin 发表于 2013-11-30 11:34
不用combn函数的，纯纯的矩阵

复制代码

这部分你是怎么想到的？怎么推出来的

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

jgchen1966

2013-12-1 01:35:48

其实，这仅是一个查询任务，可用dlply() 非常简单的现实：
dat<- matrix(c("A","B","C","rmb","zkj","rmb"),ncol=2)
dat<-as.data.frame(dat,stringsAsFactors=FALSE)
names(dat)<-c("ID","School")
datnew<-dlply(dat,.(School),f=ID)
查询结果：

> datnew$rmb
  ID School
1  A rmb
2  C rmb
> datnew$zkj
  ID School
1  B zkj
>

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2013-12-1 08:35:46

jmpamao 发表于 2013-11-30 13:48
这部分你是怎么想到的？怎么推出来的

看c_matrix矩阵总结的

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2013-12-2 10:28:40

nuomin 发表于 2013-12-1 08:35
看c_matrix矩阵总结的

不知道在R中能不能实现一个功能，就是如果A有多个学校，存在不同的行里，能不能将这些记录合并为一行，学校中间用，或其他符号区隔

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2013-12-2 13:06:42

sosmall1181 发表于 2013-12-2 10:28
不知道在R中能不能实现一个功能，就是如果A有多个学校，存在不同的行里，能不能将这些记录合并为一行，学 ...

复制代码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2013-12-2 14:42:38

jmpamao 发表于 2013-12-2 13:06

谢谢。
我现在想做的是将相同的学校找出来，数据last如下：

apply(last,1,function(i){intersect(unlist(strsplit(as.character(i[2]),split=",")),unlist(strsplit(as.character(i[6]),split=",")))})
但是学校出来的格式不对，

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝