如果有谁遇到过要算比较长时间的分析也行,最好用stata自带的例子数据

找到了,用2个大矩阵相乘
set obs 1000
g x1=(runiform())
. set matsize 10000

Current memory allocation

                    current                                 memory usage
    settable          value     description                 (1M = 1024k)
    --------------------------------------------------------------------
    set maxvar         5000     max. variables allowed           1.947M
    set memory           50M    max. data space                 50.000M
    set matsize       10000     max. RHS vars in models        763.779M
                                                            -----------
                                                               815.726M

. mkmat x, mat(m)

. c
unrecognized command:  c
r(199);

. matrix n=m'

. matrix c=m*n

di c(current_time)
set obs 1000
g x1=(runiform())
di c(current_time)
mkmat x, mat(m)
matrix n=m'
matrix c=m*n
di c(current_time)


. clear

. do "D:\app\etime.do"

. di c(current_time)
16:18:43

. set obs 1000
obs was 0, now 1000

. g x1=(runiform())

. di c(current_time)
16:18:43

. mkmat x, mat(m)

. matrix n=m'

. matrix c=m*n

. di c(current_time)
16:21:01

.
.
.
end of do-file

测试结果 10.1 32位mp 比楼上11.0 64位mp快多了
. set matsize 10000

Current memory allocation

                    current                                 memory usage
    settable          value     description                 (1M = 1024k)
    --------------------------------------------------------------------
    set maxvar         5000     max. variables allowed           1.909M
    set memory           10M    max. data space                 10.000M
    set matsize       10000     max. RHS vars in models        763.779M
                                                            -----------
                                                               775.688M

. do "D:\app\etime.do"

. di c(current_time)
16:24:39

. set obs 1000
obs was 1000, now 1000

. g x1=(runiform())
x1 already defined
r(110);

end of do-file

r(110);

. clear

. do "D:\app\etime.do"

. di c(current_time)
16:24:52

. set obs 1000
obs was 0, now 1000

. g x1=(runiform())

. di c(current_time)
16:24:52

. mkmat x, mat(m)

. matrix n=m'

. matrix c=m*n

. di c(current_time)
16:24:53

.
.
.
end of do-file

求逆矩阵不成功

. mat v=inv(c)
matrix has missing values