Creating a simple SO with PopulationEstimate

This section will walk you through the creation of a simple SO file.

First we create the data structure with the maximum likelihood parameter estimates.

MLE <- data.frame(POP_CL=2.23, POP_V=40, IVCL=0.4, IVV=0.3)

Next we add the additional columnType attribute.

attributes(MLE)$columnType <- list(c("structuralParameter"), c("structuralParameter"), c("varParameter", "variance"), c("varParameter", "variance"))

This sets POP_CL and POP_V to be structural parameters and IVCL and IVV to be variability parameters on the variance scale.

The next step is to create the PopulationEstimates object.

library(libsoc)
popest <- so_PopulationEstimates$new()

The MLE data.frame can now be attached.

popest$MLE <- MLE

Now we create the parent objects of the PopulationEstimates.

est <- so_Estimation$new()
est$PopulationEstimates <- popest
block <- so_SOBlock$new()
block$blkId <- "myEstimationRun"
block$Estimation <- est
so <- so_SO$new()
so$add_SOBlock(block)

It is possible to add messages using a convenient method

block$add_message("WARNING", "myTool", "name_of_warning", "This is the description", 2)

Another method adds a rawresults datafile from the tool.

block$add_rawresults_datafile("Description", "file", "id1")

The resulting structure can be inspected

so
## Reference class object of class "so_SO"
## Field "id":
## NULL
## Field "metadataFile":
## NULL
## Field "Description":
## NULL
## Field "PharmMLRef":
## NULL
## Field "SOBlock":
## [[1]]
## Reference class object of class "so_SOBlock"
## Field "blkId":
## [1] "myEstimationRun"
## Field "ToolSettings":
## NULL
## Field "RawResults":
## Reference class object of class "so_RawResults"
## Field "DataFile":
## [[1]]
## data frame with 0 columns and 0 rows
## 
## Field "GraphicsFile":
## NULL
## Field ".cobj":
## <pointer: 0x56123b65cca0>
## Field "TaskInformation":
## Reference class object of class "so_TaskInformation"
## Field "Message":
## [[1]]
## Reference class object of class "so_Message"
## Field "type":
## [1] "WARNING"
## Field "Toolname":
## [1] "myTool"
## Field "Name":
## [1] "name_of_warning"
## Field "Content":
## [1] "This is the description"
## Field "Severity":
## [1] 2
## Field ".cobj":
## <pointer: 0x56123be38a80>
## 
## Field "OutputFilePath":
## NULL
## Field "RunTime":
## NULL
## Field "NumberChains":
## NULL
## Field "NumberIterations":
## NULL
## Field ".cobj":
## <pointer: 0x56123ad58bc0>
## Field "Estimation":
## Reference class object of class "so_Estimation"
## Field "PopulationEstimates":
## Reference class object of class "so_PopulationEstimates"
## Field "MLE":
##   POP_CL POP_V IVCL IVV
## 1   2.23    40  0.4 0.3
## Field "Bayesian":
## NULL
## Field "OtherMethod":
## NULL
## Field ".cobj":
## <pointer: 0x56123c1ce6b0>
## Field "PrecisionPopulationEstimates":
## NULL
## Field "IndividualEstimates":
## NULL
## Field "PrecisionIndividualEstimates":
## NULL
## Field "Residuals":
## NULL
## Field "Predictions":
## NULL
## Field "OFMeasures":
## NULL
## Field "TargetToolMessages":
## NULL
## Field ".cobj":
## <pointer: 0x56123b1adec0>
## Field "Simulation":
## NULL
## Field "ModelDiagnostic":
## NULL
## Field "OptimalDesign":
## NULL
## Field ".cobj":
## <pointer: 0x56123beaafe0>
## 
## Field ".cobj":
## <pointer: 0x56123af0a2b0>

And written using so$so_write(“filename.SO.xml”)