在 R corr.test 中找不到对象'sef'

Question

我试图以我的教授在他的系统上提交和测试的代码来运行cor.test方程。不幸的是，当我运行它时，我会遇到"未找到对象sef"的错误。

这使我和我的教授都感到困惑，并且在进行了彻底的搜索之后，我们不确定如何解决这个问题。

我非常感谢您提供的任何帮助。

编辑：这是我正在使用的代码：

trendan1 <- read.table("trendan1.for.R.dat", header=TRUE, na.strings=".")
head(trendan1)
tail(trendan1)
attributes(trendan1)
is.matrix(trendan1)
id <- trendan1$id
famenv1 <- trendan1$famenv1
famenv2 <- trendan1$famenv2
famenv3 <- trendan1$famenv3
conf1 <- trendan1$conf1
conf2 <- trendan1$conf2
conf3 <- trendan1$conf3
trendan1dataset1 <- cbind(id,famenv1,famenv2,famenv3,conf1,conf2,conf3)
attributes(trendan1dataset1)
is.matrix(trendan1dataset1)
is.data.frame(trendan1dataset1)
require("psych")
describe(trendan1dataset1[,2:7])
print(describe(trendan1dataset1[,2:7]), digits=6)  

famave <-  (1*famenv1 + 1*famenv2 + 1*famenv3)/3
famlin <-  -1*famenv1 + 0*famenv2 + 1*famenv3
famquad <-  1*famenv1 - 2*famenv2 + 1*famenv3;
trendandataset2 <- cbind(famenv1,famenv2,famenv3,famave,famlin,famquad)
print(describe(trendandataset2), digits=6)
hist(famenv1)
boxplot(famenv1)
abline(h=mean(famenv1)) 
qqnorm(famenv1,ylab="famenv1")  
qqline(famenv1)                
shapiro.test(famenv1)
hist(famenv2)
boxplot(famenv2)
abline(h=mean(famenv2)) # add mean to the boxplot
qqnorm(famenv1,ylab="famenv2")  
qqline(famenv2)                
shapiro.test(famenv2)
corvars1 <- cbind(famenv1,famenv2,famenv3)

cor(corvars1,use = "everything", method = "pearson")
cov(corvars1,use = "everything")
sscp1 <- t(corvars1)%*%(corvars1) #Matrix multiplcation
sscp1
rc1 <- corr.test(corvars1, 
use="pairwise",method="pearson",adjust="holm",alpha=.05, ci=FALSE) 
attributes(rc1)
print(rc1$p, digits=6)

Answer 1

这是一个错误，当您不评估置信区间时，有时会发生。如果您将选项更改为ci = true ，则应修复它，或者简单地删除此选项，因为默认值为ci = true。