search for: myp

...ation_probesets <- scan('c://ttest.csv',skip=1,what='') dataDirectory <- system.file("data", package = "annotationTools") annotation_stuff <- read.csv(paste(dataDirectory, annotation_HGU133Plus2, sep = "/"), colClasses = "character") myPS <- c(annotation_probesets) getGENESYMBOL(myPS, annotation_stuff) getGENEONTOLOGY(myPS, annotation_HGU133Plus2) ###############error message######################## > annotation_stuff <- read.csv(paste(dataDirectory, + annotation_HGU133Plus2, sep = "/"), colClasses = "char...

Dear All I am attempting to use Sekhon's package Match to create a matched data set. I keep getting an error regarding unused arguments whether I put Y in or not. Have I missed something? myps id a propensity score (for treat) calculated from logistic regression mymatch<-match(Y=status,Tr=treat,X=myps, M=1,replace=F) Error in match(Y = status, Tr = treat, X = myps, M = 1, replace = F) : unused argument(s) (Y = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,...

...vac:inf 1 0.3884 0.3884 0.7058 0.44812 Residuals 4 2.2009 0.5502 so we get the 2.2009 on 4 df SS for redisuals with mean SS 0.5502. Then I would do: mNp<-mean(y[(vac=="N")&(inf=="+")]) mNm<-mean(y[(vac=="N")&(inf=="-")]) mYp<-mean(y[(vac=="Y")&(inf=="+")]) mYm<-mean(y[(vac=="Y")&(inf=="-")]) c( mYp, mYm, mNp, mNm ) ##[1] 2.4990492 0.5532018 2.5212655 -0.3058972 c( mYp-mYm, mNp-mNm ) ##[1] 1.945847 2.827163 after whi...

Hello! I am using package "VAR". I've fitted my model: mymodel<-VAR(mydata,myp,type="const") I can extract the Log Liklihood for THE WHOLE MODEL: logLik(mymodel) How could I calculate (other than manually) the corresponding Akaike Information Criterion (AIC)? I tried AIC - but it does not take mymodel: AIC(mymodel) # numeric(0) Thank you! -- Dimitri Liakhovits...

Hi So carrying on my use of analysis of variance to check for the effects of two factors. It's made simpler by the fact that both my factors have only two levels each, creating four unique groups. I have a highly significant interaction term. In the context of the experiment, this makes sense. I can visualise the data graphically, and sure enough I can see that both factors have different

I am writing a report using Sweave with specific R output incorporated into the text using the Sexpr{} command. I have run into two specific issues: 1) If the result inside the Sexpr{} command is an integer less than 10 I would like to print the ?word? for that number ? e.g., ?seven? for 7. 2) If a p-value is ?equal to zero? to some arbitrary number of decimal places -- e.g.,

Hi, I'm running 2.04b on two RH 5.2 boxes (A + B). The boxes are control separate domains and have a router in between them. I would like for the domain controlled by A to show up in browse list on B. I've tried using remote announce on A to broadcast on the other side of the router, but domain A does not show up.... It seems to me that WINS (nmbd) should be able to hold

I put a custom HTML file in the admin folder, but when I login to the admin interface then browse to it (e.g., http://server:8000/admin/mypage.html), I get "400 - UNRECOGNISED COMMAND". Is admin (or the admin folder) able to render only the standard Icecast files? I'm using 2.4.4 on Windows 10. Thanks, Justin -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.xiph.org/p...

...+L7fIZFIo6f6WN26uzQIxzb8CX91YnUSiUQikaQuEig9CALoIrvdbpQpJpFII6YhsLqhC4pJ OLYh97l1S6lZJpFIJBKJmdOBz6FzhC4S9vlJJyhzSKSRUrWsnnB+5uy4A6vL4TSmqm5Wh6vg RSYikUgkUj2CNhZaWuqFp8ui3s2ElE4ZQiLNblY3UX1DOytndbgEf3HoDUQTe5JIJFKFYq0r NrbEisMXdIs8qDNnmE6nA39hnwqFRFqlWH26QN2B1RHUGZ9T9SCRSKT6TLvYRBOuD1mQ59An smyHv/DRmypcQJDyikRaRVh9evsCK1ZnoE4Vg0QikYbWSkPbS1kxNKGHJ+vput0uwLnv++i4 Dn9hH47AccorEolYfaRYnXm8UK0gkUikITfU5Ls+NAGE8y9HSOY8vcM+0jvlFYk0O1hdOT3L DGV1AnUSiUSarraa8mForM7ndqfTATLnuz/YhyNwnPKKRJqtrD5q7b8hq5PPJIlEIk1XW01u MMNkdb6zY3Z1doTs6iTSTGd1edZEYe7EvNkX+W/lS+RwlIEYnsPaf/0lvAqtOux3Q7LAk0gk UrUiu/rQJEzzwvzVsYO...