循环程序中，plot命令不出图，为何？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

这段估计DCC-GARCH的代码运行没有问题。但是2.1画时序图部分的代码，Rsudio软件不显示图，单独执行式没问题的，放进循环里就不行了。请教各位高手，这是什么原因？

1. setwd("D:/R/data/Ch23")
   
2. library(tseries);library(fBasics);library(FinTS);library(rmgarch)
   
3. #第一步 数据处理
   
4. for (i in 1:5){
   
5. sampA=read.csv(paste("A",i,".csv",sep = ""),header=T);
   
6. sampH=read.csv(paste("H",i,".csv",sep = ""),header = T);
   
7. indexA=xts::as.xts(sampA[,2],as.Date(sampA[,1]));
   
8. indexH=xts::as.xts(sampH[,2],as.Date(sampH[,1]));
   
9. A=diff(log(indexA[,1]))[-1,];H=diff(log(indexH[,1]))[-1,];
   
10. AH=na.omit(merge(A,H,join = 'inner'));
    
11. colnames(AH)=c(paste("A",i,sep = ""),paste("H",i,sep = ""));
    
12. print(head(AH));
    
13. #第二部 描述性统计
    
14.    # 2.1 画时序图
    
15. opar=par(no.readonly = T);
    
16. par(mfrow=c(2,1));
    
17. plot(AH[,1],main=paste("A",i,"日收益率序列图",sep = ""),col="blue")
    
18. plot(AH[,2],main=paste("H",i,"日收益率序列图",sep = ""),col="red")
    
19. par(opar);
    
20. par(mfrow=c(1,1));
    
21.    # 2.2 描述性统计
    
22. print(basicStats(AH));
    
23. print(jarque.bera.test(AH[,1]));print(jarque.bera.test(AH[,2]));
    
24. #第三步 实证分析
    
25.    # 3.1 平稳性检验
    
26. print(adf.test(AH[,1]));print(adf.test(AH[,2]));
    
27. print(pp.test(AH[,1]));print(pp.test(AH[,2]));
    
28.    # 3.2 ARCH效应检验
    
29. fit.arA=arima(AH[,1],order=c(1,0,0));
    
30. fit.arH=arima(AH[,2],order=c(1,0,0));
    
31. print(ArchTest(fit.arA$resid,lags = 4));print(ArchTest(fit.arH$resid,lags = 4));
    
32. print(ArchTest(fit.arA$resid,lags = 8));print(ArchTest(fit.arH$resid,lags = 8));
    
33. print(ArchTest(fit.arA$resid,lags = 12));print(ArchTest(fit.arH$resid,lags = 12));
    
34.    # 3.3 设定dcc-garch模型
    
35. meanEquation=list(armaOrder=c(1,0),include.mean=F,archpow=1);
    
36. varEquation=list(model="sGARCH",garchOrder=c(1,1));
    
37. distspec=c("mvt");
    
38. spec=ugarchspec(mean.model=meanEquation,variance.model=varEquation);
    
39. Nspec=multispec(replicate(2,spec));
    
40. myspec=dccspec(Nspec,VAR=F,robust=F,
    
41.                model = ("DCC"),
    
42.               external.regressors = NULL,
    
43.                distribution=distspec);
    
44. fitDcc=dccfit(data=AH,myspec,solver="solnp");
    
45. show(fitDcc);
    
46. plot(fitDcc,which=4)
    
47. }

复制代码

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：plot distribution DCC-GARCH replicate archTest

把16行，par(mfrow=c(2,1))；改为
par(mfcol=c(2,1), mar = c(2,2,2,2))
看看能出来吗

您好，还是出不来。把数据和代码放到附件里了。麻烦帮我指点一下

Ch23.rar (54.81 KB)

2022-5-8 13:54:35 上传

llb_321 发表于 2022-5-7 21:35
把16行，par(mfrow=c(2,1))；改为
par(mfcol=c(2,1), mar = c(2,2,2,2))
看看能出来吗

您好，还是出不来。把数据和代码放到附件里了。麻烦帮我指点一下

15~20行，替换一下，我试了可以在循环中出图：

1. op<-par(mfrow=c(1,1),mar=c(3,3,5,3))
   
2. par(mfrow=c(2,1),mar=c(3,3,3,3));
   
3. print(plot(AH[,1],main=paste("A",i,"日收益率序列图",sep = ""),col="blue"));
   
4. print(plot(AH[,2],main=paste("H",i,"日收益率序列图",sep = ""),col="red",which=all));
   
5. par(op);

复制代码

llb_321 发表于 2022-5-9 00:24
15~20行，替换一下，我试了可以在循环中出图：

厉害厉害！我也出来了，十分感谢哈。能讲解下设置的原理吗？

一般plot只默认画一张图，你使用par进行参数设置后，就是把画布分成了几份，每一份上都可以画图。

似乎是个bug，base包和ggplot2里面的画图，在循环里都要print才能输出。