原文链接：http://tecdat.cn/?p=9227

原文出处：拓端数据部落公众号

数据集：行为危险因素监视系统数据

摘要：该数据集是来自全美约40万份与衰弱相干主题的问卷调查。BRFSS始于1980年代，并已通过问卷调查在美国用于监测广泛的疾病。该钻研是追溯性的，而不是设计性的试验，因而只管能够推断出相关性，但不能因果关系。

数据集中的特色既是间断的又是分类的。

指标：摸索性别，体重和年龄之间的相关性
 

第0局部：设置

library(ggplot2)library(dplyr)library(Rgraphviz)library(knitr)library(grid)library(gridExtra)
load("brfss2013.RData")
# group and count a feature with discrete valuesfeature_vcounts <- function(df, f) {  df %>%    group\_by\_at(f) %>%      count()}# method for binning valuesbin\_min\_sample <- function(p) {  n = 10  a = 10/p  b = 10/(1-p)  max(c(a,b))}# create a new df for simulating binom probability distributionbinom\_prob\_df <- function(df, f, target) {  new\_df <- feature\_vcounts(df,f)  new\_df$n\[new\_df\[f\] == target\]/sum(new_df$n)}# filtering df with subgroup valuesubgroup_df <- function(df,f, group) {  filter(df,df\[f\]==group)}# calc the vector probabilitybinom\_prob\_vec <- function(v, target) {  sum(v == target)/length(v)}# sample from dfbinom_sample <- function(s,v)  sample(v, size=s, replace=TRUE)# create the binomial sample distributionbinom\_sample\_dist <- function(df,f,target) {  sample_dist <- c()  for (i in 1:10001) {    prob <- binom\_prob\_vec(binom_sample(100,df\[,f\]),target)    sample\_dist <- append(sample\_dist,prob)}  return(sample_dist)}# convert decimal to percentto_percent <- function(pvalue) {  paste(round(pvalue*100,digits= 2),"%",sep="")}

第1局部：数据

导入和过滤数据以仅包含与糖尿病，性别，体重和年龄无关的重要特色。

# Import original file:orig_dim <- dim(brfss2013)# Select only relevant features:weight_diabetes <- brfss2013 %>%   select(sex, X_ageg5yr, weight2,diabete3)# ------------------Cleaning data------------------# 1.Weight strings -> numericweight\_diabetes$weight2 <- as.numeric(as.character(weight\_diabetes$weight2))new\_dim <- dim(weight\_diabetes)# 2. Remove Null Weights and Weights over 400weight\_diabetes <- na.omit(weight\_diabetes)weight\_diabetes <- filter(weight\_diabetes, weight2 <= 400)# 3. Remove Diabetes Responsestarget <- c("Yes", "No")weight\_diabetes <- filter(weight\_diabetes, diabete3 %in% target)# 4. Add index and reorderweight\_diabetes$index <- seq.int(nrow(weight\_diabetes))weight\_diabetes <- weight\_diabetes\[c(4,3,1,2)\]clean\_dim <- dim(weight\_diabetes)# Show data:kable(head(weight_diabetes,n=5), caption="Diabetes Data Set",padding=0, format = "markdown",align="l")

糖尿病 

分量 

性别

X_ageg5yr

没有

250

女

60至64岁

没有

127

女

50至54岁

没有

160

女

55至59岁

没有

128

女

60至64岁

没有

265

男

65至69岁

数据看起来很简略，仅蕴含该我的项目所需的性能。因为数据须要匿名，所以年龄范畴是特定年龄的平安代替计划。年龄范畴将用作此数据集的分类信息。

第2局部：钻研问题

钻研问题1：

性别，体重和年龄之间有相关性吗？（变量：性别，weight2，X_ageg5yr）

因为性别是生物辨认技术中的要害变量，因而探讨性别是否可能与其余变量相干很重要。在这种状况下，咱们正在钻研性别是否与体重相干。

钻研问题2：

性别或年龄与糖尿病相干吗？怎么样？（变量：性别，X_ageg5yr，weight2，diabete3）

该探索性我的项目的指标是查看体重/性别/年龄是否与糖尿病相干。理解任何相关性可能有助于依据患者的性别和体重告知患者患糖尿病的可能性。

钻研问题3：

年龄，体重和糖尿病之间有关系吗？（变量：性别，X_ageg5yr，weight2，diabete3）

为了进一步探讨与糖尿病的可能相关性，咱们还将钻研四个变量之间的关系。

第3局部：探索性数据分析

钻研问题1：

性别，体重和年龄之间有相关性吗？（变量：性别，weight2，X_ageg5yr）

首先检查数据的散布很重要。_性别_是二元分类的，因而咱们将用条形图形象化它的散布。

centered <- theme(plot.title = element_text(hjust = 0.5))hist\_weight <- ggplot(data=weight\_diabetes,aes(weight2,  fill=weight2))+  geom_histogram(fill='salmon',color='white') + ggtitle("Histogram \[Weight\]") + centeredweight\_diabetes$log\_weight <- log(weight_diabetes$weight2)hist\_log\_weight <- ggplot(data=weight\_diabetes,aes(log\_weight, fill=log_weight))+  geom\_histogram(fill='mediumturquoise',color='white') + ggtitle("Histogram \[Log\_Weight\]") + centeredgrid.arrange(hist\_weight, hist\_log_weight, ncol = 2)


对于_Weight2_，散布右偏，而_Weight2的_对数靠近正态。因为数据的对数标准版本简直是失常的单峰数据，因而能够将权重用于推断统计中的后续剖析。

 
女性参加者比男性参加者更多，其幅度大大超过美国的总人口。这可能表明抽样办法在性别抽样方面并非齐全随机。然而，数据样本足够大，能够持续评估衰弱危险因素。

年龄范畴仿佛在两端都偏差极其。

在比拟年龄和体重时，性别的体重散布仿佛的确存在显著差别。男性仿佛比女性重。

还应留神，仿佛还有一些成年患者体重有余50磅，许多患者体重约400磅。数据收集过程的将来剖析应摸索这些高下异样值是否是谬误，或者它们是否反映出患有重大衰弱问题的患者。

钻研问题2：

体重，年龄和/或性别与糖尿病相干吗？怎么样？（变量：性别，X_ageg5yr，weight2，diabete3）

 
当察看样本中的女性和男性参与者时，报告的糖尿病比率十分类似。

随着年龄的增长，直到80岁及更高年龄，糖尿病的报告比例也仿佛会减少。年龄可能与糖尿病有某种程度的相关性。

随着体重增加，所报告的糖尿病比例也减少。体重仿佛与报告的糖尿病有很强的相关性，应进一步探讨。

钻研问题3：

年龄，体重和糖尿病之间有关系吗？（变量：性别，X_ageg5yr，weight2，diabete3）

 
报告患有糖尿病的男性和女性的体重散布均较高。男性的体重散布比女性重。

这些箱线图描述了一个更清晰的故事。报告患有糖尿病的患者仿佛在每个年龄段都较重。报告患有糖尿病的年老患者仿佛比老年患者具备更大的体重范畴。尽管尚不分明年龄与糖尿病和体重之间的关系，但应进一步探讨这种关系。

第4局部：论断

从数据的初步摸索中能够显著看出，某些性能具备比其余性能更强的相关性。体重与性别无关。性别仿佛与体重无关。然而，糖尿病仿佛与年龄无关，而与体重密切相关。