One function to run them all… Or just eval

阳光穿透心脏的1/2处 2022-07-24 07:14 164阅读 0赞

When I was writing a [sysBio][] package, I needed a function that can be customized based on the user specified model. Such function would include user specified variables (parameters, rates, etc.) with assigned values, functions/expression that use those variables, and the ordinary differential equations that describe the model. After doing some reading, the *eval* function seemed like a logical choice. It was pretty easy to figure out how to use it and to implement the custom function I needed (if interested, you can see it in [makeODEs.R][] and [simulateModel.R][]).

This made me think about how easy it would be to use *eval* to create a general purpose function. Before I go any further, let me first say that I am aware of issues of using *eval* (security, harder readability and debugging, extra complexity, etc). If the hypothetical situation described below is a real-life problem I am working on, I would weight pros and cons before calling the presented *eval* approach the chosen approach. But this is just me being curious and we can pretend it is ok.

So, let’s say that we are working with people (users) whose goals are to create reports from different real-life data. Let’s also assume that these data are represented in tabular format and that reports represent summaries based on mathematical functions applied on the data presented in the tables. Additionally, let’s assume that our users have phobia from writing or learning any type of code (I have to add this, because my significant other had a very strong stand against *eval* and advocated for teaching the users to learn or even copy/paste the code in instead).

Clearly, using R, we could address users’ needs pretty easy – we could load and store users’ data in data frames, (check the data for consistency), use some type of *ply* function to summarize the data, and print the reports in the output files. Ignoring the read/check/print parts, we could create a summary in a single line command.

Now, let’s assume that users often want to change their summaries, for example every few months they remove one of the old summary functions and add one new . In order to keep users’ demands current, we would have to go and change our summary commands quite often.

This is where *eval* approach comes to play. Can we create a single, general function that will be able to handle different, specific users’ data summary requests?

The answer is yes.  
However, note that the presented approach does not do any data evaluation and checking, nor it performs security check on user passed strings (which, in real-life, would be pretty big security concern).

We will assume that user’s data is already in the form of data frame (with *stringsAsFactors* set to *FALSE*). Furthermore, we will assume that users provide a list of variables (that correspond to table column names) over which they want to split and summarize the data (character vector), as well as a two-column table (data frame) which will contain the user specified name of the function used for summary (first column) and the function used for summary (second column).

As data example, we can use *mtcars* data set from the [datasets][]library. We will assume that the user is interested in information about gears and for each gear wants to calculate some custom value (defined as *mpg+2\*carb*, total hp (defined as *sum(hp)*), and average mpg (defined as *mean(mpg)*):

`ex1_df <- datasets::mtcars` 
`ex1_vars <- "gear"` 
`ex1_formulas <- data.frame(labs=c("cust_val", "total_hp", "avg_mpg"), ress=c("mpg+2*carb", "sum(hp)", "mean(mpg)"), stringsAsFactors = FALSE)`

If we were coding the user's request, we could use the ddply function from the [plyr][] package as follow: 
`ex1_trueResult <- plyr::ddply(ex1_df, .variable=c("gear"), function(x) data.frame(cust_val = x$mpg+2*x$carb, total_hp = sum(x$hp), avg_mpg = mean(x$mpg)))`

The first thing I tried to do was to create a function that will recreate the *ddply* call:

[view raw][] [oneForMany\_first.R][oneForMany_first.R] hosted with ❤ by  [GitHub][]

Running this function will create a summary, but only for avg\_mpg: 
`ex1_res01 <- oneForMany_first(ex1_df, ex1_vars, ex1_formulas)`

This is because when an expression vector is passed to eval, the result of evaluating is the result of evaluating the last element. To get around this issue, we can use a *sapply* command:

[view raw][view raw 1] [oneForMany\_second.R][oneForMany_second.R] hosted with ❤ by  [GitHub][]

This time, all expressions have been evaluated and all calculations have been performed. However, due to the column names mismatches (the function did not assigned the proper names to the columns), the data is not in desired format and one would need to parse it in order to use it (which can be done, but it is not what we want). 
`ex1_res02 <- oneForMany_second(ex1_df, ex1_vars, ex1_formulas)`

Similarly, if we use a list instead of a vector, all expressions are evaluated and all values are be calculated.

[view raw][view raw 2] [oneForMany\_third.R][oneForMany_third.R] hosted with ❤ by  [GitHub][]

But the resulting data is also not in desired format and should be parsed. 
`ex1_res03 <- oneForMany_third(ex1_df, ex1_vars, ex1_formulas)`

Now, let's try something slightly different. Let's use *eval* to create a function that returns a data frame with desired columns (user specified functions):

[view raw][view raw 3] [oneForMany\_fourth.R][oneForMany_fourth.R] hosted with ❤ by  [GitHub][]

`ex1_res04 <- oneForMany_fourth(ex1_df, ex1_vars, ex1_formulas)`

This seems to work. We can test to confirm that the obtain results are the same as the hard coded ones:  
`testthat::expect_that(ex1_trueResult, testthat::equals(ex1_res04))`

Alternatively, instead of trying to replicate the form of the hard-coded function, we could use the *do.call* function. For do.call, we need to specify the function we use and provide a list of arguments.

[view raw][view raw 4] [oneForMany\_fifth.R][oneForMany_fifth.R] hosted with ❤ by  [GitHub][]

Note that we used the *mutate* instead of *summarise* function even thought we want data summary and this choice required additional filtering step. While *summarise* function perform group-wise summaries and creates a new data frame based on specified variables and functions, it fails when functions return data of different length (in our case, *sum(hp)* and *mean(mpg)* return single values, while *mpg+2\*carb* result in multiple values).The *mutate*function adds columns to the existing data frame, so unequal lengths are not the issue.Another important thing to mention is that the arguments list has to be a named list.

We can run it and test the results 
`ex1_res05 <- oneForMany_fifth(ex1_df, ex1_vars, ex1_formulas) testthat::expect_that(ex1_trueResult, testthat::equals(ex1_res05))`

This and another example are also available at:

<table style="border-collapse:collapse; border-spacing:0px; margin:0px; padding:0px; font-size:12px; line-height:1.4; color:rgb(51,51,51); border:0px"> 
 <tbody> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; padding-top:4px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding-right:10px!important; padding-bottom:1px!important; padding-left:10px!important; background:transparent"> # Using eval to create a function that will be able to produce a data summary (using the plyr::ddply function)</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> # based on the user specified data frame, variables (that correspond to table column names) over which they </td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> # want to split and summarize the data (character vector), and a two-column table (data frame) which contains </td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> # the user specified name of the function used for summary (first column) and the function used for summary (second column)</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> #</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> # For more info, see http://www.vesnam.com/Rblog/one-function-to-run-them-all/</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> #</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> &nbsp;</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> # Data, variable, and formulas used for data summary</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> ex1_df &lt;- datasets::mtcars</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> ex1_vars &lt;- "gear"</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> ex1_formulas &lt;- data.frame(labs=c("cust_val", "total_hp", "avg_mpg"), ress=c("mpg+2*carb", "sum(hp)", "mean(mpg)"), stringsAsFactors = FALSE)</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> &nbsp;</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> # This is the expected result</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> ex1_trueResult &lt;- plyr::ddply(ex1_df, .variable=c("gear"), function(x)</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> data.frame(cust_val = x$mpg+2*x$carb,</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> total_hp = sum(x$hp),</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> avg_mpg = mean(x$mpg)))</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> &nbsp;</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> # Testing version #!</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> # It passes an expression vector to eval and only the last element is evaluated.</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> ex1_res01 &lt;- oneForMany_first(ex1_df, ex1_vars, ex1_formulas)</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> &nbsp;</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> # Testing version #2</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> # Using sapply to overcome the eval expression vector issue</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> # Everything is calculated, but data is not in desired format (due to column names mismatches)</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> # and one would need to parse it</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> ex1_res02 &lt;- oneForMany_secpmd(ex1_df, ex1_vars, ex1_formulas)</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> &nbsp;</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> # Testing version #3</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> # Using a list instead of vector; everything is calculated, but data is not in desired format</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> # and one would need to parse it</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> ex1_res03 &lt;- oneForMany_third(ex1_df, ex1_vars, ex1_formulas)</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> &nbsp;</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> # Testing version #4</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> # Using eval to create a function that returns a data frame with desired columns (user specified functions)</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> ex1_res04 &lt;- oneForMany_fourth(ex1_df, ex1_vars, ex1_formulas)</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> # We can test to confirm that the results are OK</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> testthat::expect_that(ex1_trueResult, testthat::equals(ex1_res04))</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> &nbsp;</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> # Testing version #5</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> # Using do.call to call plyer</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> ex1_res05 &lt;- oneForMany_fifth(ex1_df, ex1_vars, ex1_formulas)</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> # We can test to confirm that the results are OK</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> testthat::expect_that(ex1_trueResult, testthat::equals(ex1_res05))</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> &nbsp;</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> # Different data set</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> ex2_df &lt;- datasets::quakes</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> ex2_vars &lt;- c("lat", "long")</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> ex2_formulas &lt;- data.frame(what=c("l1", "l2", "l3", "l4"),</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> how=c("lat+2*long", "depth^2/log(mag+1)", "(lat+long)/(stations*0.5)", "mean(mag)"), stringsAsFactors = FALSE)</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> &nbsp;</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> ex2_trueResult &lt;- plyr::ddply(ex2_df, .variable=c("lat", "long"), function(x)</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> data.frame(l1 = x$lat+2*x$long,</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> l2 = x$depth^2/log(x$mag+1),</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> l3 = (x$lat+x$long)/(x$stations*0.5),</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> l4 = mean(x$mag)))</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> &nbsp;</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> ex2_res04 &lt;- oneForMany_fourth(ex2_df, ex2_vars, ex2_formulas)</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> testthat::expect_that(ex2_trueResult, testthat::equals(ex2_res04))</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> &nbsp;</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> ex2_res05 &lt;- oneForMany_fifth(ex2_df, ex2_vars, ex2_formulas)</td> 
 </tr> 
 <tr> 
 <td style="">&nbsp;</td> 
 <td style="border:0px; position:relative; vertical-align:top; overflow:visible; font-family:Consolas,'Liberation Mono',Menlo,Courier,monospace; word-wrap:normal; white-space:pre; padding:1px 10px!important; background:transparent"> testthat::expect_that(ex2_trueResult, testthat::equals(ex2_res05))</td> 
 </tr> 
 </tbody> 
</table>

[view raw][view raw 5] [oneForMany.R][] hosted with ❤ by  [GitHub][]

[sysBio]: http://www.vesnam.com/Rblog/modeling-in-systems-biology/
[makeODEs.R]: https://github.com/Vessy/sysBio/blob/master/R/makeODEs.R
[simulateModel.R]: https://github.com/Vessy/sysBio/blob/master/R/simulateModel.R
[datasets]: https://stat.ethz.ch/R-manual/R-devel/library/datasets/html/00Index.html
[plyr]: https://cran.r-project.org/web/packages/plyr/index.html
[view raw]: https://gist.github.com/Vessy/225797327570867ecd23894de126bee9/raw/b9b4d8db29da2387cd5acbc77cbd737288aaff0a/oneForMany_first.R
[oneForMany_first.R]: https://gist.github.com/Vessy/225797327570867ecd23894de126bee9#file-oneformany_first-r
[GitHub]: https://github.com/
[view raw 1]: https://gist.github.com/Vessy/225797327570867ecd23894de126bee9/raw/b9b4d8db29da2387cd5acbc77cbd737288aaff0a/oneForMany_second.R
[oneForMany_second.R]: https://gist.github.com/Vessy/225797327570867ecd23894de126bee9#file-oneformany_second-r
[view raw 2]: https://gist.github.com/Vessy/225797327570867ecd23894de126bee9/raw/b9b4d8db29da2387cd5acbc77cbd737288aaff0a/oneForMany_third.R
[oneForMany_third.R]: https://gist.github.com/Vessy/225797327570867ecd23894de126bee9#file-oneformany_third-r
[view raw 3]: https://gist.github.com/Vessy/225797327570867ecd23894de126bee9/raw/b9b4d8db29da2387cd5acbc77cbd737288aaff0a/oneForMany_fourth.R
[oneForMany_fourth.R]: https://gist.github.com/Vessy/225797327570867ecd23894de126bee9#file-oneformany_fourth-r
[view raw 4]: https://gist.github.com/Vessy/225797327570867ecd23894de126bee9/raw/b9b4d8db29da2387cd5acbc77cbd737288aaff0a/oneForMany_fifth.R
[oneForMany_fifth.R]: https://gist.github.com/Vessy/225797327570867ecd23894de126bee9#file-oneformany_fifth-r
[view raw 5]: https://gist.github.com/Vessy/225797327570867ecd23894de126bee9/raw/b9b4d8db29da2387cd5acbc77cbd737288aaff0a/oneForMany.R
[oneForMany.R]: https://gist.github.com/Vessy/225797327570867ecd23894de126bee9#file-oneformany-r