r - Post-hoc pairwise fisher exact test -

-

i used have 2 factor 2 level experiment got made 3 factor 2 level experiment.

by using paste make 4 unique groups 2 factors , run fisher test outcome being whether organism lived or died.

fisher.test(mortal$alv.dead,paste(mortal$strain,mortal$capsule))

but when wanted investigate pairwise comparisions between individual groups had make inelegant filtering 2 groups entered analysis @ time. have more groups tedious hand code each paring. here fisher test test groups in 1 analysis

fisher.test(mortal$alv.dead,paste(mortal$strain,mortal$capsule,mortal$cassette))

how set method creates , tests possible pairings?

fairly easy using function combn(). thing should take account, fact combn not return names of groups correctly when put fisher.test() call inside function.

thus need adjust element in list accordingly :

some toy data:

mortal <- data.frame(       alv.dead = sample(c("alv","dead"),30,replace=t),       train = sample(letters[1:3],30,replace=t),       capsule = sample(letters[4:5],30,replace=t),       cassette = sample(letters[6:7],30,replace=t)       )

some variables

mortal$groups <- paste(mortal$train,mortal$capsule,mortal$cassette,sep="") unique.groups <- unique(mortal$groups)

and trick :

combn(unique.groups,2,function(x){     id <- mortal$groups %in% x     test <- fisher.test(table(mortal$alv.dead[id],mortal$groups[id]))     test$data.name <-       paste(         unique(           as.character(mortal$groups[id])         ),collapse="-")     return(test)}   ,simplify=false)

Comments

Post a Comment

Popular posts from this blog

ruby - When to use an ORM (Sequel, Datamapper, AR, etc.) vs. pure SQL for querying -

-
a colleague of mine designing sql queries 1 below produce reports, displayed in excel files through external data query. @ present, reporting processes on db required (no crud operations). i trying convince him better use ruby orm in order able display data in rails/sinatra app. despite obvious advantages in displaying data, advantages there him in learning use orm sequel or datamapper? the sql queries writing quite complex, , being relatively new sql, complains time-consuming , confusing. possible write extremely complex queries orm? , if so, suitable(i have heard sequel legacy dbs)? , advantages of learning ruby , using orm versus sticking plain sql, in making complex database queries? i'm datamapper maintainer, , think complex reporting should use sql. while think someday we'll have dsl provides power , conciseness of sql, i've seen far requires write more ruby code sql complex queries. rather maintain 5 line sql query 10-15 lines of ruby code desc...
Read more

php - PHPDoc: @return void necessary? -

-
is necessary this: /** * ... * * @return void */ i have quite few methods don't have return value, , seems redundant put in comment. considered bad form leave out? if makes clear documentation, leave in, isn't strictly necessary. it's entirely subjective decision. personally, leave out. edit stand corrected. after little googling, wikipedia page says: @return [type description] tag should not used constructors or methods defined void return type. the phpdoc.org website says: @return datatype description @return datatype1|datatype2 description the @return tag used document return value of functions or methods. @returns alias @return support tag formats of other automatic documentors the datatype should valid php type (int, string, bool, etc), class name type of object returned, or "mixed". if want explicitly show multiple possible return types, list them pipe-delimited without spaces (e.g. "@retur...
Read more

c++ - Convert big endian to little endian when reading from a binary file -

-
this question has answer here: how convert between big-endian , little-endian values in c++? 28 answers i've been looking around how convert big-endian little-endians. didn't find solve problem. seem there's many way can conversion. anyway following code works ok in big-endian system. how should write conversion function work on little-endian system well? this homework, since systems @ school running big-endian system. it's got curious , wanted make work on home computer also #include <iostream> #include <fstream> using namespace std; int main() { ifstream file; file.open("file.bin", ios::in | ios::binary); if(!file) cerr << "not able read" << endl; else { cout << "opened" << endl; int i_var; double d_var; while(!file.eof()) {...
Read more