pgab

High/moderate intentsity

data_pasieke1.8 <- data[Maziau1.8 == 1] %>% mutate(groupA = ifelse(high ==1,1,0),groupB= ifelse(high ==0,1,0))

# Perform Chi-squared
table <- table(data_pasieke1.8$groupA)


result <- chisq.test(table)

cat("Test statistic = ", result$statistic, "\n")

## Test statistic =  47.432

cat("p-value = ", result$p.value, "\n")

## p-value =  5.694545e-12

# Print the test statistic, p-value and alternative hypothesis
cat("Tiriamuju skaicius: ",nrow(data_pasieke1.8))

## Tiriamuju skaicius:  125

cat("High intensity pacientu: ",nrow(data_pasieke1.8[groupA == 1]),
    ". Procentas: ", nrow(data_pasieke1.8[groupA == 1])*100/nrow(data_pasieke1.8))

## High intensity pacientu:  101 . Procentas:  80.8

cat("Moderate intensity pacientu: ",nrow(data_pasieke1.8[groupA == 0]),
     ". Procentas: ", nrow(data_pasieke1.8[groupA == 0])*100/nrow(data_pasieke1.8)    )

## Moderate intensity pacientu:  24 . Procentas:  19.2

Tie jibani vaistai

Chi squared testas naudojas pratestuoti frequencies tarp dvieju variables

 duomenys <- data %>%  mutate(
                                 
                                 Simvastatinas_c =ifelse(!is.na(Simvastatinas), 1,0),
                                Rozuvastatinas_c =ifelse(!is.na(Rozuvastatinas), 1,0),
                                Atorvastatinas_c =ifelse(!is.na(Atorvastatinas), 1,0))
 
ator<-table(duomenys[Atorvastatinas_c == 1]$Maziau1.8)
rozu<-table(duomenys[Rozuvastatinas_c == 1]$Maziau1.8)
cont_table <- rbind(ator, rozu)

result <- chisq.test(cont_table)

cat("Test statistic = ", result$statistic, "\n")

## Test statistic =  0.749384

cat("p-value = ", result$p.value, "\n")

## p-value =  0.3866713

print("kiek pacientu ir koks procentas pasieke MTL<1.4 su Rozuvastatino vaistu")

## [1] "kiek pacientu ir koks procentas pasieke MTL<1.4 su Rozuvastatino vaistu"

nrow(duomenys[Maziau1.8 == 1 & Atorvastatinas_c == 1])

## [1] 102

nrow(duomenys[Maziau1.8 == 1 & Atorvastatinas_c == 1])*100/nrow(duomenys[Atorvastatinas_c == 1])

## [1] 24.75728

print("kiek pacientu ir koks procentas pasieke MTL<1.4 su Atorvastatino vaistu")

## [1] "kiek pacientu ir koks procentas pasieke MTL<1.4 su Atorvastatino vaistu"

nrow(duomenys[Maziau1.8 == 1 & Rozuvastatinas_c == 1])

## [1] 18

nrow(duomenys[Maziau1.8 == 1 & Rozuvastatinas_c == 1])*100/nrow(duomenys[Rozuvastatinas_c == 1])

## [1] 31.03448