From 8f1c6201bb5cf6f1c0432bb580ee6422e468030c Mon Sep 17 00:00:00 2001 From: Efrain Gonzalez Date: Fri, 25 Aug 2017 11:27:09 -0400 Subject: [PATCH] Updated Version --- RMarkovBlanket.r | 30 +++++++++++++++++++++--------- 1 file changed, 21 insertions(+), 9 deletions(-) diff --git a/RMarkovBlanket.r b/RMarkovBlanket.r index d601bce..475422c 100644 --- a/RMarkovBlanket.r +++ b/RMarkovBlanket.r @@ -1,5 +1,5 @@ #Efrain Gonzalez -#7/25/2017 +#8/25/2017 #Code for Markov Blanket @@ -91,7 +91,7 @@ for(j in 1:sizeDotP1){ } lrgMarkov <- dim(NewDotP2_2)[1] -MBlanky <- function(MarkovDegree = lrgMarkov, VariableStartName = "Alzheimer", VariableEndName = "Age"){ +Blanky <- function(MarkovDegree = lrgMarkov, VariableStartName = "Alzheimer", VariableEndName = "Age"){ #Finding the Parents and Children d <- 1 AllNamList1 <- vector("list",length = 3) @@ -124,9 +124,12 @@ MBlanky <- function(MarkovDegree = lrgMarkov, VariableStartName = "Alzheimer", V nc <- 1 for(nc in 1:NumofChild){ LocCOPofVar <- grep(CofVar[nc],NewDotP2_2[,2]) - COPofVar <- NewDotP2_2[LocCOPofVar,1] - if(grepl(COPofVar,VariableStartName)){ - next + COPofVar <- NewDotP2_2[LocCOPofVar,1] + if(sum(grepl(VariableStartName,COPofVar)) >= 1){ + #positions of variable start name within the vector of co parents + posoforig <- grep(VariableStartName,COPofVar) + COPofVar <- COPofVar[-posoforig] + COPlist <- append(COPlist,COPofVar) } else{ #COPlist[[nc]] <- COPofVar[COPofVar!=VariableStartName] COPlist <- append(COPlist,COPofVar) @@ -180,8 +183,11 @@ MBlanky <- function(MarkovDegree = lrgMarkov, VariableStartName = "Alzheimer", V for(np in 1:NumofVars){ LocCofVar <- grep(AllVarList[[d-1]][[ef]][np],NewDotP2_2[,1]) CofVar <- NewDotP2_2[LocCofVar,2] - #if(grepl(VariableStartName,CofVar)){ - # next + #if(sum(grepl(VariableStartName,CofVar)) >= 1){ + # #positions of variable start name within the vector of co parents + # posoforig <- grep(VariableStartName,COPofVar) + # COPofVar <- COPofVar[-posoforig] + # COPlist <- append(COPlist,COPofVar) #} else{ CofVlist <- append(CofVlist,CofVar) #} @@ -202,7 +208,14 @@ MBlanky <- function(MarkovDegree = lrgMarkov, VariableStartName = "Alzheimer", V for(ncp in 1:NumofCVars){ LocCPofCVar <- grep(CofVlist[ncp],NewDotP2_2[,2]) CPofCVar <- NewDotP2_2[LocCPofCVar,1] - CPofClist <- append(CPofClist,CPofCVar) + #if(sum(grepl(,CPofCVar)) >= 1){ + # #positions of variable start name within the vector of co parents + # posoforig <- grep(VariableStartName,COPofVar) + # COPofVar <- COPofVar[-posoforig] + # COPlist <- append(COPlist,COPofVar) + #} else{ + CPofClist <- append(CPofClist,CPofCVar) + #} ncp <- ncp + 1 } @@ -224,4 +237,3 @@ MBlanky <- function(MarkovDegree = lrgMarkov, VariableStartName = "Alzheimer", V ##The Markov Degree is that found below d } -#Now use the command MBlanky() with the appropriate settings \ No newline at end of file -- 2.29.0