simsem程序包下代码运行结果求助

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

      大家好，我今天看了一篇“介绍如何对PLS最小样本量和统计检定力计算”的文章，Aguirre-Urreta, M., & Rönkkö, M. (2015). Sample size determination and statistical power analysis in PLS using r: An annotated tutorial. Communications of the Association for Information Systems, 36(1), 33-51.
      在该文中，作者提供了一串代码用以实现上述目的，代码如下：

library(matrixpls)

model<-"! regressions
A=~0.7*x1
A=~0.7*x2
A=~0.7*x3
B=~0.7*x4
B=~0.7*x5
B=~0.8*x6
B=~0.8*x7
C=~0.6*x8
C=~0.6*x9
C=~0.6*x10
C=~0.8*x11
C=~0.8*x12
D=~0.8*x13
D=~0.8*x14
D=~0.8*x15
D ~ 0.3*A
C ~ 0.1*B
C ~ 0.5*A
D ~ 0.3*C
! residuals, variances and covariances
A ~~ 1.0*A
B ~~ 1.0*B
C ~~ 0.71*C
D ~~ 0.725*D
B ~~ 0.3*A
x1 ~~ 0.51*x1
x2 ~~ 0.51*x2
x3 ~~ 0.51*x3
x4 ~~ 0.51*x4
x5 ~~ 0.51*x5
x6 ~~ 0.36*x6
x7 ~~ 0.36*x7
x8 ~~ 0.64*x8
x9 ~~ 0.64*x9
x10 ~~ 0.64*x10
x11 ~~ 0.36*x11
x12 ~~ 0.36*x12
x13 ~~ 0.36*x13
x14 ~~ 0.36*x14
x15 ~~ 0.36*x15"

# Normal data

output <- matrixpls.sim(1000, model, n=100, multicore = TRUE, completeRep = TRUE)

summary(output)

        然后，该文提供的运行结果是：


但是我却无法得到上述结果，我的控制台所显示的所有内容为：
> library(matrixpls)

载入程辑包：‘matrixpls’

The following object is masked from ‘package:stats’:

    loadings
>
> model<-"! regressions
+ A=~0.7*x1
+ A=~0.7*x2
+ A=~0.7*x3
+ B=~0.7*x4
+ B=~0.7*x5
+ B=~0.8*x6
+ B=~0.8*x7
+ C=~0.6*x8
+ C=~0.6*x9
+ C=~0.6*x10
+ C=~0.8*x11
+ C=~0.8*x12
+ D=~0.8*x13
+ D=~0.8*x14
+ D=~0.8*x15
+ D ~ 0.3*A
+ C ~ 0.1*B
+ C ~ 0.5*A
+ D ~ 0.3*C
+ ! residuals, variances and covariances
+ A ~~ 1.0*A
+ B ~~ 1.0*B
+ C ~~ 0.71*C
+ D ~~ 0.725*D
+ B ~~ 0.3*A
+ x1 ~~ 0.51*x1
+ x2 ~~ 0.51*x2
+ x3 ~~ 0.51*x3
+ x4 ~~ 0.51*x4
+ x5 ~~ 0.51*x5
+ x6 ~~ 0.36*x6
+ x7 ~~ 0.36*x7
+ x8 ~~ 0.64*x8
+ x9 ~~ 0.64*x9
+ x10 ~~ 0.64*x10
+ x11 ~~ 0.36*x11
+ x12 ~~ 0.36*x12
+ x13 ~~ 0.36*x13
+ x14 ~~ 0.36*x14
+ x15 ~~ 0.36*x15"
>
> # Normal data
>
> output <- matrixpls.sim(1000, model, n=100, multicore = TRUE, completeRep = TRUE)
Progress tracker is not available when 'multicore' is TRUE.
>
> summary(output)
   Length     Class      Mode
        1 SimResult        S4

由于我很少用R语言，无法找出其中的问题，各位高手可否告知一二？谢谢。

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：MSEM sim SEM mse IMS 程序

你把多核心这个参数去掉也许可以，我也出现了相同的问题，但是会运行很慢，也可能跟机器有关，给你截图看一下，毕竟做1000次蒙特卡罗模拟，每次抽取100个样本，这么大规模的计算一般机器承受不了，我这配置一般，运行了一个小时才跑了三分之一，关于你说的偏最小二乘法我能找到的资料不是很多，方便把原文发我一份吗，或者上传到人大经济论坛，这样可以更好的帮助你，看一下部分资料，现在只知道偏最小二乘回归≈多元线性回归分析＋典型相关分析＋主成分分析。

嗯，刚才说的没问题，可以啦，就是去掉multicore = TRUE,不过要运行很久很久，我足足等了两个小时。

victorchan0633 发表于 2015-10-2 22:05
你把多核心这个参数去掉也许可以，我也出现了相同的问题，但是会运行很慢，也可能跟机器有关，给你截图看一 ...

好的，我试试。原文说大概需要15分钟才能运算结束，由于我的是I7 ，八线程同时运算，10分钟之内就搞定了。那我现在试试单线程处理。原文附给你： Sample Size Determination and Statistical Power Analysis in PLS U.pdf (924.53 KB)

2015-10-3 08:23:53 上传

jinairenshi 发表于 2015-10-3 08:25
好的，我试试。原文说大概需要15分钟才能运算结束，由于我的是I7 ，八线程同时运算，10分钟之内就搞定了。 ...

嗯，好的，我下载了，看好后有新发现再交流