如何将Stata生成的矩阵转换为对角矩阵？现有代码失效求助

你的核心需求是把生成的x1到x10变量转换成对角矩阵形式——也就是第i列（x_i）仅在第i行保留原始值，其余行全部设为0。先看看你原来代码的问题：

1. 类型不匹配错误：_n是数值型的行号，而"i'"是把循环变量i转成了字符串，用~=`比较数值和字符串会导致逻辑判断失效。
2. 循环逻辑冗余：你嵌套了foreach遍历所有变量，但其实每个x变量只需要针对自己的行号做处理，不需要遍历所有列。

修正后的代码方案一（基于你的原始思路调整） #

// 先生成原始变量
forv i = 1/10 {
    set obs 10
    gen x`i' = (_n-1)*10+`i'
}

// 转换为对角矩阵形式
forv i = 1/10 {
    // 对x`i'变量，仅保留第i行的值，其余行设为0
    replace x`i' = 0 if _n != `i'
}

更简洁的合并写法 #

如果想一步到位生成对角矩阵变量，可以把生成和转换合并成一个循环：

clear
set obs 10
forv i = 1/10 {
    gen x`i' = cond(_n == `i', (_n-1)*10+`i', 0)
}

这里用cond()函数直接判断：如果当前行号等于i，就生成原始值，否则设为0，代码更紧凑。

验证结果 #

运行完代码后，你可以用list x1-x10查看结果，会看到每一列只有对应行的位置有非0值，完全符合对角矩阵的要求。

内容的提问来源于stack exchange，提问作者user9539005