plot.new has not been called yet报错

问题遇到的现象和发生背景

执行WGCNA代码时报错，本段代码拟计算Beta值

问题相关代码，请勿粘贴截图

powers = c(c(1:10), seq(from = 12, to = 20, by = 2))
sft = pickSoftThreshold(datExpr, powerVector = powers, verbose = 5)

pdf("./WGCNA/1_beta.pdf",width = 8,height = 5)
par(mfrow = c(1, 2))
cex1 = 0.9
plot(
    sft$fitIndices[, 1],
    -sign(sft$fitIndices[, 3]) * sft$fitIndices[, 2],
    xlab = "Soft Threshold (power)",
    ylab = "Scale Free Topology Model Fit,signed R^2",
    type = "n",
    main = paste("Scale independence")
)

text(
    sft$fitIndices[, 1],
    -sign(sft$fitIndices[, 3]) * sft$fitIndices[, 2],
    labels = powers,
    cex = cex1,
    col = "red"
)

运行结果及报错内容

text(

• sft$fitIndices[, 1],
  

• -sign(sft$fitIndices[, 3]) * sft$fitIndices[, 2],
  

• labels = powers,
  

• cex = cex1,
  

• col = "red"
  

• )
  Error in text.default(sft$fitIndices[, 1], -sign(sft$fitIndices[, 3]) * :
  plot.new has not been called yet

我的解答思路和尝试过的方法

在网上检索，包括本站内，多数提示plot.new has not been called yet会在旧文件被删除的情况下出现，导致新文件无法建立。看到好几例问答都关于plot绘图，或者abline函数。目前报错的位置，下一行才是abline函数

我想要达到的结果

学习报错出现的原因，以及解决方法。
谢谢

你好，我是有问必答小助手，非常抱歉，本次您提出的有问必答问题，技术专家团超时未为您做出解答


本次提问扣除的有问必答次数，将会以问答VIP体验卡（1次有问必答机会、商城购买实体图书享受95折优惠）的形式为您补发到账户。


因为有问必答VIP体验卡有效期仅有1天，您在需要使用的时候【私信】联系我，我会为您补发。