代码运行后显示十二个图，要求将12个图合并成三个图，每个图有四个小图，求大神

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

library(sensitivity)



#####------------------------- Plot Morris test results -------------------------#####

output_names  <- c("opt-out","imitation","optimization","repetition") # Corresponding to c(1,2,3,4)


####---------- VII. Analysis of test results (postprocessing)  ----------####

paraNames <- c("ref_income", "tol_income", "tol_activity", "lambda",
               "alpha_plus", "alpha_minus","phi_plus", "phi_minus")

### plots for mu* against mu/sigma
for (i in (1:length(sim_results_morris_org[,1])))
{
  # add simulation results (as vector) to morris object
  tell(mo, sim_results_morris_org[i,])
  mu <- apply(mo$ee, 2, mean)
  mu_star <- apply(mo$ee, 2, function(x) mean(abs(x)))
  sigma <- apply(mo$ee, 2, sd)

  ## mu_star against mu
  # ask user for pressing <Enter> for entering the plot
  plot(mu_star, mu)
  title(main=output_names)
  text(mu_star, mu, labels=paraNames, cex= 0.8, pos=3)

  # mu_star against sigma
  plot(mu_star, sigma)
  title(main=output_names)
  text(mu_star, sigma, labels=paraNames, cex= 0.8, pos=3)
  # mu against sigma
  plot(mu, sigma)
  title(main=output_names)
  text(mu, sigma, labels=paraNames, cex= 0.8, pos=3)
}

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：十二个 Sensitivity Library results output

plot前面加一个op=par(mfrow=c(2,2))

zhou1_20 发表于 2019-6-27 09:50
plot前面加一个op=par(mfrow=c(2,2))

library(sensitivity)



#####------------------------- Plot Morris test results -------------------------#####

output_names  <- c("opt-out","imitation","optimization","repetition") # Corresponding to c(1,2,3,4)


####---------- VII. Analysis of test results (postprocessing)  ----------####

paraNames <- c("ref_income", "tol_income", "tol_activity", "lambda",
               "alpha_plus", "alpha_minus","phi_plus", "phi_minus")

### plots for mu* against mu/sigma
for (i in (1:length(sim_results_morris_org[,1])))
{
  # add simulation results (as vector) to morris object
  tell(mo, sim_results_morris_org[i,])
  mu <- apply(mo$ee, 2, mean)
  mu_star <- apply(mo$ee, 2, function(x) mean(abs(x)))
  sigma <- apply(mo$ee, 2, sd)
  
  ## mu_star against mu
  # ask user for pressing <Enter> for entering the plot
  op=par(mfrow=c(2,2))
  plot(mu_star, mu)
  title(main=output_names)
  text(mu_star, mu, labels=paraNames, cex= 0.8, pos=3)
  
  # mu_star against sigma
  op=par(mfrow=c(2,2))
  plot(mu_star, sigma)
  title(main=output_names)
  text(mu_star, sigma, labels=paraNames, cex= 0.8, pos=3)
  # mu against sigma
  op=par(mfrow=c(2,2))
  plot(mu, sigma)
  title(main=output_names)
  text(mu, sigma, labels=paraNames, cex= 0.8, pos=3)

是这样吗，怎么显示不出来图形结果呀

zhou1_20 发表于 2019-6-27 09:50
plot前面加一个op=par(mfrow=c(2,2))

library(sensitivity)



#####------------------------- Plot Morris test results -------------------------#####

output_names  <- c("opt-out","imitation","optimization","repetition") # Corresponding to c(1,2,3,4)


####---------- VII. Analysis of test results (postprocessing)  ----------####

paraNames <- c("ref_income", "tol_income", "tol_activity", "lambda",
               "alpha_plus", "alpha_minus","phi_plus", "phi_minus")

### plots for mu* against mu/sigma
for (i in (1:length(sim_results_morris_org[,1])))
{
  # add simulation results (as vector) to morris object
  tell(mo, sim_results_morris_org[i,])
  mu <- apply(mo$ee, 2, mean)
  mu_star <- apply(mo$ee, 2, function(x) mean(abs(x)))
  sigma <- apply(mo$ee, 2, sd)
  
  ## mu_star against mu
  # ask user for pressing <Enter> for entering the plot
  op=par(mfrow=c(2,2))
  plot(mu_star, mu)
  title(main=output_names)
  text(mu_star, mu, labels=paraNames, cex= 0.8, pos=3)
  
  # mu_star against sigma
  op=par(mfrow=c(2,2))
  plot(mu_star, sigma)
  title(main=output_names)
  text(mu_star, sigma, labels=paraNames, cex= 0.8, pos=3)
  # mu against sigma
  op=par(mfrow=c(2,2))
  plot(mu, sigma)
  title(main=output_names)
  text(mu, sigma, labels=paraNames, cex= 0.8, pos=3)

是这样吗，怎么显示不出来图形结果呀

加到for前面看下op=par(mfrow=c(2,2))

HHLDPPH 发表于 2019-6-27 10:29
library(sensitivity)

显示出来的不是2*2的形式

zhou1_20 发表于 2019-6-27 12:04
加到for前面看下op=par(mfrow=c(2,2))

嗯嗯，我按照你说的，已经做出来三张2*2的图，但是横纵坐标对应的有问题,是混乱的。按说是有三组，每个类型对应一张图：1.mu_star/mu   2.mu_star/sigma  3. mu/sigma  不知道问题出在哪

我做出来就是这样的