search for: x14

...grams, 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+C17)A18*(X18*x1+X28*x2)+1*sqrt(B18*(X18*x1+X28*x2)^2+C18)A19*(X19*x1+X29*x2)+-1*sqrt(B19*(X19*x1+X2...

...,1),x3=c(5,3,4,1,5,5,5,5,5,5),x4=c(5,3,4,5,5,5,5,5,5,5),x5=c(5,4,3,3,4,4,4,5,NA,5),x6=c(3,5,2,1,4,NA,NA,5,3,4),x7=c(5,4,3,4,4,3,4,NA,NA,4),x8=c(5,5,3,4,5,4,5,4,5,5),x9=c(5,1,1,1,NA,3,3,2,2,3),x10=c(5,5,2,3,4,3,1,3,2,5),x11=c(5,4,4,5,5,4,5,4,5,5),x12=c(5,4,5,5,5,4,5,5,5,5),x13=c(5,4,4,5,5,5,5,5,5,5),x14=c(5,3,3,NA,5,4,5,3,NA,5),x15=c(4,4,4,3,4,5,2,4,5,5),x16=c(5,4,3,3,4,4,3,3,4,5)) and I want to carry out a confirmatory factor analysis using the "cfa" function (lavaan). library(lavaan) then I specify my model cfa.model<-'      f1=~x1+x2+x3+x4+x5                 f2=~x6+x7+x8+x...

Hi All, I'm attempting to use eps-regression or nu-regression SVM to compute probabilities but the predict function applied to an svm model object returns values outside [0, 1]: Variable Data looks like: Present X02 X03 X05 X06 X07 X13 X14 X15 X18 1 0 1634 48 2245.469 -1122.0750 3367.544 11105.013 2017.306 40 23227 2 0 1402 40 2611.519 -811.2500 3422.769 10499.425 1800.475 40 13822 3 0 1379 40 2576.150 -842.8500 3419.000 10166.237 2328.756 37 14200 4 0 1869 51 2645.794 -982.2938 3628.088 9610.037 1699.656 43 20762 ... and bgEnv look...

...olumns need to be ordered. The first column X is at the right position, but the remaining columns X1-Xn should be ordered like this: X1, X2, X3 etc instead of like below. > colnames(pos1) [1] "X" "X1" "X10" "X11" "X12" "X13" "X14" "X15" "X16" "X17" "X18" "X19" "X2" "X20" "X3" "X4" "X5" "X6" "X7" "X8" "X9" > pos1[1:5,1:5] X X1 X10 X11 X12 1...

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

Hi, May be this helps (Assuming that there are only '0's and '1's in the dataset) dat1<-read.table(text=" ??????? ID X0 X1 X2 X3 X4 X5 X6 X7 X8 X9 X10 X11 X12 X13 X14 X15 1?? 5184??? 0??? 0??? 0??? 0??? 0?? 0?? 0??? 0??? 0??? 1???? 0????? 0????? 0????? 0????? 0????? 0 2?? 6884??? 0??? 0??? 1??? 0??? 0?? 1?? 0??? 0??? 0??? 0???? 0????? 0????? 0????? 0????? 0????? 0 3? 11329?? 0??? 0??? 0??? 0??? 0?? 1?? 0??? 0??? 0??? 0???? 0????? 0????? 0????? 0????? 0????? 0 4?...

..., X11 = c(295.869476318359, 297.010070800781, 299.330169677734, 297.627593994141, 297.964935302734), X12 = c(295.986236572266, 297.675567626953, 299.056671142578, 297.598907470703, 298.481842041016 ), X13 = c(295.947601318359, 297.934448242188, 298.745391845703, 297.704925537109, 298.819091796875), X14 = c(294.654327392578, 296.722717285156, 297.0986328125, 296.508239746094, 297.822021484375 ), X15 = c(292.176361083984, 294.49658203125, 293.888305664062, 294.172149658203, 296.117095947266), X16 = c(288.400726318359, 291.029113769531, 289.361907958984, 290.566772460938, 293.554016113281 ), X17 = c...

How about: x <- data.frame(matrix(rnorm(1550),c(50,31))) model <- step(lm(x[,1] ~ as.matrix(x[,2:31]))) --Matt -----Original Message----- From: r-help-bounces at stat.math.ethz.ch [mailto:r-help-bounces at stat.math.ethz.ch]On Behalf Of rongguiwong Sent: Monday, September 20, 2004 20:52 PM To: r-help at stat.math.ethz.ch Subject: [R] how to take this experiment with R? This message uses

Dear all, I had a R segmentation fault, and then invoked debug mode and ran step by step. When I reached "terms(Y~X1*X2*...*X16)", I would then have "segmentation" fault. However, if I just ran this under regular "R interactive" mode, it would be fine though taking long time. My questions are: 1. Is there a known limit of terms for a formula? 2. Why does the

Dear all, I had a R segmentation fault, and then invoked debug mode and ran step by step. When I reached "terms(Y~X1*X2*...*X16)", I would then have "segmentation" fault. However, if I just ran this under regular "R interactive" mode, it would be fine though taking long time. My questions are: 1. Is there a known limit of terms for a formula? 2. Why does the

...9; [ldap://localhost] 'DC=@,DC=filial1.company.local,CN=MicrosoftDNS,DC=DomainDnsZones,DC=filial1,DC=company,DC=local' [ldap://dc02] ??? Difference in attribute values: ??????? dnsRecord => ['\x04\x00\x01\x00\x05\xf0\x00\x00:\x00\x00\x00\x00\x00\x02X\x00\x00\x00\x00\xa8k7\x00\xc0\xa8\x14\x01', '\x04\x00\x01\x00\x05\xf0\x00\x00:\x00\x00\x00\x00\x00\x03\x84\x00\x00\x00\x00\xb4k7\x00\xc0\xa8\x14\x02', '\x1e\x00\x02\x00\x05\xf0\x00\x00:\x00\x00\x00\x00\x00\x0e\x10\x00\x00\x00\x00\x00\x00\x00\x00\x1c\x04\x04dc02\x07filial1\x07company\x05local\x00', '\x1e\x00\x02\...

Hi everybody, If I've a data frame like this: dataframe a X0 X2 X4 X6 X8 X10 X12 X14 X16 1957 0 0 0 0 0 0 0 0 0 1958 0 0 0 0 0 0 0 0 0 1959 831 0 0 0 0 0 0 0 0 1960 544 282 0 0 0 0 0 0 0 1961 446 365 0 0 0 0 0 0 0 1962 442 473 0 0 0 0 0 0 0 1963 595 468 0 0 0 0 0 0 0 1964 400...

Hi, Suppose I've got a matrix, and the first few elements look like "x1 + x3 + x4 + x5 + x1:x3 + x1:x4" "x1 + x2 + x3 + x5 + x1:x2 + x1:x5" "x1 + x3 + x4 + x5 + x1:x3 + x1:x5" and so on (have got terms from x1 ~ x14). If I want to replace all the x1 with i7, all x2 with i14, all x3 with i13, for example. Is there an easy way? I tried to put what I want to replace in a vector, like: repl = c("i7", "i14", "i13", "d2", "i8", "i5", "i6&...

...timate SE z P(>|z|) 0.539 0.556 0.969 0.333 AIC: 114.2012 > (occu2 <- occu(~x3 ~x1, umf)) Call: occu(formula = ~x3 ~ x1, data = umf) Occupancy: Estimate SE z P(>|z|) (Intercept) -1.729 1.11 -1.5505 0.121 x13 -0.020 1.16 -0.0172 0.986 x14 -0.163 1.35 -0.1207 0.904 Detection: Estimate SE z P(>|z|) (Intercept) -1.0434 1.0762 -0.969 0.332 x3 0.0168 0.0106 1.581 0.114 AIC: 117.5561 > (occu3 <- occu(~x3+x1 ~x1, umf)) Call: occu(formula = ~x3 + x1 ~ x1, data = umf) Occupancy...

...0781, 299.330169677734, > > 297.627593994141, 297.964935302734), X12 = c(295.986236572266, > > 297.675567626953, 299.056671142578, 297.598907470703, 298.481842041016 > > ), X13 = c(295.947601318359, 297.934448242188, 298.745391845703, > > 297.704925537109, 298.819091796875), X14 = c(294.654327392578, > > 296.722717285156, 297.0986328125, 296.508239746094, 297.822021484375 > > ), X15 = c(292.176361083984, 294.49658203125, 293.888305664062, > > 294.172149658203, 296.117095947266), X16 = c(288.400726318359, > > 291.029113769531, 289.361907958984, 290.5...

...,DC=DomainDnsZones,DC= > > filial1,DC=company,DC=local' [ldap://dc02] > > > Difference in attribute values: > > > dnsRecord => > > > ['\x04\x00\x01\x00\x05\xf0\x00\x00:\x00\x00\x00\x00\x00\x02X\x00\x00 > > \x00\x00\xa8k7\x00\xc0\xa8\x14\x01', > > '\x04\x00\x01\x00\x05\xf0\x00\x00:\x00\x00\x00\x00\x00\x03\x84\x00\x > > 00\x00\x00\xb4k7\x00\xc0\xa8\x14\x02', > > '\x1e\x00\x02\x00\x05\xf0\x00\x00:\x00\x00\x00\x00\x00\x0e\x10\x00\x > > 00\x00\x00\x00\x00\x00\x00\x1c\x04\x04dc02\x07filial1\x07com...

Hi.. i have an 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, s...

...s to be compared: 19 Comparing: 'DC=@,DC=_msdcs.mytld.at,CN=MicrosoftDNS,DC=ForestDnsZones,DC=mytld,DC=at' [ldap://dc] 'DC=@,DC=_msdcs.mytld.at,CN=MicrosoftDNS,DC=ForestDnsZones,DC=mytld,DC=at' [ldap://pre01svdeb03] Difference in attribute values: dnsRecord => ['\x14\x00\x02\x00\x05\xf0\x00\x00n\x00\x00\x00\x00\x00\x03\x84\x00\x00\x00\x00\x00\x00\x00\x00\x12\x03\x02dc\nmytld\x02at\x00', '\x1e\x00\x02\x00\x05\xf0\x00\x00n\x00\x00\x00\x00\x00\x03\x84\x00\x00\x00\x00\x00\x00\x00\x00\x1c\x03\x0cpre01svdeb03\nmytld\x02at\x00', 'D\x00\x06\x00\x05\xf0\...

...318359, 297.010070800781, 299.330169677734, > 297.627593994141, 297.964935302734), X12 = c(295.986236572266, > 297.675567626953, 299.056671142578, 297.598907470703, 298.481842041016 > ), X13 = c(295.947601318359, 297.934448242188, 298.745391845703, > 297.704925537109, 298.819091796875), X14 = c(294.654327392578, > 296.722717285156, 297.0986328125, 296.508239746094, 297.822021484375 > ), X15 = c(292.176361083984, 294.49658203125, 293.888305664062, > 294.172149658203, 296.117095947266), X16 = c(288.400726318359, > 291.029113769531, 289.361907958984, 290.566772460938, 293.554...

Dear experts, I have in all 14 matrices which stands for gene expression divergence and 14 matrices which stands for gene sequence divergence. I have tried joining them by using the concatanation function giving SequenceDivergence <- c(X1,X2,X3,X4,X5,X6,X7,X8,X9,X10,X11,X12,X13,X14) ExpressionDivergence <- c(Y1,Y2,Y3,Y4,Y5,Y6,Y7,Y8,Y9,Y10,Y11,Y12,Y13,Y14) where X1,X2..X14 are the expression matrices containing r-values and Y1,Y2..Y14 are the ones with patristic distances Now, I want to plot SequenceDivergence vs. Expression Divergence Tried doing that using plot (Sequ...