Question

我试图在R中运行corr.test等式，我的教授在他的系统上提交并测试了代码。不幸的是，当我运行它时，我收到一个“找不到对象的错误”。

这让我和我的教授感到困惑，并且做了彻底的搜索，我们不确定如何解决这个问题。

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

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

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。

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

1 个答案: