排序中的R

Question

我建立逻辑回归模型在x轴（称为“mylogit”）使用GLM功能中的R，如下所示：

mylogit <- glm(answer ~ as.factor(gender) + age, data = mydata, family = "binomial") 

其中年龄是数字和性别分类（男性和女性） 。

然后我开始用建好的模型做出预测。

pred <- predict(mylogit, type = "response") 

我可以很容易地通过做使预测的时间序列图：

plot.ts(ts(pred)) 

给一个情节，看起来像这样：

Plot of Time against Predictions

这给出了一个阴谋的预测。

我的问题是这样的： 根据glm中指定的性别（男性或女性）是否可以将x轴放入段中？换句话说，我可以预测y轴，并且在x轴上有性别（分为男性和女性）吗？

我想从绘制数据的样本与此类似：

我所做的：

bind = cbind(mydata, pred) 

'捆绑' 看起来是这样的：

pred   age  gender 
0.9461198  32   male 
0.9463577  45   female 
0.9461198  45   female 
0.9461198  37   female 
0.9477645  40   male 
0.8304513  32   female 

Answer 1

退房＃ 4，this blog post，“4.如何创建两个不同的X轴或Y轴”。

我给你的建议是，你看看一些专用的R绘图工具，如ggplot2。

Answer 2

我不认为你需要使用ts和plot.ts，因为你拥有的数据不是时间序列，对不对？在绘图之前对pred进行排序。

# Get data 
str <- "pred,age,gender 
0.9461198,32,male 
0.9463577,45,female 
0.9461198,45,female 
0.9461198,37,female 
0.9477645,40,male 
0.8304513,32,female" 
bind <- read.csv(textConnection(str)) 

# Plot 
bind <- bind[order(bind$gender),] 
plot(bind$pred, col = bind$gender) 

library(ggplot2) 
ggplot(bind, aes(x = gender, y = pred)) + 
    geom_point(position = position_jitter(width = .3)) 

或者不创建bind你可以做plot(pred[order(mydata$gender)])。