r - plot lmer result in faceted ggplot -


i'm analysing repeated measures drug trials data , i'm not sure how plot lmer results when using faceted ggplots. have made initial plot of individual slopes master dataset, i'm doing lmer analyses separately sex.

using publicly available data, has 2 treatment groups compared 4 have, replicable example below. uses reshape2, lme4, , ggplot2 packages.

catanx <- read.fwf(file=("http://www.stat.ufl.edu/~winner/data/cats_anxiety1.dat"),                widths=c(-6,2,-5,3,-5,3,-7,1,-7,1,-7,1,-7,1,-7,1,-6,2,-6,2,-6,2,-6,2,-6,2)) colnames(catanx) <- c('id','weight','age_months','gender','environment','origin','treatment','result','emotime1','emotime2',                   'emotime3','emotime4','emotime5') library("reshape2") catanxrm <- melt(catanx, id.vars=c("id", "gender", "treatment"), measure.vars=c("emotime1", "emotime2", "emotime3",                                                                             "emotime4", "emotime5")) catanxrm$sex <- with(catanxrm, ifelse(gender==1, "neut female", ifelse(gender==2, "neut male", "whole female"))) catanxrm$time <- with(catanxrm, ifelse(variable=="emotime1", 1, ifelse(variable=="emotime2", 2, ifelse(variable=="emotime3", 3,                                   ifelse(variable=="emotime4", 4,5))))) catanxrm.male <- subset(catanxrm, gender=="2") library("lme4") male.lmer <- lmer(value ~ treatment * time + (time + 1|id), data=catanxrm.male) library("ggplot2") anxscores<-ggplot(catanxrm, aes(time, value, colour=sex))+ geom_line(aes(group = id))+ labs(x="time anxiety measured", y="anxiety score", title="effect of zylkene on anxiety")+  facet_grid(. ~ treatment) anxscores

information dataset is here.

how plot correct summary line lmer in both facets, differ on basis of treatment?

in real life example, i'll analysing females there 2 sets of lines plot per facet.

create data frame (e.g. lines.df) intercept (e.g. int) , slope (slo) variables each line of df corresponds 1 facet, plot on top:

+ geom_abline(aes(intercept = int, slope = slo), data = lines.df)

Comments

Post a Comment

Popular posts from this blog

jquery - How do you format the date used in the popover widget title of FullCalendar? -

Image
in fullcalendar, how format date used in popover widget title when clicking '+n more'? code: <div class="fc-popover fc-more-popover"> <div class="fc-header fc-widget-header"> <span class="fc-close fc-icon fc-icon-x"></span> <span class="fc-title">tuesday, april 7</span> ... </div> ... </div> you use daypopoverformat that. determines date format of title of popover created eventlimitclick option. must date formatting string. default value "dddd, mmmm d" english , "ll" other languages.
Read more

Bubble Sort Manually a Linked List in Java -

Image
this first question here. trying manually sort linked list of integers in java , can not figure out wrong code. suggestions? don't error, still have output unordered. tried few different ways, nothing worked. appreciate if can me that. public class node { int data; node nextnode; public node(int data) { this.data = data; this.nextnode = null; } public int getdata() { return this.data; } } // node class public class datalinkedlist implements datainterface { private node head; private int size; public datalinkedlist(){ this.head = null; this.size = 0; } public void add(int data) { node node = new node(data); if (head == null) { head = node; } else { node currentnode = head; while(currentnode.nextnode != null) { currentnode = currentnode.nextnode; } currentnode.nextnode = node; ...
Read more

asp.net mvc - SSO between MVCForum and Umbraco7 -

i need integrate sso between mvcforum (forum of company) , umbraco7 (manages cms company). however struggling shared authentication work. tried several ways suggested on other threads still not running expected. here did: 1. in mvcforum (forum.mywebsite.com) 1.1. machine key in web.xml: <machinekey validationkey="the same key umbraco7" decryptionkey="the same key umbraco7" validation="sha1" decryption="aes" /> 1.2. authenticate forms in web.xml <forms name="my-sso-auth" protection="all" path="/" timeout="43200" domain="forum.mywebsite.com" loginurl="/members/logon" enablecrossappredirects="true" /> 2. in umbraco7 (news.mywebsite.com) 2.1. machine key in web.xml: <machinekey validationkey="the same key mvcforum" decryptionkey="the same key mvcforum" validation="sha1" decryption="aes" /> 2.2. a...
Read more