SUBAYNO

knitr::opts_chunk$set(echo = TRUE)

setwd("D:/ATENIN")

library(stats)
library(dplyr)

## 
## Attaching package: 'dplyr'

## The following objects are masked from 'package:stats':
## 
##     filter, lag

## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union

library(rmarkdown)
library(tidyverse)

## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ forcats   1.0.0     ✔ readr     2.1.4
## ✔ ggplot2   3.4.1     ✔ stringr   1.5.0
## ✔ lubridate 1.9.2     ✔ tibble    3.1.8
## ✔ purrr     1.0.1     ✔ tidyr     1.3.0

## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## ℹ Use the ]8;;http://conflicted.r-lib.org/conflicted package]8;; to force all conflicts to become errors

library(readxl)
library(rstatix)

## 
## Attaching package: 'rstatix'
## 
## The following object is masked from 'package:stats':
## 
##     filter

library(car)

## Loading required package: carData
## 
## Attaching package: 'car'
## 
## The following object is masked from 'package:purrr':
## 
##     some
## 
## The following object is masked from 'package:dplyr':
## 
##     recode

library(ggpubr)

library(readxl)
DATA<- read_excel("D:/ATENIN/Book1.xlsx")

a <- DATA%>%
  group_by(Age)%>%
  summarise(count=n())%>%
  mutate(Percentage =round((count/sum(count)*100),2))
paged_table(a)

b <- DATA%>%
  group_by(`Nature of Work`)%>%
  summarise(count=n())%>%
  mutate(Percentage =round((count/sum(count)*100),2))
paged_table(b)

c <- DATA%>%
  group_by(`Years in Service`)%>%
  summarise(count=n())%>%
  mutate(Percentage =round((count/sum(count)*100),2))
paged_table(c)

Dat2<-DATA%>%
  mutate(Scale = ifelse(`Average` <= 1.8, "Strongly Disagree",
                                     ifelse(`Average` <= 2.6, "Disagree",
                                            ifelse(`Average` <= 3.4, "Neutral",
                                                   ifelse(`Average` <=4.2, "Agree", "Strongly Disagree")))))%>%
  group_by(Average)
Dat2

    Shapiro-Wilk normality test

data:  Dat2$Average
W = 0.92805, p-value = 0.004665

Data4<-Dat2%>%
  group_by(Age == "28-35 y.o", Average)%>%
  filter(`Age == "28-35 y.o"`== "TRUE")
Data4

Data3<-Dat2%>%
  group_by(Age == "18-27 y.o", Average)%>%
  filter(`Age == "18-27 y.o"`== "TRUE")
Data3

median(Data3$Average)

## [1] 3.7

median(Data4$Average)

## [1] 3.85

Data5<-Dat2%>%
  group_by(Age == "36-45 y.o", Average)%>%
  filter(`Age == "36-45 y.o"`== "TRUE")
Data5

median(Data5$Average)

## [1] 3.675

Data6<-Dat2%>%
  group_by(Age == "46 y.o & above", Average)%>%
  filter(`Age == "46 y.o & above"`== "TRUE")
Data6

median(Data3$Average)

## [1] 3.7

median(Data4$Average)

## [1] 3.85

median(Data5$Average)

## [1] 3.675

median(Data6$Average)

## [1] 3.375