2013-04-20 145 views
0

我想绘制y轴上的连续BMI和x轴上家庭收入的分类变量,我希望图形绘制每个类别的平均BMI。但是,我不确定如何找到每个家庭收入因素的平均BMI。如何找到每个分类变量的连续变量的均值

Dataset nh (5994 total IDs with Observations) (Parts of the 2009-2010 NHANES Dataset) 
> dput(head(nh)) 
structure(list(SeqN = c(51624L, 51628L, 51629L, 51630L, 51633L, 51635L), 
Gender = c(1L, 2L, 1L, 2L, 1L, 1L), Age = c(34L, 60L, 26L, 49L, 80L, 80L), 
Ethnicity = c(3L, 4L, 1L, 3L, 3L, 3L), FamSize = c(4L, 2L, 5L, 3L, 2L, 1L), 
RatioIncomePoverty = c(1.36, 0.69, 1.01, 1.91, 1.27, 1.69), 
MECWgt2 = c(81528.77201, 21000.33872, 22633.58187, 74112.48684, 12381.11532, 22502.50666), 
BMI = c(32.22, 42.39, 32.61, 30.57, 26.04, 27.62), 
LengthUS = c(NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_,NA_integer_), 
Education = c(3L, 3L, 2L, 4L, 4L, 2L), LocationBorn = c(1L, 1L, 1L, 1L, 1L, 1L), 
FamIncome = c(6L, 3L, 6L, 7L, 4L, 4L)), .Names = c("SeqN", 
"Gender", "Age", "Ethnicity", "FamSize", "RatioIncomePoverty", 
"MECWgt2", "BMI", "LengthUS", "Education", "LocationBorn", "FamIncome"), 
row.names = c(NA, 6L), class = "data.frame") 

faminc <- as.character(nhanes$FamIncome) 
faminc 

有关如何建模数据以实现此目标的任何建议将不胜感激。

+2

搞什么名堂你试过这么远吗?你可以发布你已经尝试过的东西吗?你有什么困难的时候找出答案?如果不是的话,我只能假设这是一个为我而设的问题。 – Arun 2013-04-20 16:14:02

+2

另外,通过粘贴'dput(head(nh))'的输出来发布你的数据,这样我们可以将它复制/粘贴到我们的R会话中,并且它将具有相同的结构(并且我们可以告诉你是否有'FamInc'作为一个因素或数字)。 – Gregor 2013-04-20 16:20:16

+0

@Arun我不是在寻找让我成为整个图形或“为我而做”的人,而是一些方向,因为我不确定我能够适当地说出我想要做什么基本的谷歌搜索。我想我最大的问题是如何建模我的数据或者什么功能可以为每类家庭收入创建一个平均BMI的数字变量。我编辑了原稿以包含更多信息。 – MryElln 2013-04-20 18:31:49

回答

2

下面是使用aggregate底液:

a <- aggregate(BMI ~ FamIncome, data=nh, FUN=mean) 
barplot(a$BMI, names.arg=a$FamIncome) 

enter image description here

2

这可能工作:

library(plyr) 
    nhh<-ddply(nh,.(famIncome), summarise, mean.bmi=mean(bmi)) # find mean bmi 
    with(nhh, plot(famIncome,mean.bmi)) # simple plot 
+0

我正在使用ddply,但我想方式混淆它,我想。这似乎已经完成了诀窍,给出了每个家庭收入因素的平均体重指数。希望这将是我需要制作我想要的情节。谢谢。 – MryElln 2013-04-20 19:23:22

+0

我已经更新了情节答案(这只适用于简单情节)。如果您确信,请检查接受的答案。 – Metrics 2013-04-20 19:25:12

+0

这有效。感谢附加的情节,它的工作原理虽然我不需要它!很高兴看到你是如何做到的,所以我可以将它实现到我的ggplot2图中。我会将其标记为已回答。 – MryElln 2013-04-20 19:30:57