Final.knit

---
title: "Analysis of PMN and Treatment Practices"
output: html_document
---

## Setup and Libraries




``` r
# Import the dataset
df <- read.csv('data.csv')

# View the first few rows of the dataset
head(df)

##   No_AI No_NI AI_NI Breed Age Parity BCS Feed_Supply HUD Pract_treatment PMN
## 1     4     2     1     2   3      1   4           2   2               1   4
## 2     6     1     1     2   7      3   4           3   2               1   7
## 3     7     1     1     2   5      1   4           3   2               1  11
## 4     5     1     1     2   2      0   1           3   2               2   7
## 5     5     0     0     2   5      1   2           3   2               5  12
## 6     5     1     1     2   7      3   1           2   2               2   8
##              IMVic Bacterial
## 1 Escherichia coli       Yes
## 2 Escherichia coli       Yes
## 3 Escherichia coli       Yes
## 4 Escherichia coli       Yes
## 5 Escherichia coli       Yes
## 6 Escherichia coli       Yes

# Fix Breed column: replace 'ho' with 2 and convert to numeric
df$Breed[df$Breed == 'ho'] <- 2
df$Breed <- as.numeric(df$Breed)

# Extract numeric values from Age column
df$Age_years <- as.numeric(str_extract(df$Age, "\\d+"))

# Create bacterial infection status based on IMVic column
df <- df %>% 
  mutate(bacterial_infection = ifelse(IMVic == "No isolation", 'No', 'Yes')) %>% 
  mutate(bacterial_infection = factor(bacterial_infection))

# Create PMN category based on median
median_pmn <- median(df$PMN, na.rm = TRUE)
df <- df %>% 
  mutate(PMN_category = factor(ifelse(PMN >= median_pmn, "High", "Low")))

# Create Age Groups
df <- df %>% 
  mutate(Age_group = cut(Age_years,
                         breaks = c(2, 4, 6, 8, Inf),
                         labels = c('2-3', '4-5', '6-7', '8+')))

# Create Parity Groups
df <- df %>% 
  mutate(Parity_group = case_when(
    Parity %in% c("0", "1") ~ as.character(Parity),
    TRUE ~ "2+"
  )) %>% 
  mutate(Parity_group = factor(Parity_group))

# Convert categorical variables to factors
df <- df %>% 
  mutate(
    BCS = factor(BCS),
    Pract_treatment = factor(Pract_treatment),
    Feed_Supply = factor(Feed_Supply)
  )

# Perform Shapiro-Wilk test for normality on PMN
shapiro.test(df$PMN)  # Result suggests normality

## 
##  Shapiro-Wilk normality test
## 
## data:  df$PMN
## W = 0.96342, p-value = 0.006027

# Perform Chi-Square test for PMN Category vs Bacterial Infection
chi_bact <- chisq.test(table(df$PMN_category, df$bacterial_infection))
print(chi_bact)

## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  table(df$PMN_category, df$bacterial_infection)
## X-squared = 8.6816, df = 1, p-value = 0.003214

# Levene's Test for Homogeneity of Variance: PMN vs Age Group
df <- df %>% filter(!is.na(Age_group))
leveneTest(PMN ~ Age_group, data = df)

## Levene's Test for Homogeneity of Variance (center = median)
##       Df F value Pr(>F)
## group  2  0.3417 0.7115
##       85

# Perform Fisher's Exact Test for PMN vs Parity
fisher_parity <- fisher.test(table(df$PMN_category, df$Parity_group))
print(fisher_parity)

## 
##  Fisher's Exact Test for Count Data
## 
## data:  table(df$PMN_category, df$Parity_group)
## p-value = 1
## alternative hypothesis: two.sided

# Modify Parity Group and perform Chi-Square test again
df$Parity_group <- factor(ifelse(df$Parity_group == "0", "1", as.character(df$Parity_group)))
chi_parity <- chisq.test(table(df$PMN_category, df$Parity_group))
print(chi_parity)

## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  table(df$PMN_category, df$Parity_group)
## X-squared = 0, df = 1, p-value = 1

# Chi-Square Test with Simulation
chi_parity <- chisq.test(table(df$PMN_category, df$Parity_group), simulate.p.value = TRUE, B = 2000)
print(chi_parity)

## 
##  Pearson's Chi-squared test with simulated p-value (based on 2000
##  replicates)
## 
## data:  table(df$PMN_category, df$Parity_group)
## X-squared = 0.0034521, df = NA, p-value = 1

# Kruskal-Wallis test for PMN vs BCS
kw_bcs <- kruskal.test(PMN ~ BCS, data = df)
print(kw_bcs)

## 
##  Kruskal-Wallis rank sum test
## 
## data:  PMN by BCS
## Kruskal-Wallis chi-squared = 11.226, df = 4, p-value = 0.02414

# Perform Chi-Square test for PMN Category vs Practical Treatment
chi_treatment <- chisq.test(table(df$PMN_category, df$Pract_treatment))
print(chi_treatment)

## 
##  Pearson's Chi-squared test
## 
## data:  table(df$PMN_category, df$Pract_treatment)
## X-squared = 8.0563, df = 4, p-value = 0.08954

# Group treatments and perform Chi-Square test
df <- df %>%
  mutate(Pract_treatment_grouped = case_when(
    Pract_treatment %in% c("Uterine Lavage", "Mineral Mixture") ~ "Group 1",
    Pract_treatment %in% c("Hormonal Therapy", "All") ~ "Group 2",
    TRUE ~ "Nothing"
  ))

df$Pract_treatment_grouped <- factor(df$Pract_treatment_grouped)

# Perform Chi-Square Test again for grouped treatments
chi_treatment <- chisq.test(table(df$PMN_category, df$Pract_treatment_grouped))
print(chi_treatment)

## 
##  Chi-squared test for given probabilities
## 
## data:  table(df$PMN_category, df$Pract_treatment_grouped)
## X-squared = 0.18182, df = 1, p-value = 0.6698

# Create binary treatment variable and perform Chi-Square test
df <- df %>%
  mutate(Treatment_binary = ifelse(Pract_treatment == "Nothing", "No Treatment", "Received Treatment")) %>%
  mutate(Treatment_binary = factor(Treatment_binary))

# Perform Chi-Square Test for binary treatment
chi_treatment <- chisq.test(table(df$PMN_category, df$Treatment_binary))
print(chi_treatment)

## 
##  Chi-squared test for given probabilities
## 
## data:  table(df$PMN_category, df$Treatment_binary)
## X-squared = 0.18182, df = 1, p-value = 0.6698

# Fisher's Exact Test for PMN vs Pract_treatment
fisher_treatment <- fisher.test(table(df$PMN_category, df$Pract_treatment))
print(fisher_treatment)

## 
##  Fisher's Exact Test for Count Data
## 
## data:  table(df$PMN_category, df$Pract_treatment)
## p-value = 0.09138
## alternative hypothesis: two.sided

# Fisher's Exact Test with Simulation
chi_treatment <- chisq.test(table(df$PMN_category, df$Pract_treatment), simulate.p.value = TRUE, B = 2000)
print(chi_treatment)

## 
##  Pearson's Chi-squared test with simulated p-value (based on 2000
##  replicates)
## 
## data:  table(df$PMN_category, df$Pract_treatment)
## X-squared = 8.0563, df = NA, p-value = 0.08946

# Loop through treatments and perform Chi-Square tests for each
for (treatment in unique(df$Pract_treatment)) {
  df_subset <- df %>% mutate(Pract_treatment_binary = ifelse(Pract_treatment == treatment, "Yes", "No"))
  
  chi_test <- chisq.test(table(df_subset$PMN_category, df_subset$Pract_treatment_binary))
  
  cat("\nChi-Square Test for:", treatment, "\n")
  print(chi_test)
}

## 
## Chi-Square Test for: 1 
## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  table(df_subset$PMN_category, df_subset$Pract_treatment_binary)
## X-squared = 1.833e-31, df = 1, p-value = 1
## 
## 
## Chi-Square Test for: 5 
## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  table(df_subset$PMN_category, df_subset$Pract_treatment_binary)
## X-squared = 3.1494, df = 1, p-value = 0.07596
## 
## 
## Chi-Square Test for: 2 
## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  table(df_subset$PMN_category, df_subset$Pract_treatment_binary)
## X-squared = 0.27112, df = 1, p-value = 0.6026
## 
## 
## Chi-Square Test for: 3 
## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  table(df_subset$PMN_category, df_subset$Pract_treatment_binary)
## X-squared = 2.7041, df = 1, p-value = 0.1001
## 
## 
## Chi-Square Test for: 4 
## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  table(df_subset$PMN_category, df_subset$Pract_treatment_binary)
## X-squared = 0.39539, df = 1, p-value = 0.5295

# Loop through treatments and perform Fisher's Exact tests for each
for (treatment in unique(df$Pract_treatment)) {
  df_subset <- df %>% mutate(Pract_treatment_binary = ifelse(Pract_treatment == treatment, "Yes", "No"))
  
  fisher_test <- fisher.test(table(df_subset$PMN_category, df_subset$Pract_treatment_binary))
  
  cat("\nFisher's Exact Test for:", treatment, "\n")
  print(fisher_test)
}

## 
## Fisher's Exact Test for: 1 
## 
##  Fisher's Exact Test for Count Data
## 
## data:  table(df_subset$PMN_category, df_subset$Pract_treatment_binary)
## p-value = 1
## alternative hypothesis: true odds ratio is not equal to 1
## 95 percent confidence interval:
##  0.3102923 3.0695599
## sample estimates:
## odds ratio 
##    0.98202 
## 
## 
## Fisher's Exact Test for: 5 
## 
##  Fisher's Exact Test for Count Data
## 
## data:  table(df_subset$PMN_category, df_subset$Pract_treatment_binary)
## p-value = 0.05207
## alternative hypothesis: true odds ratio is not equal to 1
## 95 percent confidence interval:
##  0.02067101 1.10338048
## sample estimates:
## odds ratio 
##   0.208894 
## 
## 
## Fisher's Exact Test for: 2 
## 
##  Fisher's Exact Test for Count Data
## 
## data:  table(df_subset$PMN_category, df_subset$Pract_treatment_binary)
## p-value = 0.498
## alternative hypothesis: true odds ratio is not equal to 1
## 95 percent confidence interval:
##  0.5213948 3.8323739
## sample estimates:
## odds ratio 
##   1.404703 
## 
## 
## Fisher's Exact Test for: 3 
## 
##  Fisher's Exact Test for Count Data
## 
## data:  table(df_subset$PMN_category, df_subset$Pract_treatment_binary)
## p-value = 0.07913
## alternative hypothesis: true odds ratio is not equal to 1
## 95 percent confidence interval:
##   0.8376056 15.4594142
## sample estimates:
## odds ratio 
##   3.237622 
## 
## 
## Fisher's Exact Test for: 4 
## 
##  Fisher's Exact Test for Count Data
## 
## data:  table(df_subset$PMN_category, df_subset$Pract_treatment_binary)
## p-value = 0.4171
## alternative hypothesis: true odds ratio is not equal to 1
## 95 percent confidence interval:
##  0.1620197 2.0635608
## sample estimates:
## odds ratio 
##  0.6034639

# Logistic regression model for PMN category vs Pract_treatment
logit_model <- glm(PMN_category ~ Pract_treatment, data = df, family = binomial)

# Display summary of the logistic regression model
summary(logit_model)

## 
## Call:
## glm(formula = PMN_category ~ Pract_treatment, family = binomial, 
##     data = df)
## 
## Coefficients:
##                  Estimate Std. Error z value Pr(>|z|)
## (Intercept)       -0.1054     0.4595  -0.229    0.819
## Pract_treatment2   0.2485     0.5956   0.417    0.677
## Pract_treatment3   1.0217     0.7491   1.364    0.173
## Pract_treatment4  -0.4055     0.6912  -0.587    0.557
## Pract_treatment5  -1.3987     0.9068  -1.543    0.123
## 
## (Dispersion parameter for binomial family taken to be 1)
## 
##     Null deviance: 121.81  on 87  degrees of freedom
## Residual deviance: 113.31  on 83  degrees of freedom
## AIC: 123.31
## 
## Number of Fisher Scoring iterations: 4

# Group treatments further and perform Chi-Square test
df <- df %>%
  mutate(Pract_treatment_grouped = case_when(
    Pract_treatment %in% c("Uterine Lavage", "Hormonal Therapy") ~ "Treatment A",
    Pract_treatment %in% c("Mineral Mixture", "All") ~ "Treatment B",
    TRUE ~ "No Treatment"
  ))

chi_grouped <- chisq.test(table(df$PMN_category, df$Pract_treatment_grouped))
print(chi_grouped)

## 
##  Chi-squared test for given probabilities
## 
## data:  table(df$PMN_category, df$Pract_treatment_grouped)
## X-squared = 0.18182, df = 1, p-value = 0.6698

# Define function for bootstrap Chi-Square test
boot_func <- function(data, index) {
  sampled_data <- data[index, ]
  chisq.test(table(sampled_data$PMN_category, sampled_data$Pract_treatment))$statistic
}

set.seed(123)
boot_result <- boot(df, boot_func, R = 1000)

## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect
## Warning in chisq.test(table(sampled_data$PMN_category,
## sampled_data$Pract_treatment)): Chi-squared approximation may be incorrect

# Display bootstrap confidence interval
boot.ci(boot_result, type = "perc")

## BOOTSTRAP CONFIDENCE INTERVAL CALCULATIONS
## Based on 1000 bootstrap replicates
## 
## CALL : 
## boot.ci(boot.out = boot_result, type = "perc")
## 
## Intervals : 
## Level     Percentile     
## 95%   ( 2.713, 24.731 )  
## Calculations and Intervals on Original Scale

# Boxplot: PMN Distribution by Bacterial Infection
p1 <- ggplot(df, aes(bacterial_infection, PMN)) +
  geom_boxplot(aes(fill = factor(..middle..)), show.legend = FALSE) +
  labs(title = 'PMN Distribution by Bacterial Infection', x = 'Bacterial Infection', y = 'PMN (numbers)') +
  theme_bw() +
  theme(plot.title = element_text(hjust = 0.5, face = "bold", size = 14, family = "Times")) +
  scale_fill_brewer(palette = 'Blues')
print(p1)

## Warning: The dot-dot notation (`..middle..`) was deprecated in ggplot2 3.4.0.
## ℹ Please use `after_stat(middle)` instead.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.

## Warning in grid.Call(C_stringMetric, as.graphicsAnnot(x$label)): font family
## not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database

## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

# Bar plot: PMN Count by Age Group
p2 <- ggplot(df, aes(x = factor(Age_group), y = PMN, fill = factor(Age_group))) +
  geom_bar(stat = "identity") +
  labs(title = "PMN Count by Age Group", x = "Age Group", y = "PMN Count") +
  scale_fill_viridis_d(option = 'C') +
  theme_minimal()
print(p2)

# Boxplot: PMN vs Parity Group
p3 <- ggplot(df, aes(x = Parity_group, y = PMN)) +
  geom_boxplot(aes(fill = factor(..middle..))) +
  labs(x = "Parity Group", y = "PMN Count") +
  theme_minimal()
print(p3)

# Boxplot: PMN vs Body Condition Score (BCS)
p4 <- ggplot(df, aes(x = BCS, y = PMN)) +
  geom_boxplot(aes(fill = factor(..middle..)), show.legend = FALSE) +
  labs(title = 'Effect of Body Condition Score on PMN', x = "Body Condition Score", y = "PMN Count") +
  theme_minimal()
print(p4)

# Boxplot: PMN vs Treatment Practice
custom_label <- c('1' = 'Uterine Lavage', '2' = 'Mineral Mixture', '3' = 'Hormonal Therapy', '4' = 'All', '5' = 'Nothing')
p5 <- ggplot(df, aes(x = Pract_treatment, y = PMN)) +
  geom_boxplot(aes(fill = factor(..middle..)), show.legend = FALSE) +
  scale_x_discrete(labels = custom_label) +
  labs(title = 'Effect of Treatment Practice on PMN count', x = "Treatment Practice", y = "PMN Count") +
  theme_bw()
print(p5)