关于一个算法的实现

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

问题有些复杂。有三个变量 a b c
a           b             c
1          m            c1
1          n             c2
1          s             c3
2          m            c4
2          n             c5
2          s             c6
2          t             c7
3         m           c8
3         n            c9
4
。
。
5


。。。。。。

c列是一些数值。想完成的是：
对于a列中任何两个观测值，（比如，1&2, 2&3, 1&3......）需要计算b列值相同的c列的乘积之和（拗口。。。）
比如对于a列中的1&2
需要计算：c1*c4+c2*c5+c3*c6 （如果没有匹配的就忽略，如c7,而且数据能够保证如果存在的话对应是一一的）
最后还想把这个乘积之和放到一个矩阵的（1，2）这个位置上。

弄了很久搞不定，望高手指点，谢谢！！

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：高手指点 存在的 观测值 而且

死办法，循环

如何写这个循环呢 还是嵌套的循环
还有如何搜索b列match的记录呢
俺是初级菜鸟 谢谢了

doreenjiao 发表于 2010-9-8 17:18
需要计算：c1*c4+c2*c5+c3*c6 （如果没有匹配的就忽略，如c7,而且数据能够保证如果存在的话对应是一一的）

*没明白红线部分,姑且认为不作处理， 如果b为string

*试试

reshape wide c,i(a) j(b) string
sum a
loca lm=r(max)
matr x=J(`lm',`lm',0)
forvalue i=1/`lm' {
forvalue j=1/`lm' {
if `i'~=`j'{  //对于a列中任何两个观测值，如果包含相同，则去掉该条件
   foreach vnm of varlist c* {
    if `vnm'[`i']~=. & `vnm'[`j']~=.{
     mat x[`i',`j']=`vnm'[`i']*`vnm'[`j']+x[`i',`j']
    }
   }
  }
}
}
matlist x
*等于0视为不合条件
*有一种情况没有考虑：真正的积和等于0

1# doreenjiao
* 对积尔和为0的克服
tempfile d1
reshape wide c,i(a) j(b) string
save `d1',replace
renpfix "" t
levelsof ta,loc(lvl)
expand _N,
gsort ta
egen a=fill(`lvl' `lvl' )
merg m:1 a using `d1'
foreach v of varlist c*{
gen j`v'=t`v'*`v'
}
egen jsum=rowtotal(j*)
keep a ta jsum
reshape wide jsum, i(a) j(ta)
mkmat jsum*,mat(x)
matlist x

Try this:

preserve
encode b, gen(d)
summ d
local maxd= r(max)
drop b

reshape wide c, i(a) j(d)
foreach x of varlist * {
        replace `x' = 0 if `x' == .
}

gen value = 0
forvalue i = 1 / `maxd' {
        replace value = value + c`i' * c`i'[_n-1] if _n>1
}

mkmat a value, matrix(data)
matrix list data

restore

多谢! 我来慢慢学习一下.