search for: x3c

...he factor encoding changes for a higher level interaction. Consider the following minimal reproducible example: -------------- > runmatrix = expand.grid(X1=c(1,-1),X2=c(1,-1),X3=c("A","B","C"))> model.matrix(~(X1+X2+X3)^3,data=runmatrix) (Intercept) X1 X2 X3B X3C X1:X2 X1:X3B X1:X3C X2:X3B X2:X3C X1:X2:X3B X1:X2:X3C 1 1 1 1 0 0 1 0 0 0 0 0 0 2 1 -1 1 0 0 -1 0 0 0 0 0 0 3 1 1 -1 0 0 -1 0 0 0 0 0 0 4 1 -1 -1...

...> for a higher level interaction. Consider the following minimal reproducible > example: > > -------------- > >> runmatrix = expand.grid(X1=c(1,-1),X2=c(1,-1),X3=c("A","B","C"))> model.matrix(~(X1+X2+X3)^3,data=runmatrix) (Intercept) X1 X2 X3B X3C X1:X2 X1:X3B X1:X3C X2:X3B X2:X3C X1:X2:X3B X1:X2:X3C > 1 1 1 1 0 0 1 0 0 0 0 > 0 0 > 2 1 -1 1 0 0 -1 0 0 0 0 > 0 0 > 3 1 1 -1 0 0 -1 0 0 0 0 > 0...

...Below are the three results, similar to your first mail. (The first two are basically the same, of course.) Please pick one result which you think is not consistent with the heuristic and please give what you think is the correct result: model.matrix(~(X1+X2+X3)^3-X1:X3) (Intercept) X1 X2 X3B X3C X1:X2 X2:X3B X2:X3C X1:X2:X3B X1:X2:X3C model.matrix(~(X1+X2+X3)^3-X2:X3) (Intercept) X1 X2 X3B X3C X1:X2 X1:X3B X1:X3C X1:X2:X3B X1:X2:X3C model.matrix(~(X1+X2+X3)^3-X1:X2) (Intercept) X1 X2 X3B X3C X1:X3B X1:X3C X2:X3B X2:X3C X1:X2:X3A X1:X2:X3B X1:X2:X3C (I take it that th...

...lar to your first mail. (The first > two are basically the same, of course.) Please pick one result which > you think is not consistent with the heuristic and please give what > you think is the correct result: > > model.matrix(~(X1+X2+X3)^3-X1:X3) > (Intercept) > X1 X2 X3B X3C > X1:X2 X2:X3B X2:X3C > X1:X2:X3B X1:X2:X3C > > model.matrix(~(X1+X2+X3)^3-X2:X3) > (Intercept) > X1 X2 X3B X3C > X1:X2 X1:X3B X1:X3C > X1:X2:X3B X1:X2:X3C > > model.matrix(~(X1+X2+X3)^3-X1:X2) > (Intercept) > X1 X2 X3B X3C > X1:X3B X1:X3C X2...

...----------- >> >> >> > >> >> >> >> runmatrix = expand.grid(X1=c(1,-1),X2=c(1,-1),X3=c("A","B","C"))> >> >> >> model.matrix(~(X1+X2+X3)^3,data=runmatrix) (Intercept) X1 X2 X3B >> >> >> X3C >> >> >> X1:X2 X1:X3B X1:X3C X2:X3B X2:X3C X1:X2:X3B X1:X2:X3C >> >> >> > 1 1 1 1 0 0 1 0 0 0 0 >> >> >> > 0 0 >> >> >> > 2 1 -1 1 0 0 -1 0...

...imal >> reproducible >> > example: >> > >> > -------------- >> > >> >> runmatrix = expand.grid(X1=c(1,-1),X2=c(1,-1),X3=c("A","B","C"))> >> model.matrix(~(X1+X2+X3)^3,data=runmatrix) (Intercept) X1 X2 X3B X3C >> X1:X2 X1:X3B X1:X3C X2:X3B X2:X3C X1:X2:X3B X1:X2:X3C >> > 1 1 1 1 0 0 1 0 0 0 0 >> > 0 0 >> > 2 1 -1 1 0 0 -1 0 0 0 0 >> > 0 0 >> > 3 1...

...ple: >> >> > >> >> > -------------- >> >> > >> >> >> runmatrix = expand.grid(X1=c(1,-1),X2=c(1,-1),X3=c("A","B","C"))> >> >> model.matrix(~(X1+X2+X3)^3,data=runmatrix) (Intercept) X1 X2 X3B X3C >> >> X1:X2 X1:X3B X1:X3C X2:X3B X2:X3C X1:X2:X3B X1:X2:X3C >> >> > 1 1 1 1 0 0 1 0 0 0 0 >> >> > 0 0 >> >> > 2 1 -1 1 0 0 -1 0 0 0 0 >> >> &g...

...; >> > example: > >> > > >> > -------------- > >> > > >> >> runmatrix = expand.grid(X1=c(1,-1),X2=c(1,-1),X3=c("A","B","C"))> > >> model.matrix(~(X1+X2+X3)^3,data=runmatrix) (Intercept) X1 X2 X3B X3C > >> X1:X2 X1:X3B X1:X3C X2:X3B X2:X3C X1:X2:X3B X1:X2:X3C > >> > 1 1 1 1 0 0 1 0 0 0 0 > >> > 0 0 > >> > 2 1 -1 1 0 0 -1 0 0 0 0 > >> > 0 0 >...

...t; > >> >> > -------------- > >> >> > > >> >> >> runmatrix = expand.grid(X1=c(1,-1),X2=c(1,-1),X3=c("A","B","C"))> > >> >> model.matrix(~(X1+X2+X3)^3,data=runmatrix) (Intercept) X1 X2 X3B > X3C > >> >> X1:X2 X1:X3B X1:X3C X2:X3B X2:X3C X1:X2:X3B X1:X2:X3C > >> >> > 1 1 1 1 0 0 1 0 0 0 0 > >> >> > 0 0 > >> >> > 2 1 -1 1 0 0 -1 0 0 0 0...

...gt; >> > > >> >> >> >> runmatrix = expand.grid(X1=c(1,-1),X2=c(1, > -1),X3=c("A","B","C"))> > >> >> >> model.matrix(~(X1+X2+X3)^3,data=runmatrix) (Intercept) X1 X2 > X3B > >> >> >> X3C > >> >> >> X1:X2 X1:X3B X1:X3C X2:X3B X2:X3C X1:X2:X3B X1:X2:X3C > >> >> >> > 1 1 1 1 0 0 1 0 0 0 0 > >> >> >> > 0 0 > >> >> >> > 2 1 -1 1 0...

...puts below for your reference. Thanks in advanced. regards, Jei ********************************** #1 > library(survival) > library(boot) > filename <- type2a # this file has 100 records > > r <- 50 > > data.cox <- coxph(Surv(time,cens==1)~x1+x2+x3a+x3b+x3c, data=filename) > data.surv <- survfit(data.cox) > data.cens <- survfit(Surv(time-0.001*(cens==1),cens!=1),data=filename) > > beta.fun <- function(filename) { + cox <- coef(coxph(Surv(filename$time,filename$cens==1) ~ filename$x1+filename$x2+filename$x3a+filename$x3b+filen...

I am trying to do something in R and would appreciate a push into the right direction. I hope some of you experts can help. I have two distributions obtrained from 10000 datapoints each (about 10000 datapoints each, non-normal with multi-model shape (when eye-balling densities) but other then that I know little about its distribution). When plotting the two distributions together I can see that

From: Matt Fleming <matt.fleming at linux.intel.com> These patches are based on the elflink branch. This set of patches is my attempt at moving the command-line interface functionality out of the core and into an ELF module to reduce the size of the core. The most interesting patch is [PATCH 4/4] which moves the cli code out of core/elflink and into com32/elflink/modules. [PATCH 4/4] is

From: Matt Fleming <matt.fleming at linux.intel.com> This is a series of patches that, * shrink the core by moving things into an ldlinux ELF module * begin wiring up some of the C versions of various functions The core now only contains essential code and loads the ldlinux module to do everything else, like providing a command line interface and loading kernels. The config file parsing

...xF9\x51\x29\xE5\xE5\x46\x70\xF1" -+ (uint8_t*)"abcdbcdecdefdefgefghfghighijhijkijkljklmklmnlmnomnopnopq", 1, -+ (uint8_t*)"\x84\x98\x3E\x44\x1C\x3B\xD2\x6E\xBA\xAE\x4A\xA1\xF9\x51\x29\xE5\xE5\x46\x70\xF1" - }, { -- "a", 1000000, -- "\x34\xAA\x97\x3C\xD4\xC4\xDA\xA4\xF6\x1E\xEB\x2B\xDB\xAD\x27\x31\x65\x34\x01\x6F" -+ (uint8_t*)"a", 1000000, -+ (uint8_t*)"\x34\xAA\x97\x3C\xD4\xC4\xDA\xA4\xF6\x1E\xEB\x2B\xDB\xAD\x27\x31\x65\x34\x01\x6F" - }, { -- "01234567012345670123456701234567012345670123456701234567012...

From: Matt Fleming <matt.fleming at linux.intel.com> This series fixes some bugs and switches the elflink branch to be entirely ELF modules. It applies on top of, http://syslinux.zytor.com/archives/2011-April/016369.html The deletions in the diff stat below are mainly from deleting com32/elflink/modules (finally!). Now there should be no duplicate code because we don't need COM32 and