Body Shape Data

require(devtools)

## Loading required package: devtools

install_url("http://mosaic-web.org/go/Repository/DannyKaplan/DCF_0.11.tar.gz")

## Installing DCF_0.11.tar.gz from
## http://mosaic-web.org/go/Repository/DannyKaplan/DCF_0.11.tar.gz

## Installing DCF

## pdflatex not found. Not building PDF vignettes.

## '/usr/lib/R/bin/R' --vanilla CMD build '/tmp/RtmpWTAA8N/DCF' --no-manual \
## --no-resave-data --no-vignettes

## 

## '/usr/lib/R/bin/R' --vanilla CMD INSTALL \
## '/tmp/RtmpWTAA8N/DCF_0.11.tar.gz' --library='/home/azaidi/R/library' \
## --with-keep.source

## 

require(DCF)

## Loading required package: DCF

## Loading required package: mosaic

## Loading required package: lattice

## Loading required package: grid

## Loading required package: survival

## Loading required package: splines

## Loading required package: Hmisc

## Hmisc library by Frank E Harrell Jr
## 
## Type library(help='Hmisc'), ?Overview, or ?Hmisc.Overview') to see overall
## documentation.
## 
## NOTE:Hmisc no longer redefines [.factor to drop unused levels when
## subsetting.  To get the old behavior of Hmisc type dropUnusedLevels().

## Attaching package: 'Hmisc'

## The following object(s) are masked from 'package:survival':
## 
## untangle.specials

## The following object(s) are masked from 'package:base':
## 
## format.pval, round.POSIXt, trunc.POSIXt, units

## Attaching package: 'mosaic'

## The following object(s) are masked from 'package:Hmisc':
## 
## do

## The following object(s) are masked from 'package:stats':
## 
## D, binom.test, median, prop.test, sd, var

## The following object(s) are masked from 'package:base':
## 
## max, mean, min, print, sample

## Loading required package: ggplot2

## Loading required package: reshape2

## Loading required package: plyr

## Attaching package: 'plyr'

## The following object(s) are masked from 'package:mosaic':
## 
## count

## The following object(s) are masked from 'package:Hmisc':
## 
## is.discrete, summarize

data(nhanes)
small = sample(nhanes, 2000)

Describe the relationship between height and weight. Is there reason to think that it's different for the two sexes?

ggplot(data = nhanes) + geom_point(aes(x = hgt, y = wgt)) + aes(colour = sex)

## Warning: Removed 4967 rows containing missing values (geom_point).

Describe the relationship between weight and BMI. Is it different from the two sexes? Where do the people with diabetes show up?

ggplot(data = nhanes) + geom_point(aes(x = wgt, y = bmi)) + aes(colour = sex) + 
    aes(size = diab)

## Warning: Removed 5176 rows containing missing values (geom_point).

Is there a relationship between BMI and age? Where do the people with diabetes show up?

ggplot(data = nhanes) + geom_point(aes(x = age, y = bmi)) + aes(colour = diab)

## Warning: Removed 4967 rows containing missing values (geom_point).

Is cholesterol level a good predictor of the development of diabetes?

ggplot(data = nhanes) + geom_point(aes(x = chol, y = diab))

## Warning: Removed 7839 rows containing missing values (geom_point).