similar to: Contrasts in mixed effects model: difference between differences

To whom it may concern: The help page for ?legend refers to a `title.cex` parameter, which suggests that the function has such a parameter. As far as I can tell, though, it doesn't; here's an example: > plot(1,1) > legend("topright",pch=1, legend="something", title="my legend", title.cex=2) Error in legend("topright", pch = 1, legend =

Hello everyone, Does anyone know what would cause the skewness() function (from e1071), as well as skew() from psych, to return a value of NaN? I have a vector of positively-skewed data (https://docs.google.com/file/d/0B6-m45Jvl3ZmYzlHRVRHRURzbVk/edit?usp=sharing) which these functions return a value for like normal: > skewness( data ) # returns 1.400405 but when I instead give those

Right, that's my point. The help page mentions a `title.cex`, like I said; saying that `cex` sets the default `title.cex` sure implies to me (and presumably to the other people whose discussion I linked) that a `title.cex` parameter exists. Since no such parameter exists, this bit in the documentation is misleading (suggesting that there is a `title.cex` parameter which can be set, when there

Hello everyone, I've conducted a Type III repeated-measures ANOVA using Anova() from the car package, based on the suggestions at http://blog.gribblelab.org/2009/03/09/repeated-measures-anova-using-r/(option 3) and http://languagescience.umd.edu/wiki/EEG#ERP_ANOVA_in_R. My ANOVA has two factors: Condition (3 levels) and Region (6 levels) and their interaction. Below is code to run the Anova

Hello everyone, I am running a mixed effects model where I have two fixed factors, one with 2 levels and one with 4, and their interaction. Let's say these are my factors and their levels: FirstFactor: 1, 2 SecondFactor: A, B, C, D For the interaction, I am interested in the four two-way comparisons, not the two four-way comparisons. In other words, I want to test whether 1A is

Hello, duplicated() does not seem to work for a long vector. For example, if you download the data from https://docs.google.com/open?id=0B6-m45Jvl3ZmNmpaSlJWMXo5bmc (a vector with about 12,000 numbers) and then run the following code which does duplicated() over the whole vector but just shows the last 30 elements: data.frame( tail(verylong, 30), tail(duplicated(verylong), 30) ) you'll see

Goodmorning everybody, I need to run an R program via SSH. Usually I open R, I run the program and I stay logged-in, waiting for the output. As a matter of fact, if I close the connection with SSH I loose the calculations and the output of my R program. What command I have to use in order to preseve the calculations and the output without staying logged-in a SSH connection? thanks in advance

Hello, I have a list of data with multiple elements, and each element in the list has multiple variables in it. Here's an example: ### Make the fake data dv <- c(1,3,4,2,2,3,2,5,6,3,4,4,3,5,6) subject <- factor(c("s1","s1","s1","s2","s2","s2","s3","s3","s3",

Great R people, I have noticed a strange behaviour in read.delim() and friends in the R 2.7.0 version. I will describe you the problem and also the solution I already found, just to be sure it is an expected behaviour and also to tell people, who may experience the same difficulty, a way to overcome it. And also to see if it is a proper behaviour or maybe a correction is needed. Here is the

I have some doubts about the validity of my procedure to estimeate linear contrasts ina a factorial design. For sake of semplicity, let's imagine a one way ANOVA with three levels. I am interested to test the significance of the difference between the first and third level (called here contrast C1) and between the first and the seconda level (called here contrast C2). I used the following

Dear all, I studied in tank prey fish behavior. Using the design described below (and R code), I want to test the effects of both habitat and predator (and interaction) on prey fish's vertical distribution, which was recorded (with repeated measures) as a categorical variable. I found that package mlogit might fit to my need but I don't know how to specify my complex design in the

>>>>> POLITZER-AHLES, Stephen [CBS] <stephen.politzerahles at polyu.edu.hk> >>>>> on Sat, 25 Mar 2017 13:25:32 +0000 writes: > Right, that's my point. The help page mentions a > `title.cex`, like I said; saying that `cex` sets the > default `title.cex` sure implies to me (and presumably to > the other people whose discussion

Hi, I'd want to configure persistent static routes, ie in config files, but I can't configure static routes, I tested: [root at centos7 ~]# cat /etc/sysconfig/network-scripts/routes-enp0s3 ADDRESS0=10.10.10.0 NETMASK0=255.255.255.0 GATEWAY0=192.168.1.1 logs says nothing: ago 27 15:38:29 centos7.vhgroup.corp NetworkManager[2299]: <info> Activation (enp0s3) Stage 5 of 5 (IPv4

> On 25 Mar 2017, at 00:39 , POLITZER-AHLES, Stephen [CBS] <stephen.politzerahles at polyu.edu.hk> wrote: > > To whom it may concern: > > > The help page for ?legend refers to a `title.cex` parameter, which suggests that the function has such a parameter. No it does not. All arguments are listed and documented, none of them is title.cex, and there's no

Hello everyone, I'm currently writing an SMB client program but I ran into a nasty bug. The problem occurs if I send an trans2 request (find first2 for example) after an andx create request; the returned error in the trans2 request is STATUS_INVALID_INFO_CLASS. I have no idea what that means...did not find any information regarding this. The error/problem does not occur when I don't

Dear developers, I found a problem that 'smbclient' loses some files when it's creating a tar ball. A patch attached this message is for samba 2.2.8a and fixes this problem. Would you please check it? I think there is a protocol mismatch with Windows 2000 at cli_list_new() in libsmb/clilist.c (FIND FIRST2/FIND NEXT2). Samba assumes that 'Last Name Offset' response field is

Dear R users, I have a reproducible data and try to create new variable "clo" is 1 if know variable is equal to "very well" or "fairly well" and getalong is 4 or 5 otherwise it is 0. rep_data<- read.table(header=TRUE, text=" id1 id2 know getalong 100000016_a1 100000016_a2 very well 4 100000035_a1 100000035_a2 fairly

Hello, The attached l32or64.c implements a com32 module that boots two different kernels with different initrds depending on whether the cpu has long mode support or not. I stumbled upon two problems while developing it with current git version (last commit 595705ffad4f63cfeb84e9bb1243df03808c2fff). The first was that syslinux_boot_linux didn't work for me. Both the command line and initrd

I have a fully working cyrus+sendmail mail server hosting all users for foo.com. Company foo has been acquired by company bar and it has been decided that bar.com is the new domain of choice. The administrator of bar.com has added a pile of redirects to send various user at bar.com addresses so they come up to my server at foo.com. My users have changed their From: addresses to user at bar.com and

Hello Amelie, I don't have an answer to your question, but I just wanted to point out this page I noticed recently ( http://hlplab.wordpress.com/2009/05/07/multinomial-random-effects-models-in-r/), which might be helpful. I'm also interested in figuring out how to do a multinomial glmm, so if you find out anything I'd be happy to hear more about it! Based on what I've found so