简单回归但出现错误

如何解决简单回归但出现错误

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

我的脚本是：

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”

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