我不能按性别绘制我的qq图,但可以将它们组合起来。它一直在说:“ y为空或仅具有NA”,而且我似乎无法弄清楚应该如何更改代码。您可以在下面看到我的一些代码:
## 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不会让我插入我的整个代码。谢谢