我无法为单个性别绘制我的 qq 图,但我可以将它们组合起来。它一直在说:“y 是空的或只有 NAs”,我似乎无法弄清楚我应该如何更改我的代码。您可以在下面看到我的一些代码:
## Read data
D<-read.table("bmi1_data.csv",header=TRUE,sep=";",as.is=TRUE)
## Divide data into two subsets according to gender
Dfemale<-subset(D,gender == 0)
Dmale<-subset(D,gender==1)
我在代码的前面已经计算了 BMI。
## New variable 'logbmi' with log-transformed BMI
D$logbmi <- log(D$bmi)
它可以绘制这个:
## qq-plot of log-transformed BMI
qqnorm(D$logbmi)
qqline(D$logbmi)
但它不能绘制这个:
## qq-plot of log-transformed BMI for females
qqnorm(Dfemale$logbmi)
qqline(Dfemale$logbmi)
这也不起作用:
## qq-plot of log-transformed BMI for males
qqnorm(Dmale$logbmi)
qqline(Dmale$logbmi)
我希望这是有道理的,因为 stackoverflow 不允许我插入整个代码。谢谢