regression - Julia contrast coding -


having trouble assigning custom contrasts categorical variables in regression. looks 1 can assign effectscoding or dummycoding using contrast parameter

    model = fit!(lmm(@formula(response ~ 1 + factor + (1|sub)), data,               contrasts = dict(:factor => effectscoding()) ))

but how assign custom contrasts?

for instance, in r can do

contrasts(data$factor) <- cbind("a_vs_b"= c(0.5, -0.5, 0), "ab_vs_c"= c(-0.25, -0.25, 0.5))  > contrasts(data$factor)        a_vs_b   ab_vs_c        0.5     -0.25 b       -0.5     -0.25 c        0.0      0.50

in addition effectscoding(), dataframes supports other contrast specifications. , in general, can supply contrast matrix (of right size, k-by-(k-1) k categories) in r. example:

using rdatasets        # install pkg.add("rdatasets") using mixedmodels      # install pkg.add("mixedmodels")  iris = dataset("datasets", "iris")  contrast_matrix = [0.5 -0.25; -0.5 -0.25; 0.0 0.5] fit!(lmm(@formula(sepallength ~ 1 + species + (1|petalwidth)), iris,    contrasts = dict(:species => contrastscoding(contrast_matrix)) ))

the specific contrast matrix in question helmert coding matrix ordering of categories. in case, simpler do:

fit!(lmm(@formula(sepallength ~ 1 + species + (1|petalwidth)), iris,    contrasts = dict(:species => helmertcoding()) ))

helmertcoding takes optional arguments base , levels (not named, see docs) switch around category orderings.

hope helps, actual code gave few problems on v0.7 julia, on v0.5 should work.


Comments

Post a Comment

Popular posts from this blog

PHP and MySQL WP -

i trying accomplish 2 things here code. writing php directly wp page using insert php plugin. 1) query mysql database , return results 1 column 4 columns. 2) make each result own hyperlink directs user new wp page runs new query on page. i able query show results , turn results dynamic hyperlink, cannot formatting down. right returns result 1 column. here have: [insert_php] $servername = "localhost"; $username = "login"; //edited privacy $password = "password"; //edited privacy $database = "database"; //edited privacy // create connection $conn = new mysqli($servername, $username, $password, $database); // check connection if ($conn->connect_error) { die("connection failed: " . $conn->connect_error); } echo "connected successfully"; //get results database $result = mysqli_query($conn,"select mine mines order mine"); $all_property = array(); //declare array saving property //showing proper...
Read more

android - InAppBilling registering BroadcastReceiver in AndroidManifest -

i'm implementing inapp billing in android app. works well, however, i'm trying decouple broadcast receiver activity manifest. suggestion in android's trivialdrive sample: // important: dynamically register broadcast messages updated purchases. // register receiver here instead of <receiver> in manifest // because call getpurchases() @ startup, therefore can ignore // broadcasts sent while app isn't running. // note: registering listener in activity bad idea, done here // because sample. regardless, receiver must registered after // iabhelper setup, before first call getpurchases(). currently there's class extends broadcastreceiver : public class iabbroadcastreceiver extends broadcastreceiver { /** * intent action receiver should filter for. */ public static final string action = "com.android.vending.billing.purchases_updated"; private final iabbroadcastlistener mlistener; public iabbroadcastreceiver(iabbroadcastlistener listener) { m...
Read more

nginx - phpPgAdmin - log in works but I have to login again after clicking on any links -

i have setup postgres 9.6 , phppgadmin on centos 7 server nginx, i can log in phppgadmin, postgresql server on right of page remains crossed out , link try browse (databases, roles etc) requires log in again. i'm thinking seems sessions problem, unsure how resolve it. any pointers on how fix or should looking @ in order track down problem? solved. turns out owner of /var/lib/php/session apache - needed nginx. sorted changing ownership nginx below: chown -r nginx:nginx /var/lib/php/session/
Read more