Skip to content

Commit 8fa4b5f

Browse files
authored
Merge pull request #93 from bd2kccd/v1.2.0
Update README.md
2 parents 87e627b + 4a15753 commit 8fa4b5f

File tree

1 file changed

+22
-6
lines changed

1 file changed

+22
-6
lines changed

README.md

Lines changed: 22 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -27,26 +27,42 @@ library(rcausal)
2727
data("charity") #Load the charity dataset
2828

2929
tetradrunner.getAlgorithmDescription(algoId = 'fges')
30-
tetradrunner.getAlgorithmParameters(algoId = 'fges',scoreId = 'fisher-z')
3130
#Compute FGES search
32-
tetradrunner <- tetradrunner(algoId = 'fges',df = charity,scoreId = 'fisher-z',
33-
dataType = 'continuous',alpha=0.1,faithfulnessAssumed=TRUE,maxDegree=-1,verbose=TRUE)
31+
tetradrunner <- tetradrunner(algoId = 'fges',df = charity,scoreId = 'sem-bic',
32+
dataType = 'continuous',faithfulnessAssumed=TRUE,maxDegree=-1,verbose=TRUE)
3433

3534
tetradrunner$nodes #Show the result's nodes
3635
tetradrunner$edges #Show the result's edges
36+
37+
graph <- tetradrunner$graph
38+
graph$getAttribute('BIC')
39+
40+
nodes <- graph$getNodes()
41+
for(i in 0:as.integer(nodes$size()-1)){
42+
node <- nodes$get(i)
43+
cat(node$getName(),": ",node$getAttribute('BIC'),"\n")
44+
}
3745
```
3846
### Discrete Dataset
3947
```R
4048
library(rcausal)
4149
data("audiology") #Load the charity dataset
4250

43-
tetradrunner.getAlgorithmParameters(algoId = 'fges',scoreId = 'bdeu')
4451
#Compute FGES search
45-
tetradrunner <- tetradrunner(algoId = 'fges',df = audiology,scoreId = 'bdeu',dataType = 'discrete',
46-
alpha=0.1,faithfulnessAssumed=TRUE,maxDegree=-1,verbose=TRUE)
52+
tetradrunner <- tetradrunner(algoId = 'fges',df = audiology,scoreId = 'cg-bic-score',dataType = 'discrete',
53+
faithfulnessAssumed=TRUE,maxDegree=-1,verbose=TRUE)
4754

4855
tetradrunner$nodes #Show the result's nodes
4956
tetradrunner$edges #Show the result's edges
57+
58+
graph <- tetradrunner$graph
59+
graph$getAttribute('BIC')
60+
61+
nodes <- graph$getNodes()
62+
for(i in 0:as.integer(nodes$size()-1)){
63+
node <- nodes$get(i)
64+
cat(node$getName(),": ",node$getAttribute('BIC'),"\n")
65+
}
5066
```
5167

5268
### Prior Knowledge

0 commit comments

Comments
 (0)