search for: x22

...ed as.name() but it did not work (because of size?). These are creating trouble with subsequent programs, which I tested with strings that for some reason do not have these double quotes (see very bottom). > cum1 [1] "A11*(X11*x1+X21*x2)+1*sqrt(B11*(X11*x1+X21*x2)^2+C11)A12*(X12*x1+X22*x2)+1*sqrt(B12*(X12*x1+X22*x2)^2+C12)A13*(X13*x1+X23*x2)+-1*sqrt(B13*(X13*x1+X23*x2)^2+C13)A14*(X14*x1+X24*x2)+-1*sqrt(B14*(X14*x1+X24*x2)^2+C14)A15*(X15*x1+X25*x2)+1*sqrt(B15*(X15*x1+X25*x2)^2+C15)A16*(X16*x1+X26*x2)+1*sqrt(B16*(X16*x1+X26*x2)^2+C16)A17*(X17*x1+X27*x2)+1*sqrt(B17*(X17*x1+X27*x2)^2...

Hi R users! Does anyone know command similar to cbind for adding a column to a object of the list. For example on this list: > c1 $`1` x11 x22 1 1 1 2 1 2 3 1 3 4 1 4 5 1 5 $`2` x11 x22 6 2 6 7 2 7 8 2 8 9 2 9 10 2 10 11 2 11 i would like to add column, named random created with - runif(5)-, to the first object in order to obtain something like this: > c1 $`1` x11 x22 random 1 1 1 0.5624256 2 1 2 0.6752097 3 1 3 0.100...

Dear all, I have encountered a problem when perform stepwise regression. The dataset have more 9 independent variables, but 7 observation. In R, before performing stepwise, a lm object should be given. fm <- lm(y ~ X1 + X2 + X3 + X11 + X22 + X33 + X12 + X13 + X23) However, summary(fm) will give: Residual standard error: NaN on 0 degrees of freedom Multiple R-Squared: 1, Adjusted R-squared: NaN F-statistic: NaN on 6 and 0 DF, p-value: NA In this situation, step() or stepAIC() will not give any useful information....

...0.63 -0.82 [2,] 0.18 0.49 [3,] -0.84 0.74 [4,] 1.60 0.58 [5,] 0.33 -0.31 I need to form a matrix X such that X = x11 0 0 x21 0 0 0 x11 0 0 x21 0 0 0 x11 0 0 x21 x12 0 0 x22 0 0 0 x12 0 0 x22 0 0 0 x12 0 0 x22 ... x15 0 0 x25 0 0 0 x15 0 0 x25 0 0 0 x15 0 0 x25 whe...

I have a dataset that is basically structureless. Its dimension varies from row to row and sep(s) are a mixture of tab and semi colon (;) and example is HEADER1 HEADER2 HEADER3 HEADER3 A1 B1 C1 X11;X12;X13 A2 B2 C2 X21;X22;X23;X24;X25 A3 B3 C3 A4 B4 C4 X41;X42;X43 A5 B5 C5 X51 etc., say. Note that a blank under HEADER3 corresponds to non occurance and all semi colon (;) delimited variables are under HEADER3. These values run into tens of thousands. I want to give...

...order to extract columns, but R console gives me the "+" sign at the end of the code, pointing out incomplete command, so it is not working: m<-read.table("C:/backup/symmetrical.csv") n<-subset(m, select=c("X1", "X7", "X12", X15", "X22", "X26", "X31", "X34", "X39", "X44", "x51", "X58) Therefore, I have no tried with row names yet. Any suggestions? Sorry for the inconvenience. I have read some information about this but always have the same problem with "...

...0.33749086 0.60395258 > x2 0.72017894 1.11932077 > > > > DF2 <- with(DF, data.frame(y=c(y1,y2))) > > DF2$x11 <- with(DF, c(x1, rep(0,100))) > > DF2$x21 <- with(DF, c(x2, rep(0,100))) > > DF2$x12 <- with(DF, c(rep(0,100), x1)) > > DF2$x22 <- with(DF, c(rep(0,100), x2)) > > DF2$x1 <- with(DF, c(x1, x1)) > > DF2$wh <- rep(c(0,1), each=100) > > fm1 <- lm(y~wh + x11 + x21 + x12 + x22, DF2) > > fm1 > > Call: > lm(formula = y ~ wh + x11 + x21 + x12 + x22, data = DF2) > > Coefficients: &...

Hi all, I try to represent a multiple curve graphic where the x-axis is the temperature and the different y-axes are the different X (X22,X43,X44...) some X corresponds to the same molecule (22 and 44 are for CO2 for instance) so I use the same colour for them. I wanna mix the linetype with the colour to be able to visually see the difference between X43 and X45 The best I have done up to now is this code but it crashes with :&quot...

...t function called the horizontal direct product (*~), which is some sort of variant on the Kronecker product. For example if x is 2x2 and y is 2x2 the horizontal direct product, z, of x and y is defined (in the Gauss manual) as: row 1 = x11*y11 x11*y12 x12*y11 x12*y12 row 2 = x21*y21 x21*y22 x22*y21 x22*y22 Or in R code if: x <- matrix(seq(1,4,by=1),2,2, byrow=TRUE) y <- matrix(seq(5,8,by=1),2,2, byrow=TRUE) The resulting matrix, if I had an operator, would be the following matrix z, here formed in a contrived manner: z.1 <- c(5, 6, 10, 12) z.2 <- c(21,24,28,32) z <- r...

Dear R-Users, I used the box.tidwell () function of the car Package. When I used the following formula: semi.sub.in.mi1.boxtidwell_h<-box.tidwell(RENT_LG ~ I(age+1)+I(age2+1)+X06A + I(X08B+1) + I(X22+1) + I(X24+1) + X31A, ~B_YEAR + C_X01 + C_X14 + C_X19 + C_X29A +C_X21 + C_X23 + D_X12 + D_X17 + D_X18 + D_X25 + D_X27 + D_X30 + D_X32 + D_X35, data = semi.sub.in.mi1) everything is fine. However, when I replaced the time dummy variable: semi.sub.in.mi1.boxtidwell_h<-box.tidwell(RENT_...

Hi I get a matrix with the columns named like "X13 X22 X1 X14 ...", i.e. not successively. That has it's good reasons, but now how can I rename the columns to "X1 X2 X3 ..."? Thanks a lot, David

Hello, suppose I have a multivariate multiple regression model such as the following: > DF<-data.frame(x1=rep(c(0,1),each=50),x2=rep(c(0,1),50)) > tmp<-rnorm(100) > DF$y1<-tmp+DF$x1*.5+DF$x2*.3+rnorm(100,0,.5) > DF$y2<-tmp+DF$x1*.5+DF$x2*.7+rnorm(100,0,.5) > x.mlm<-lm(cbind(y1,y2)~x1+x2,data=DF) > coef(x.mlm) y1 y2 (Intercept)

...56 ), X19 = c(288.79296875, 290.357543945312, 289.657928466797, 291.449066162109, 293.095275878906), X20 = c(291.999877929688, 292.838348388672, 293.840362548828, 294.412322998047, 294.941253662109 ), X21 = c(293.615447998047, 294.028106689453, 296.072296142578, 296.447387695312, 295.824615478516), X22 = c(294.719848632812, 295.392028808594, 297.453216552734, 297.114288330078, 296.883209228516 ), X23 = c(295.634429931641, 296.783294677734, 298.592346191406, 297.469512939453, 297.832122802734)), .Names = c("X0", "X1", "X2", "X3", "X4", "X5&quo...

Hi, Checking the llvm test-suite SingleSource/Regression/C++/EH/class_hierarchy testcase on aarch64 I noted something strange: Dump of assembler code for function _Z4funcj: 0x0000000000400650 <+0>: stp x22, x21, [sp,#-48]! 0x0000000000400654 <+4>: stp x20, x19, [sp,#16] 0x0000000000400658 <+8>: stp x29, x30, [sp,#32] 0x000000000040065c <+12>: add x29, sp, #0x20 => 0x0000000000400660 <+16>: .inst 0x2bfffff7 ; undefined Where it should be...

...: // %entry sub sp, sp, #224 // =224 stp x28, x27, [sp, #128] // 8-byte Folded Spill stp x26, x25, [sp, #144] // 8-byte Folded Spill stp x24, x23, [sp, #160] // 8-byte Folded Spill stp x22, x21, [sp, #176] // 8-byte Folded Spill stp x20, x19, [sp, #192] // 8-byte Folded Spill stp x29, x30, [sp, #208] // 8-byte Folded Spill ldrsw x8, [x0, #4424] sxtw x10, w2 <------------- w2 is the use of spilled value b...

...t; sub sp, sp, #224 // =224 > > stp x28, x27, [sp, #128] // 8-byte Folded Spill > > stp x26, x25, [sp, #144] // 8-byte Folded Spill > > stp x24, x23, [sp, #160] // 8-byte Folded Spill > > stp x22, x21, [sp, #176] // 8-byte Folded Spill > > stp x20, x19, [sp, #192] // 8-byte Folded Spill > > stp x29, x30, [sp, #208] // 8-byte Folded Spill > > ldrsw x8, [x0, #4424] > > sxtw x10, w2 <-------...

...5312, 289.657928466797, > > 291.449066162109, 293.095275878906), X20 = c(291.999877929688, > > 292.838348388672, 293.840362548828, 294.412322998047, 294.941253662109 > > ), X21 = c(293.615447998047, 294.028106689453, 296.072296142578, > > 296.447387695312, 295.824615478516), X22 = c(294.719848632812, > > 295.392028808594, 297.453216552734, 297.114288330078, 296.883209228516 > > ), X23 = c(295.634429931641, 296.783294677734, 298.592346191406, > > 297.469512939453, 297.832122802734)), .Names = c("X0", "X1", > > "X2", &quo...

...nd the cost of this distribution plan is $1150. However, I would like to add the following two constraints: # 1. weighted sums by column # w is a vector of known constants, i.e., w = c(1.2, .9, .7, 2.3) # r is also known, i.e., r = 4 w1*x11 + w2*x21 + w3*x21 + w4*x41 == r # col 1 w1*x12 + w2*x22 + w3*x32 + w4*x42 == r # col 2 w1*x13 + w2*x23 + w3*x33 + w4*x43 == r # col 3 # 2. By column, the number of X's greater than zero should be two or greater. In this small example, this condition is satisfied, but I would like to make sure that it is also satisfied in my problem. # 3....

...expression of the form: model1<-nls(y~beta1*(x1+(k1*x2)+(k1*k1*x3)+(k2*x4)+(k2*k1*x5)+(k2*k2*x6)+(k3*x7)+(k3*k4*x8)+(k3*k2*x9)+(k3*k3*x10)+ (k4*x11)+(k4*k1*x12)+(k4*k2*x13)+(k4*k3*x14)+(k4*k4*x15)+(k5*x16)+(k5*k1*x17)+(k5*k2*x18)+(k5*k3*x19)+ (k5*k4*x20)+(k5*k5*x21)+(k6*x22)+(k6*k1*x23)+(k6*k2*x24)+(k6*k3*x25)+(k6*k4*x26)+(k6*k5*x27)+(k6*k6*x28)+ (k7*x29)+(k7*k1*x30)+(k7*k2*x31)+(k7*k3*x32)+(k7*k4*x33)+(k7*k5*x34)+(k7*k6*x35)+(k7*k7*x36)),x, start=c(beta1=-0.001480981,k1=0.001,k2=0.001,k3=0.001,k4=0.001,k5=0.001,k6=0.001,k7=0.001),control=c(max...

...rg>> wrote: > > Hi, > > Checking the llvm test-suite SingleSource/Regression/C++/EH/class_hierarchy > testcase on aarch64 I noted something strange: > > Dump of assembler code for function _Z4funcj: > 0x0000000000400650 <+0>: stp x22, x21, [sp,#-48]! > 0x0000000000400654 <+4>: stp x20, x19, [sp,#16] > 0x0000000000400658 <+8>: stp x29, x30, [sp,#32] > 0x000000000040065c <+12>: add x29, sp, #0x20 > => 0x0000000000400660 <+16>: .inst 0x2bff...