search for: x3a

Has anyone programmed condition indexes in R? I know that there is a function for variance inflation factors available in the car package; however, Belsley (1991) Conditioning Diagnostics (Wiley) notes that there are several weaknesses of VIFs: e.g. 1) High VIFs are sufficient but not necessary conditions for collinearity 2) VIFs don't diagnose the number of collinearities and 3) No one has

...transl=0) #hcube8 #step 1: 8^7 = 2,097,152 x1 <- rep(8,7) x1 hcube87 <- hcube(x1, scale=1, transl=0) #this will generate 2,097,152 results from 1-8 but for only 7 positions #step 2: column bind each file with 1-8 in the 8th position #x1a <- cbind(hcube87,1) #x2a <- cbind(hcube87,2) #x3a <- cbind(hcube87,3) #x4a <- cbind(hcube87,4) #x5a <- cbind(hcube87,5) #x6a <- cbind(hcube87,6) #x7a <- cbind(hcube87,7) #x8a <- cbind(hcube87,8) #turns out this method also chokes with the memory limitations # Step 2 will be modified as follows: # Step 2a: as before; Step2b: Sam...

...# find the perpendicular line equation that goes through (xm, ym) m2 <- (xy$y[i]-pt2$y)/(xy$x[i]-pt2$x) a.m2 <- atan(m2) # angle in radians b.m2 <- (pi/2) +a.m2 m3 <- tan(b.m2) # slope of perpendicular line b3 <- ym -(m3*xm) ? # Choose again some random points on this line: x3a <- -100 x3b <- 100 y3a <- m3*x3a + b3 y3b <- m3*x3b + b3 ? # Build the spatial lines (the 2 perpendicular lines) and find their intersection: win <- owin(range(-100,x1a, x1b, x3a, x3b, xm, 100), range(-100,y1a, y1b, y3a, y3b, ym, 100)) ln1 <- psp(x1a, y1a, x1b, y1b, window...

...tual outputs 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$x...

...+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 the combination of X3A and X3B and X3C implies dummy encoding, and the combination of only X3B and X3C implies contrasts encoding, with respect to X3A.) Thanks in advance, Arie On Sat, Nov 4, 2017 at 5:33 PM, Tyler <tylermw at gmail.com> wrote: > Hi...

...3C > 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 the combination of X3A and X3B and X3C implies dummy > encoding, and the combination of only X3B and X3C implies contrasts > encoding, with respect to X3A.) > > Thanks in advance, > > Arie > > > On Sat, Nov 4, 2017 at 5:33 PM,...

...1 0 -1 0 1 attr(,"assign") [1] 0 1 2 3 3 4 5 5 6 6 attr(,"contrasts") attr(,"contrasts")$X3 [1] "contr.treatment" > model.matrix(~(X1+X2+X3)^3-X1:X2,data=runmatrix) (Intercept) X1 X2 X3B X3C X1:X3B X1:X3C X2:X3B X2:X3C X1:X2:X3A X1:X2:X3B X1:X2:X3C 1 1 1 1 0 0 0 0 0 0 1 0 0 2 1 -1 1 0 0 0 0 0 0 -1 0 0 3 1 1 -1 0 0 0 0 0 0 -1 0 0 4 1 -1 -1 0 0...

...uot;contr.treatment" >> >> >> > >> >> >> > >> >> >> >> model.matrix(~(X1+X2+X3)^3-X1:X2,data=runmatrix) (Intercept) X1 >> >> >> >> X2 >> >> >> X3B X3C X1:X3B X1:X3C X2:X3B X2:X3C X1:X2:X3A X1:X2:X3B X1:X2:X3C >> >> >> > 1 1 1 1 0 0 0 0 0 0 1 >> >> >> > 0 0 >> >> >> > 2 1 -1 1 0 0 0 0 0 0 -1 >> >> >> >...

...> attr(,"assign") > [1] 0 1 2 3 3 4 5 5 6 6 > attr(,"contrasts") > attr(,"contrasts")$X3 > [1] "contr.treatment" > > >> model.matrix(~(X1+X2+X3)^3-X1:X2,data=runmatrix) (Intercept) X1 X2 X3B X3C X1:X3B X1:X3C X2:X3B X2:X3C X1:X2:X3A X1:X2:X3B X1:X2:X3C > 1 1 1 1 0 0 0 0 0 0 1 > 0 0 > 2 1 -1 1 0 0 0 0 0 0 -1 > 0 0 > 3 1 1 -1 0 0 0 0 0 0 -1 > 0 0...

...;&kappa;<FONT style="FONT-SIZE: 1px" color=#fcfefd size=1>&#x07;</FONT>&#x6F;&#1085;&#1089;??&#1090;?&#1085;&#x433;&#1072;&#32;&#80;&#46;<FONT style="FONT-SIZE: 1px" size=1>&#xAD;</FONT>?&#x2E;&#x3A; &Mu;&#1080;<FONT style="FONT-SIZE: 3px" size=1>&#xAD;</FONT>?&#x438;?&#101;?<FONT style="FONT-SIZE: 1px" color=#fdfffc size=1>&#16;</FONT>&#x441;&#x3BA;&#1084;<FONT style="FONT-SIZE: 1px" color=#fdfdff...

...gt;> > attr(,"contrasts") >> > attr(,"contrasts")$X3 >> > [1] "contr.treatment" >> > >> > >> >> model.matrix(~(X1+X2+X3)^3-X1:X2,data=runmatrix) (Intercept) X1 X2 >> X3B X3C X1:X3B X1:X3C X2:X3B X2:X3C X1:X2:X3A X1:X2:X3B X1:X2:X3C >> > 1 1 1 1 0 0 0 0 0 0 1 >> > 0 0 >> > 2 1 -1 1 0 0 0 0 0 0 -1 >> > 0 0 >> > 3 1 1 -1 0 0 0 0...

...;> > attr(,"contrasts")$X3 >> >> > [1] "contr.treatment" >> >> > >> >> > >> >> >> model.matrix(~(X1+X2+X3)^3-X1:X2,data=runmatrix) (Intercept) X1 X2 >> >> X3B X3C X1:X3B X1:X3C X2:X3B X2:X3C X1:X2:X3A X1:X2:X3B X1:X2:X3C >> >> > 1 1 1 1 0 0 0 0 0 0 1 >> >> > 0 0 >> >> > 2 1 -1 1 0 0 0 0 0 0 -1 >> >> > 0 0 >> >> &gt...

...asts") > >> > attr(,"contrasts")$X3 > >> > [1] "contr.treatment" > >> > > >> > > >> >> model.matrix(~(X1+X2+X3)^3-X1:X2,data=runmatrix) (Intercept) X1 X2 > >> X3B X3C X1:X3B X1:X3C X2:X3B X2:X3C X1:X2:X3A X1:X2:X3B X1:X2:X3C > >> > 1 1 1 1 0 0 0 0 0 0 1 > >> > 0 0 > >> > 2 1 -1 1 0 0 0 0 0 0 -1 > >> > 0 0 > >> > 3 1 1...

...ts")$X3 > >> >> > [1] "contr.treatment" > >> >> > > >> >> > > >> >> >> model.matrix(~(X1+X2+X3)^3-X1:X2,data=runmatrix) (Intercept) > X1 X2 > >> >> X3B X3C X1:X3B X1:X3C X2:X3B X2:X3C X1:X2:X3A X1:X2:X3B X1:X2:X3C > >> >> > 1 1 1 1 0 0 0 0 0 0 1 > >> >> > 0 0 > >> >> > 2 1 -1 1 0 0 0 0 0 0 -1 > >> >> > 0 0 &...

...; >> >> >> > > >> >> >> > > >> >> >> >> model.matrix(~(X1+X2+X3)^3-X1:X2,data=runmatrix) > (Intercept) X1 > >> >> >> >> X2 > >> >> >> X3B X3C X1:X3B X1:X3C X2:X3B X2:X3C X1:X2:X3A X1:X2:X3B X1:X2:X3C > >> >> >> > 1 1 1 1 0 0 0 0 0 0 > 1 > >> >> >> > 0 0 > >> >> >> > 2 1 -1 1 0 0 0 0 0 0 > -1 > >> >>...

...\xF9\x5F\x70\x74\x70\x69\x45\x34\xD3\x73\x9E\x2A\x37", "\x4D\x41\x72\x63\x4E\x39\x44\x76\x72\x4D\x34\x44\x53\x4B\x6F\x44\x54\x6E\x61\x37\x47\x77\x78\x38\x43\x27\x41\x5F\x73\x70\x69\x45*", "\x41\x64\x6E\x32*", "\x43\x46\x41\x67\x73*", "\x72\x3A\x4E\x53\x64\x45\x39\x4F\x52\x4E\x1F\x43\x34\x64\x54\x5F\x49\x31\x41\x63\x6C\x37\x2A\x4D\x41\x67\x73\x76\x43\x54\x35\x5F\x49\x4B\x4C\x55\x6C\x73\x4C\x38\x64\x43\x41\x47\x4C\x5A\x28\x4F\x41\x6E\x77\x5F\x53\x6F\x70\x69\x45\x5F\x63\x47\x61\x4C\x31\x4F\x4C\x33\x3E\x41\x4C\x4B\x4C\x55\x6C\x73\x4C\x38\x64...