Aids Data Analysis

Improting Data Set

data1<-read.csv("D:/Amoud University/Semester 3/APS 8115 Statistical Programming & Data Mgt/aids_data.csv")
age<-data1$age

Descriptive Statistics

library(EnvStats)
## Warning: package 'EnvStats' was built under R version 4.1.3
## 
## Attaching package: 'EnvStats'
## The following objects are masked from 'package:stats':
## 
##     predict, predict.lm
## The following object is masked from 'package:base':
## 
##     print.default
library(moments)
## Warning: package 'moments' was built under R version 4.1.3
## 
## Attaching package: 'moments'
## The following objects are masked from 'package:EnvStats':
## 
##     kurtosis, skewness
library(AdequacyModel)
library(psycho)
## Warning: package 'psycho' was built under R version 4.1.3
mean(age)
## [1] 3.503608
median(age)
## [1] 3.610918
mean(age,0.30)
## [1] 3.571432
weighted.mean(age)
## [1] 3.503608
geoMean(age)
## [1] 3.457378
var(age)
## [1] 0.2891845
sd(age)
## [1] 0.5377588
skewness(age)
## [1] -0.7327695
kurtosis(age)
## [1] 3.342799
moment(age,order=3)
## [1] 45.91942
moment(age,order=5)
## [1] 642.345
cv(age)
## [1] 0.1534872
geoSD(age)
## [1] 1.184711
range(age)
## [1] 1.386294 4.382027
moment(age,order=2)
## [1] 12.56305
summaryFull(age)
##                                   age
## N                            206.0000
## Mean                           3.5040
## Median                         3.6110
## 10% Trimmed Mean               3.5460
## Geometric Mean                 3.4570
## Skew                          -0.7382
## Kurtosis                       0.3810
## Min                            1.3860
## Max                            4.3820
## Range                          2.9960
## 1st Quartile                   3.2190
## 3rd Quartile                   3.8920
## Standard Deviation             0.5378
## Geometric Standard Deviation   1.1850
## Interquartile Range            0.6730
## Median Absolute Deviation      0.5231
## Coefficient of Variation       0.1535
## attr(,"class")
## [1] "summaryStats"
## attr(,"stats.in.rows")
## [1] TRUE
## attr(,"drop0trailing")
## [1] TRUE
descriptive(age)

## $mean
## [1] 3.50361
## 
## $median
## [1] 3.61092
## 
## $mode
## [1] 3.75
## 
## $variance
## [1] 0.28918
## 
## $Skewness
## [1] -0.73277
## 
## $Kurtosis
## [1] 0.3428
## 
## $minimum
## [1] 1.38629
## 
## $maximum
## [1] 4.38203
## 
## $n
## [1] 206

Data Visualisation

Histogram

hist(age,col="orange",main="Age in Aids Data",
     xlab="Age in Years",ylab="Cumulative Frequency")

Box Plot

boxplot(age,col="light green",main="Age in Aids Data",
     ylab="Age in Years")

VIOPLOT

library(vioplot)
## Warning: package 'vioplot' was built under R version 4.1.3
## Loading required package: sm
## Warning: package 'sm' was built under R version 4.1.3
## Package 'sm', version 2.2-5.7: type help(sm) for summary information
## Loading required package: zoo
## Warning: package 'zoo' was built under R version 4.1.3
## 
## Attaching package: 'zoo'
## The following objects are masked from 'package:base':
## 
##     as.Date, as.Date.numeric
vioplot(age,col="light Blue",main="Age in Aids Data",
        ylab="Age in Years")

Kernel Density

plot(density(age),col="Red",main="Age in Aids Data",
     xlab="Age in Years",ylab="Cumulative Frequency")
polygon(density(age),col="Yellow",border="black",main="Age in Aids Data",
        xlab="Age in Years",ylab="Cumulative Frequency")