问题标签 [sas-iml]

希望有一个简单的答案。我正在做一项模拟研究，我需要在大约一千次重复中从均匀分布 U(25,200) 中抽取随机数量的个体 N。一次复制的代码如下所示：

我在 DATA 步骤之外创建了这两个宏变量，因为我需要在随后的 DATA 步骤和 SAS 和 IML 中的 DO 循环中重复调用 N 变量。

问题是每次我在复制中调用 N 时，它都会重新采样 U，这必然会修改 N。因此，N 在复制中不会保持不变。这个问题显示在下面的代码中，我首先将 N 创建为一个变量（在个体之间是恒定的），并使用 DO 循环为每个个体创建 X 的样本预测值。注意N中的值与个体总数不一样，这也是个问题。

我猜我需要做的是在整个复制过程中以某种方式保持 U 不变，然后允许对复制 2 重新采样，依此类推。通过保持 U 不变，N 将必然保持不变。

有没有办法使用宏变量来做到这一点？

在 SAS IML 中，如果两个矩阵包含缺失值，则不能相乘以获得两个矩阵的乘积，这对吗？

SAS 中没有解决这个问题的方法？？？

在 SAS 9.3 中，我们可以使用函数轻松计算 lat、lon 之间的欧几里得距离矩阵distance

我在 SAS 9.2 中找不到此功能我如何在 9.2 中执行此操作？

谢谢你。

我是 SAS 新手，我有这个例子：

如何仅使用 Z=1 和仅 Z=2 创建新数据？

谢谢你。

我正在尝试以下代码：

所以我希望在完成 do 循环后获得 20 个数据集。SAS有可能吗？我可以使用宏来做到这一点，但我不喜欢在内部使用宏PROC IML！

提前致谢。

我有一个名为 m1 的 6*6 矩阵，我想在 SAS 中使用 Do Loop 创建矩阵，使得 m2=m1*m1; m3=m2*m1；m4=m3*m1 ... mi=m(i-1)*m1。

这是我写的：

而且它不起作用，因为 cat(m,1) 可能不正确。我该如何使用 Do Loop 呢？非常感谢您的时间和帮助！

如何在 SAS 中使用排列进行方差分析？这样做的目的是缩短运行时间并保存输出而不是打印出来。我曾经使用 iml 进行了配对 t 检验和排列。进展顺利。但我有点担心 iml 是否也可以在 ANOVA 中进行排列？我还尝试了 multtest 来生成排列。multtest 排列很长，而 iml 很长。哪个效率更高？谢谢！

我正在将一个包从 R 翻译成 IML，完成后将在网上免费 :)。我从一个大矩阵的分解中得到了不同的结果，当你看它们时，这两个结果看起来是一样的，但是，例如，如果我取 U 的前 2 列并做 U'*U，我的 2x2 矩阵将完全不同（ U_11 = 1.1e-17 和 U_11 =1.4e-17)。差异非常小（3e-18），这让我认为这可能与每个软件使用的小数位数有关，SAS IML 和 R。有人对这个主题有更多了解吗？我该如何测试呢？谢谢你。

不久，我正在将一个 R 包翻译成 IML 语言，我完全为 R 和 IML 之间的 SVD 分解结果而苦苦挣扎。

代码：

SAS/IML 代码：

V 是矩阵的两边，包含来自 SVD 分解的右奇异向量，但是在这里，请看一下 R 的这些结果：

这些来自 SAS/IML：

如您所见，这些值相同，但它们的符号有时相同，有时不同。我找不到发生这种情况的任何原因，这完全让我发疯。有人可以给我一个解决方案，或者至少是一个提示或任何东西。任何帮助将不胜感激，即使它不能解决这个问题。

非常感谢。

编辑：显然，发布的结果只是总矩阵的一部分。

我需要为使用 SAS相关的两个beta 分布变量生成随机值。感兴趣的两个变量的特征如下：

X1有mean = 0.896和variance = 0.001。

X2有mean = 0.206和variance = 0.004。

对于X1和X2，p = 0.5，其中p是相关系数。

使用 SAS，我了解如何使用函数生成指定 beta 分布的随机数X = RAND('BETA', a, b)，其中a和b是变量X的两个形状参数，可以从均值和方差计算得出。但是，我想同时为两者生成值，X1同时X2指定它们在p = 0.5 时相关。