Health Data Analysis in R

Shakeel Ahmed

2024-12-26

R Markdown

-R Markdown see http://rmarkdown.rstudio.com.

-When you click the Knit button a document will be generated that includes both content as well as the output of any embedded R code chunks within the document.

In this Lecture we see how to analyze data in R (NHANES dataset)

How to import data from stata

library(haven)
 NHANES_2005_18_with_mortality <- read_dta("C:/Users/Admin/Downloads/Dr Shakeel/NHANES_2005-18_with_mortality.dta")

#View Data

## tibble [70,190 × 958] (S3: tbl_df/tbl/data.frame)
##  $ seqn        : num [1:70190] 31127 31128 31129 31130 31131 ...
##   ..- attr(*, "label")= chr "Respondent sequence number"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ sddsrvyr    : num [1:70190] 4 4 4 4 4 4 4 4 4 4 ...
##   ..- attr(*, "label")= chr "Data Release Number"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ ridstatr    : num [1:70190] 2 2 2 2 2 2 2 2 2 1 ...
##   ..- attr(*, "label")= chr "Interview/Examination Status"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ ridexmon    : num [1:70190] 2 1 2 2 2 2 2 2 2 NA ...
##   ..- attr(*, "label")= chr "Six month time period"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ riagendr    : num [1:70190] 1 2 1 2 2 1 2 1 1 2 ...
##   ..- attr(*, "label")= chr "Gender"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ ridageyr    : num [1:70190] 0 11 15 85 44 70 16 73 0 41 ...
##   ..- attr(*, "label")= chr "Age at Screening Adjudicated - Recode"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ ridagemn    : num [1:70190] 11 132 189 NA 535 842 193 882 10 493 ...
##   ..- attr(*, "label")= chr "Age in Months - Recode"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ ridageex    : num [1:70190] 12 132 190 NA 536 843 194 883 11 NA ...
##   ..- attr(*, "label")= chr "Exam Age in Months - Recode"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ ridreth1    : num [1:70190] 3 4 4 3 4 3 4 3 5 4 ...
##   ..- attr(*, "label")= chr "Race/Ethnicity - Recode"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ dmqmilit    : num [1:70190] NA NA NA 2 2 1 NA 1 NA 2 ...
##   ..- attr(*, "label")= chr "Veteran/Military Status"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ dmdborn     : num [1:70190] 1 1 1 1 1 1 1 1 1 1 ...
##   ..- attr(*, "label")= chr "Country of Birth - Recode"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ dmdcitzn    : num [1:70190] 1 1 1 1 1 1 1 1 1 1 ...
##   ..- attr(*, "label")= chr "Citizenship Status"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ dmdyrsus    : num [1:70190] NA NA NA NA NA NA NA NA NA NA ...
##   ..- attr(*, "label")= chr "Length of time in US"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ dmdeduc3    : num [1:70190] NA 4 10 NA NA NA 9 NA NA NA ...
##   ..- attr(*, "label")= chr "Education Level - Children/Youth 6-19"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ dmdeduc2    : num [1:70190] NA NA NA 4 4 5 NA 3 NA 4 ...
##   ..- attr(*, "label")= chr "Education Level - Adults 20+"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ dmdschol    : num [1:70190] NA 1 1 NA NA NA 1 NA NA NA ...
##   ..- attr(*, "label")= chr "Now attending school?"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ dmdmartl    : num [1:70190] NA NA 5 2 1 1 5 1 NA 5 ...
##   ..- attr(*, "label")= chr "Marital Status"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ dmdhhsiz    : num [1:70190] 4 7 6 1 4 2 3 2 3 1 ...
##   ..- attr(*, "label")= chr "Total number of people in the Household"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ dmdfmsiz    : num [1:70190] 4 6 6 1 4 2 3 2 3 1 ...
##   ..- attr(*, "label")= chr "Total number of people in the Family"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ indhhinc    : num [1:70190] 4 8 10 4 11 11 11 12 11 7 ...
##   ..- attr(*, "label")= chr "Annual Household Income"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ indfminc    : num [1:70190] 4 5 10 4 11 11 11 12 11 7 ...
##   ..- attr(*, "label")= chr "Annual Family Income"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ indfmpir    : num [1:70190] 0.75 0.77 2.71 1.99 4.65 5 5 NA 5 3.57 ...
##   ..- attr(*, "label")= chr "Family PIR"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ ridexprg    : num [1:70190] NA 2 NA NA 2 NA 2 NA NA 3 ...
##   ..- attr(*, "label")= chr "Pregnancy Status at Exam - Recode"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ dmdhrgnd    : num [1:70190] 2 1 1 2 1 1 1 1 1 2 ...
##   ..- attr(*, "label")= chr "HH Ref Person Gender"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ dmdhrage    : num [1:70190] 21 47 41 85 36 70 47 73 31 41 ...
##   ..- attr(*, "label")= chr "HH Ref Person Age"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ dmdhrbrn    : num [1:70190] 1 1 1 1 1 1 1 1 1 1 ...
##   ..- attr(*, "label")= chr "HH Ref Person Country of Birth"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ dmdhredu    : num [1:70190] 3 2 4 4 5 5 2 3 5 4 ...
##   ..- attr(*, "label")= chr "HH Ref Person Education Level"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ dmdhrmar    : num [1:70190] 1 NA 1 2 1 1 1 1 1 5 ...
##   ..- attr(*, "label")= chr "HH Ref Person Marital Status"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ dmdhsedu    : num [1:70190] 2 NA 4 NA 4 5 4 3 5 NA ...
##   ..- attr(*, "label")= chr "HH Ref Person's Spouse Education Level"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ sialang     : num [1:70190] 1 1 1 1 1 1 1 1 1 1 ...
##   ..- attr(*, "label")= chr "Language of SP Interview"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ siaproxy    : num [1:70190] 1 1 1 2 2 2 2 2 1 2 ...
##   ..- attr(*, "label")= chr "Proxy used in SP Interview?"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ siaintrp    : num [1:70190] 2 2 2 2 2 2 2 2 2 2 ...
##   ..- attr(*, "label")= chr "Interpreter used in SP Interview?"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ fialang     : num [1:70190] 1 1 1 1 1 1 1 1 1 1 ...
##   ..- attr(*, "label")= chr "Language of Family Interview"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ fiaproxy    : num [1:70190] 2 2 2 2 2 2 2 2 2 2 ...
##   ..- attr(*, "label")= chr "Proxy used in Family Interview?"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ fiaintrp    : num [1:70190] 2 2 2 2 2 2 2 2 2 2 ...
##   ..- attr(*, "label")= chr "Interpreter used in Family Interview?"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ mialang     : num [1:70190] NA 1 1 NA 1 1 1 1 NA NA ...
##   ..- attr(*, "label")= chr "Language of MEC Interview"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ miaproxy    : num [1:70190] NA 2 2 NA 2 2 2 2 NA NA ...
##   ..- attr(*, "label")= chr "Proxy used in MEC Interview?"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ miaintrp    : num [1:70190] NA 2 2 NA 2 2 2 2 NA NA ...
##   ..- attr(*, "label")= chr "Interpreter used in MEC Interview?"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ aialang     : num [1:70190] NA 1 1 NA 1 1 1 1 NA NA ...
##   ..- attr(*, "label")= chr "Language of ACASI Interview"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ wtint2yr    : num [1:70190] 6435 9082 5317 29961 26458 ...
##   ..- attr(*, "label")= chr "Full Sample 2 Year Interview Weight"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ wtmec2yr    : num [1:70190] 6571 8987 5587 34031 26771 ...
##   ..- attr(*, "label")= chr "Full Sample 2 Year MEC Exam Weight"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ sdmvpsu     : num [1:70190] 2 1 1 2 1 2 1 2 2 2 ...
##   ..- attr(*, "label")= chr "Masked Variance Pseudo-PSU"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ sdmvstra    : num [1:70190] 44 52 51 46 48 52 51 48 52 51 ...
##   ..- attr(*, "label")= chr "Masked Variance Pseudo-Stratum"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ urxuma      : num [1:70190] NA 82.1 17.8 NA 18 6.5 25.2 10.1 NA NA ...
##   ..- attr(*, "label")= chr "Albumin, urine (ug/mL)"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ urxums      : num [1:70190] NA 82.1 17.8 NA 18 6.5 25.2 10.1 NA NA ...
##   ..- attr(*, "label")= chr "Albumin, urine (mg/L)"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ urxucr      : num [1:70190] NA 291 288 NA 202 162 275 140 NA NA ...
##   ..- attr(*, "label")= chr "Creatinine, urine (mg/dL)"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ urxcrs      : num [1:70190] NA 25724 25459 NA 17857 ...
##   ..- attr(*, "label")= chr "Creatinine, urine (umol/L)"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ alq101      : num [1:70190] NA NA NA NA 2 1 NA 1 NA NA ...
##   ..- attr(*, "label")= chr "Had at least 12 alcohol drinks/1 yr?"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ alq110      : num [1:70190] NA NA NA NA 1 NA NA NA NA NA ...
##   ..- attr(*, "label")= chr "Had at least 12 alcohol drinks/lifetime?"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ alq120q     : num [1:70190] NA NA NA NA 0 4 NA 2 NA NA ...
##   ..- attr(*, "label")= chr "How often drink alcohol over past 12 mos"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ alq120u     : num [1:70190] NA NA NA NA NA 1 NA 1 NA NA ...
##   ..- attr(*, "label")= chr "# days drink alcohol per wk, mo, yr"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ alq130      : num [1:70190] NA NA NA NA NA 1 NA 2 NA NA ...
##   ..- attr(*, "label")= chr "Avg # alcoholic drinks/day -past 12 mos"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ alq140q     : num [1:70190] NA NA NA NA NA 0 NA 0 NA NA ...
##   ..- attr(*, "label")= chr "#days have 5 or more drinks/past 12 mos"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ alq140u     : num [1:70190] NA NA NA NA NA NA NA NA NA NA ...
##   ..- attr(*, "label")= chr "# days per week, month, year?"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ alq150      : num [1:70190] NA NA NA NA 2 2 NA 2 NA NA ...
##   ..- attr(*, "label")= chr "Ever have 5 or more drinks every day?"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmdstats    : num [1:70190] 1 1 1 4 3 1 1 1 1 NA ...
##   ..- attr(*, "label")= chr "Body Measures Component Status Code"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmxwt       : num [1:70190] 10.2 40.1 74.6 NA 75.2 ...
##   ..- attr(*, "label")= chr "Weight (kg)"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmiwt       : num [1:70190] NA NA NA NA NA NA NA NA NA NA ...
##   ..- attr(*, "label")= chr "Weight Comment"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmxrecum    : num [1:70190] 73.6 NA NA NA NA NA NA NA 74.8 NA ...
##   ..- attr(*, "label")= chr "Recumbent Length (cm)"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmirecum    : num [1:70190] NA NA NA NA NA NA NA NA NA NA ...
##   ..- attr(*, "label")= chr "Recumbent Length Comment"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmxhead     : num [1:70190] NA NA NA NA NA NA NA NA NA NA ...
##   ..- attr(*, "label")= chr "Head Circumference (cm)"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmihead     : num [1:70190] NA NA NA NA NA NA NA NA NA NA ...
##   ..- attr(*, "label")= chr "Head Circumference Comment"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmxht       : num [1:70190] NA 152 168 NA 156 ...
##   ..- attr(*, "label")= chr "Standing Height (cm)"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmiht       : num [1:70190] NA NA NA NA NA NA NA NA NA NA ...
##   ..- attr(*, "label")= chr "Standing Height Comment"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmxbmi      : num [1:70190] NA 17.4 26.5 NA 30.9 ...
##   ..- attr(*, "label")= chr "Body Mass Index (kg/m**2)"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmxleg      : num [1:70190] NA 37.6 42.7 NA 38 40.4 39.2 41.5 NA NA ...
##   ..- attr(*, "label")= chr "Upper Leg Length (cm)"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmileg      : num [1:70190] NA NA NA NA NA NA NA NA NA NA ...
##   ..- attr(*, "label")= chr "Upper Leg Length Comment"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmxcalf     : num [1:70190] NA 29.3 40.6 NA 36.6 35.6 31.7 42.6 NA NA ...
##   ..- attr(*, "label")= chr "Maximal Calf Circumference (cm)"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmicalf     : num [1:70190] NA NA NA NA NA NA NA NA NA NA ...
##   ..- attr(*, "label")= chr "Maximal Calf Comment"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmxarml     : num [1:70190] 16 34.3 36.5 NA 35 37.5 35.8 42.7 15.6 NA ...
##   ..- attr(*, "label")= chr "Upper Arm Length (cm)"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmiarml     : num [1:70190] NA NA NA NA NA NA NA NA NA NA ...
##   ..- attr(*, "label")= chr "Upper Arm Length Comment"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmxarmc     : num [1:70190] 15.5 21.7 32.6 NA 35.8 31.2 22.3 33 16.7 NA ...
##   ..- attr(*, "label")= chr "Arm Circumference (cm)"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmiarmc     : num [1:70190] NA NA NA NA NA NA NA NA NA NA ...
##   ..- attr(*, "label")= chr "Arm Circumference Comment"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmxwaist    : num [1:70190] NA 62.8 97.8 NA 96 ...
##   ..- attr(*, "label")= chr "Waist Circumference (cm)"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmiwaist    : num [1:70190] NA NA NA NA NA NA NA NA NA NA ...
##   ..- attr(*, "label")= chr "Waist Circumference Comment"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmxthicr    : num [1:70190] NA 39.5 55.9 NA 53.7 48 41.3 50.5 NA NA ...
##   ..- attr(*, "label")= chr "Thigh Circumference (cm)"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmithicr    : num [1:70190] NA NA NA NA NA NA NA NA NA NA ...
##   ..- attr(*, "label")= chr "Thigh Circumference Comment"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmxtri      : num [1:70190] 12.8 10.4 18.8 NA NA 10.4 6.8 13 10.4 NA ...
##   ..- attr(*, "label")= chr "Triceps Skinfold (mm)"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmitri      : num [1:70190] NA NA NA NA 2 NA NA NA NA NA ...
##   ..- attr(*, "label")= chr "Triceps Skinfold Comment"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmxsub      : num [1:70190] 10 8.4 17.6 NA NA 22.2 8.2 16 10.4 NA ...
##   ..- attr(*, "label")= chr "Subscapular Skinfold (mm)"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bmisub      : num [1:70190] NA NA NA NA 1 NA NA NA NA NA ...
##   ..- attr(*, "label")= chr "Subscapular Skinfold Comment"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bpq020      : num [1:70190] NA NA NA 2 1 1 2 1 NA 2 ...
##   ..- attr(*, "label")= chr "Ever told you had high blood pressure"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bpq030      : num [1:70190] NA NA NA NA 1 1 NA 1 NA NA ...
##   ..- attr(*, "label")= chr "Told had high blood pressure - 2+ times"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bpq040a     : num [1:70190] NA NA NA NA 1 1 NA 1 NA NA ...
##   ..- attr(*, "label")= chr "Taking prescription for hypertension"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bpq050a     : num [1:70190] NA NA NA NA 1 1 NA 1 NA NA ...
##   ..- attr(*, "label")= chr "Now taking prescribed medicine for HBP"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bpq060      : num [1:70190] NA NA NA 1 1 1 NA 1 NA 1 ...
##   ..- attr(*, "label")= chr "Ever had blood cholesterol checked"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bpq070      : num [1:70190] NA NA NA 9 1 1 NA 1 NA 1 ...
##   ..- attr(*, "label")= chr "When blood cholesterol last checked"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bpq080      : num [1:70190] NA NA NA 2 2 1 NA 1 NA 2 ...
##   ..- attr(*, "label")= chr "Doctor told you - high cholesterol level"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bpq090a     : num [1:70190] NA NA NA NA NA 1 NA 1 NA NA ...
##   ..- attr(*, "label")= chr "Told to eat less fat for cholesterol"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bpq090b     : num [1:70190] NA NA NA NA NA 1 NA 1 NA NA ...
##   ..- attr(*, "label")= chr "Told to reduce weight for cholesterol"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bpq090c     : num [1:70190] NA NA NA NA NA 1 NA 1 NA NA ...
##   ..- attr(*, "label")= chr "Told to exercise more for cholesterol"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bpq090d     : num [1:70190] NA NA NA NA NA 1 NA 1 NA NA ...
##   ..- attr(*, "label")= chr "Told to take prescriptn for cholesterol"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bpq100a     : num [1:70190] NA NA NA NA NA 1 NA 1 NA NA ...
##   ..- attr(*, "label")= chr "Now eating fewer high fat foods"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bpq100b     : num [1:70190] NA NA NA NA NA 1 NA 1 NA NA ...
##   ..- attr(*, "label")= chr "Now controlling weight"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bpq100c     : num [1:70190] NA NA NA NA NA 1 NA 1 NA NA ...
##   ..- attr(*, "label")= chr "Now increasing exercise"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ bpq100d     : num [1:70190] NA NA NA NA NA 1 NA 1 NA NA ...
##   ..- attr(*, "label")= chr "Now taking prescribed medicine"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ cdq001      : num [1:70190] NA NA NA 2 2 2 NA 2 NA 1 ...
##   ..- attr(*, "label")= chr "SP ever had pain or discomfort in chest"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ cdq002      : num [1:70190] NA NA NA NA NA NA NA NA NA 2 ...
##   ..- attr(*, "label")= chr "SP get it walking uphill or in a hurry"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##  $ cdq003      : num [1:70190] NA NA NA NA NA NA NA NA NA NA ...
##   ..- attr(*, "label")= chr "During an ordinary pace on level ground"
##   ..- attr(*, "format.stata")= chr "%10.0g"
##   [list output truncated]

Slide with Plot