library("UsingR")
## Warning: package 'UsingR' was built under R version 4.0.5
## Loading required package: MASS
## Loading required package: HistData
## Warning: package 'HistData' was built under R version 4.0.5
## Loading required package: Hmisc
## Warning: package 'Hmisc' was built under R version 4.0.5
## Loading required package: lattice
## Loading required package: survival
## Loading required package: Formula
## Loading required package: ggplot2
## Warning: package 'ggplot2' was built under R version 4.0.5
## 
## Attaching package: 'Hmisc'
## The following objects are masked from 'package:base':
## 
##     format.pval, units
## 
## Attaching package: 'UsingR'
## The following object is masked from 'package:survival':
## 
##     cancer
reaction.time
##      age gender control     time
## 1  16-24      F       T 1.360075
## 2  16-24      M       T 1.467939
## 3  16-24      M       T 1.512036
## 4  16-24      F       T 1.390647
## 5  16-24      M       T 1.384208
## 6  16-24      M       C 1.393875
## 7  16-24      M       T 1.419043
## 8  16-24      F       T 1.461187
## 9  16-24      F       T 1.382461
## 10 16-24      M       C 1.260152
## 11 16-24      M       T 1.292917
## 12 16-24      F       T 1.331092
## 13 16-24      M       T 1.369171
## 14 16-24      F       C 1.330471
## 15 16-24      F       T 1.344506
## 16 16-24      M       C 1.357869
## 17 16-24      F       T 1.520244
## 18 16-24      M       T 1.352397
## 19 16-24      F       T 1.348420
## 20 16-24      F       T 1.373000
## 21   25+      M       C 1.306813
## 22   25+      F       C 1.384192
## 23   25+      M       T 1.345643
## 24   25+      M       C 1.594901
## 25   25+      M       T 1.341972
## 26   25+      F       T 1.472212
## 27   25+      F       T 1.440685
## 28   25+      M       T 1.502757
## 29   25+      F       T 1.501541
## 30   25+      F       T 1.562813
## 31   25+      F       C 1.475851
## 32   25+      F       T 1.475456
## 33   25+      M       C 1.312376
## 34   25+      F       T 1.447160
## 35   25+      M       T 1.463118
## 36   25+      M       C 1.293425
## 37   25+      F       C 1.402860
## 38   25+      M       C 1.245497
## 39   25+      F       T 1.475806
## 40   25+      M       T 1.520106
## 41   25+      F       T 1.437079
## 42   25+      F       C 1.583610
## 43   25+      F       T 1.499246
## 44   25+      M       C 1.436701
## 45   25+      F       T 1.524056
## 46   25+      F       T 1.510929
## 47   25+      F       C 1.522563
## 48   25+      M       C 1.462412
## 49   25+      F       C 1.313935
## 50   25+      F       T 1.440382
## 51   25+      M       T 1.444503
## 52   25+      F       C 1.316956
## 53   25+      M       T 1.517435
## 54   25+      M       T 1.536446
## 55   25+      M       T 1.444865
## 56   25+      M       C 1.442599
## 57   25+      M       C 1.355206
## 58   25+      F       T 1.614979
## 59   25+      M       T 1.527699
## 60   25+      M       T 1.466591
total_MT=sum(reaction.time$time[reaction.time$gender == "M" & reaction.time$control == "T"])
total_MT
## [1] 25.90884
total_MC=sum(reaction.time$time[reaction.time$gender == "M" & reaction.time$control == "C"])
total_MC
## [1] 16.46183
total_FT=sum(reaction.time$time[reaction.time$gender == "F" & reaction.time$control == "T"])
total_FT
## [1] 31.91398
total_FC=sum(reaction.time$time[reaction.time$gender == "F" & reaction.time$control == "C"])
total_FC
## [1] 11.33044

##segunda##

barplot(table(reaction.time$gender, reaction.time$control), main = "Relacion entre genero y control ")

mayor_MT=max(reaction.time$time[reaction.time$gender == "M" & reaction.time$control == "T"])
mayor_MT
## [1] 1.536446
mayor_MC=max(reaction.time$time[reaction.time$gender == "M" & reaction.time$control == "C"])
mayor_MC
## [1] 1.594901
mayor_FT=max(reaction.time$time[reaction.time$gender == "F" & reaction.time$control == "T"])
mayor_FT
## [1] 1.614979
mayor_FC=max(reaction.time$time[reaction.time$gender == "F" & reaction.time$control == "C"])
mayor_FC
## [1] 1.58361

##cuarta##

barplot(table(reaction.time$gender, reaction.time$control), main = "Cantidad de muestras por año", ylab="No. Gender", xlab="control")