I have EPIC data from 8 samples.

Sample sheet: Sorry I am new to posting here, so please forgive my bad formatting:

Plate   Bead chip # Well Position   Well Position   Sentrix_ID  Array   Complete Barcode    Sample_Name
1447    1   A01 01A 201172520042    R01C01  201172520042_R01C01 MCF10A N 2d
1447    1   B01 01B 201172520042    R02C01  201172520042_R02C01 MCF10A H 0.5% 2d
1447    1   C01 01C 201172520042    R03C01  201172520042_R03C01 68 N 2d
1447    1   D01 01D 201172520042    R04C01  201172520042_R04C01 68 H 4% 2d
1447    1   E01 01E 201172520042    R05C01  201172520042_R05C01 68 H 0.5% 2d
1447    1   F01 01F 201172520042    R06C01  201172520042_R06C01 68 N 8 wk
1447    1   G01 01G 201172520042    R07C01  201172520042_R07C01 68 N 3d
1447    1   H01 01H 201172520042    R08C01  201172520042_R08C01 68 H 4% 8wk

I am trying to find DMPs from the samples N 8 wk and H 4% 8wk. I did the following for normalizing and processing:

setwd("~/Desktop/HYPNOR/")

baseDir <- file.path("~/Desktop/HYPNOR/")
list.files(file.path(baseDir, "201172520042"))
list.files(file.path(baseDir))

targets <- read.metharray.sheet(baseDir)

RGset <- read.metharray.exp(targets = targets)
RGset
phenoData <- pData(RGset)
phenoData[,1:6]

pd <- pData(RGset)
pd
pd[,1:4]

qcReport(RGset, sampNames = pd$Sample_Name, pdf = "qcReport.pdf")

#Preprocessing/Normalizing/DYE correction
noob <- preprocessNoob(RGset)
RGset
GMset_map <- mapToGenome(noob)
snp_add <- addSnpInfo(GMset_map)

remove_snps <- dropLociWithSnps(snp_add, snps=c("SBE","CpG"), maf=0.05)

dmpFinder(remove_snps, pheno = pd$Sample_Group, type = c("categorical"))

What I have tried hasn't worked and I don't know how to single those two samples only for dmpfinder. I made a new sample sheet with just the two samples, but that didn't work and also created a subset of the data, and that didn't work.

I was wondering if there is a better solution for this? I am quite new to EPIC methylation analysis, so my apologies if I sound naive.

Any help would be awesome, please. Thank you.