简单回归但出现错误

问题描述

我需要做一个简单的二元回归。

我的脚本是:

staff_data<-read.table("staff_data (2).txt",header = T,sep = "\t",dec=".")

install.packages("Hmisc")
library(Hmisc)
library(psych)

developing<-subset(staff_data,(developing==1))
income<-subset(staff_data,(incomecategory==2))

#subsets,mean,median,and standard deviation - Part 1
Q1<-subset(staff_data,(mobile>=21.2999913817616 & mobile<=85.6747730239075))


mean(Q1$measlessimunization)
sd(Q1$measlessimunization)
median(Q1$measlessimunization)


integrate.xy(r1density$x,r1density$y,60,110)

#t-test - Part 2
Q2<-subset(staff_data,(mobile>=17.6560769700962 & mobile<=92.6889722787938))
developed<-subset(Q2,developing==0)
mean(developed$measlessimunization)
developing<-subset(Q2,developing==1)
mean(developing$measlessimunization)
x=developed$measlessimunization
y=developing$measlessimunization


t.test(x,y)

#create variable and regression - Part 3

staff_data$logstaffwomen<-log10(staff_data$staff_women_share)
staff_data$loggnipc<-log10(staff_data$gnipc)

Q3<-subset(staff_data,(mobile>=33.8045645560072 & mobile<=100.091284733476))

hist(staff_data$logstaffwomen)
hist(staff_data$loggnipc)

plot(staff_data$loggnipc,staff_data$logstaffwomen,main="Staff,given by population (all logs)",ylab="Logarithm of professional FEMALE staff",xlab="Logarithm of GNIPC")
text(staff_data$loggnipc,staff_data$countrycode,cex=0.6,pos=4,col="red") 

abline(lm(logstaffwomen~loggnipc,data=staff_data),col="darkblue") # regression line (y~x)
abline(lm(staff_data$logstaffwomen~staff_data$loggnipc),col="darkblue") # regression line (y~x)

mybeautifulmodel<-lm(loggnipc~staff_women_share,data=staff_data)
summary(mybeautifulmodel)

错误

abline(lm(logstaffwomen~loggnipc,col="darkblue") # 回归线 (y~x) lm.fit(x,y,offset = offset,single.ok = single.ok,...) 中的错误: 'y' 中的 NA/NaN/Inf abline(lm(staff_data$logstaffwomen~staff_data$loggnipc),...) 中的错误: 'y' 中的 NA/NaN/Inf 总结(我的美丽模型) 总结错误(mybeautifulmodel):找不到对象“mybeautifulmodel”

请问,有人能解释一下我做错了什么吗?数据集错了吗?

解决方法

暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!

如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@)