Retrieve a value from S4 class output in R -

-

i writing simulation in r. have decided use s4 class return 2 values in function. when run simulation, wondering how can retrieve values output calculate moments of distribution them such mean?

setclass(class="coalescent",           representation(            tmrca="numeric",            t_total="numeric"            )          )

the output looks below:

> tmrca_sim <-replicate(10000, coalescent_fun(n, ne)) > head(tmrca_sim) [[1]] object of class "coalescent" slot "tmrca": [1] 6.723592  slot "t_total": [1] 9.693661   [[2]] object of class "coalescent" slot "tmrca": [1] 1.592346  slot "t_total": [1] 11.50406

what want extract values of "tmrca" , "t_total" , calculate mean. of course use many other ways simulation want learn use of classes @ same time.

you can extract data matrix:

mx <- sapply(tmrca_sim, function(x) sapply(slotnames(x), slot, object=x))

with made data mx looks like:

             [,1]      [,2]      [,3]      [,4]      [,5]  tmrca   0.3823880 0.3403490 0.5995658 0.1862176 0.6684667  t_total 0.8696908 0.4820801 0.4935413 0.8273733 0.7942399

and can use rowmeans or apply:

rowmeans(mx) apply(mx, 1, mean)  # equivalent, though slower , more flexible apply(mx, 1, var)

though note in comments, slow way of doing things in r. want coalescent_fun produce objects 2 vectors many entries each, not 1 object per simulation.


Comments

Post a Comment

Popular posts from this blog

Java 3D LWJGL collision -

-
i making 3d java game lwjgl library, , wondering how add collision detection, player not go through models. i using obj models. here objloader class, loads models: package renderengine; import java.io.bufferedreader; import java.io.file; import java.io.filenotfoundexception; import java.io.filereader; import java.util.arraylist; import java.util.list; import models.rawmodel; import org.lwjgl.util.vector.vector2f; import org.lwjgl.util.vector.vector3f; public class objloader { public static rawmodel loadobjmodel(string filename, loader loader){ filereader fr = null; try { fr = new filereader(new file("res/"+filename+".obj")); } catch (filenotfoundexception e) { system.err.println("couldn't load file!"); e.printstacktrace(); } bufferedreader reader = new bufferedreader(fr); string line; list<vector3f> vertices = new arraylist<vector3f...
Read more

spring - SubProtocolWebSocketHandler - No handlers -

-
i have ugly error during deploying spring app on jboss. 18:11:16,025 error [org.apache.catalina.core.containerbase.[jboss.web].[default-host].[/consumer]] (msc service thread 1-7) exception sending context initialized event listener instance of class org.springframework.web.context.contextloaderlistener: org.springframework.context.applicationcontextexception: failed start bean 'subprotocolwebsockethandler'; nested exception java.lang.illegalargumentexception: no handlers @ org.springframework.context.support.defaultlifecycleprocessor.dostart(defaultlifecycleprocessor.java:176) [spring-context-4.1.3.release.jar:4.1.3.release] @ org.springframework.context.support.defaultlifecycleprocessor.access$200(defaultlifecycleprocessor.java:51) [spring-context-4.1.3.release.jar:4.1.3.release] @ org.springframework.context.support.defaultlifecycleprocessor$lifecyclegroup.start(defaultlifecycleprocessor.java:346) [spring-context-4.1.3.release.jar:4.1.3.release] @ org.s...
Read more

methods - python can't use function in submodule -

-
my folder structure this: pythonstuff/ program.py modulea/ __init__.py foo.py bar.py here's code bar.py : def hello(): print("hello, world!") here's code program.py : #!/usr/bin/env python3 modulea import bar bar.hello() i run $ python3 program.py somehow error: file "program.py", line 3, in <module> bar.hello() attributeerror: 'module' object has no attribute 'hello' edit: __init__.py file empty. edit2: after trying realized had bar.py in root directory contained hello() method. bar.py in modulea/ directory empty. add it import os __all__ = [] module in os.listdir(os.path.dirname(__file__)): if module != '__init__.py' , module[-3:] == '.py': __all__.append(module[:-3]) the init .py files required make python treat directories containing packages; done prevent directories common name, such string, unintentionally hiding v...
Read more