如何从R子集中删除异常值？

我对R非常陌生，正试图从一个子集中移除异常值以改善GLM。我使用的代码是：如何从R子集中删除异常值？

data$sel <- ifelse(data$chol==8.3 & data$whr==1.14 ,(0), (1)) 
data 
dim(data) 
data2 <- subset(data, !(chol==8.3 & whr==1.14)) 
dim(data2)

我这样做，但是，当我试图绘制新的数据图表，问题依然出现。我用来绘制新图的代码是：

dataF2 <- subset(data2, sex=="F") 
dataF2 
dataM2 <- subset(data2, sex=="M") 
dataM2 

plot(chol ~ whr, data=data2, type="n", ylab="Cholesterol (mM/L)", xlab="Waist-hip Ratio") 

#### add each group of points: 
points(dataF2$chol ~ dataF2$whr, pch=1, col="red") 
points(dataM2$chol ~ dataM2$whr, pch=2, col="blue") 

Females2 <- lm(chol ~ whr, data=dataF2) 
summary(Females) 
abline(-8.053, 14.801, col="red") 
Males2 <- lm(chol ~ whr, data=dataM2) 
summary(Males2) 
abline(-3.5896, 8.5617, col="blue")

我不知道我在做什么错，或者我错过了一些重要的步骤。我所要做的只是删除一个或多个异常值，以便改进GLM，然后生成没有异常值的图。

数据

'data.frame': 100 obs. of 7 variables: 
$ age : int 42 41 40 43 61 44 39 70 47 55 ... 
$ sex : Factor w/ 2 levels "F","M": 1 1 1 1 1 1 2 1 1 1 ... 
$ bmi : num 28.5 26.2 25.2 29.3 34.6 26.7 25.5 33.2 28.8 31.4 ... 
$ whr : num 0.83 0.82 0.84 0.82 0.89 0.81 0.93 0.97 0.86 0.88 ... 
$ sysbp: int 132 135 141 142 167 145 137 182 150 160 ... 
$ chol : num 4.2 3.3 4.3 4.1 5.5 4.4 3.9 6.7 4.7 5.2 ... 
$ smoke: Factor w/ 2 levels "N","Y": 1 2 2 2 1 2 1 2 2 2 ...

数据2：

'data.frame': 99 obs. of 8 variables: 
$ age : int 42 41 40 43 61 44 39 70 47 55 ... 
$ sex : Factor w/ 2 levels "F","M": 1 1 1 1 1 1 2 1 1 1 ... 
$ bmi : num 28.5 26.2 25.2 29.3 34.6 26.7 25.5 33.2 28.8 31.4 ... 
$ whr : num 0.83 0.82 0.84 0.82 0.89 0.81 0.93 0.97 0.86 0.88 ... 
$ sysbp: int 132 135 141 142 167 145 137 182 150 160 ... 
$ chol : num 4.2 3.3 4.3 4.1 5.5 4.4 3.9 6.7 4.7 5.2 ... 
$ smoke: Factor w/ 2 levels "N","Y": 1 2 2 2 1 2 1 2 2 2 ... 
$ sel : num 1 1 1 1 1 1 1 1 1 1 ...

来源

2017-11-25 Anna

您的代码不是[重现]（https://stackoverflow.com/help/mcve）。请提供您的数据或其中的一些最小形式。 – SirSaleh

您的代码仍然不可重现。塞德里克的回答有帮助吗？ – AkselA

到子集数据，您可以通过TRUE FALSE的矢量线条。

你的数据帧是数据2 [行，列]

data2 <- data2[data2$chol!=8.3 & data2$whr!=1.14,]

读this r blogger example更好地了解如何从您的数据帧删除一些行。

你也刚使用

data2 <- subset(data, (chol!=8.3 & whr!=1.14))

或

data2 <- subset(data, !(chol==8.3 | whr==1.14))

我的猜测是，你有没有线，其中既哲和腰臀比有那些价值观，要删除两个不同行。

最后建议，请您传递的东西，利用第一配方可以让你检查的代码位，

data2$chol!=8.3 & data2$whr!=1.14

，并检查是否返回TRUE FALSE的矢量你希望

来源

2017-11-25 20:58:09 Cedric

如何从R子集中删除异常值？

回答

相关问题