search for: b21

Is there an element-by-element multiplication in R, like the .* operator in Matlab? eg: A (2x3) B (2x3) C=A.*B C (2x3) C = [[a11*b11 a12*b12 a13*b13]; [a21*b21 a22*b22 a23*b23]] I can't find one... Thanks -Mike Gates

....30,2416.20,2484.80,2608.50,2744.10,2729.30,2695.00,2826.70,2958.60,3115.20,3192.40,3187.10,3248.80,3166.00,3277.70,3492.00,3573.50) money<-data.frame(r,M,Y) attach(money) ols1<-lm(log(M)~log(r)+log(Y)) output1<-summary(ols1) coef1<-ols1$coefficients a1<-coef1[[1]] b11<-coef1[[2]] b21<-coef1[[3]] money.m1<-nls(log(M)~a+b*r^g+c*Y^g,data=money,start=list(a=a1,b=b11,g=1,c=b21)) summary(money.m1) money.m2<-boxcox(money.m1) Prof. Ikerne del Valle Erkiaga Department of Applied Economics V Faculty of Economic and Business Sciences University of the Basque Country A...

...Panel. R is installed from my limited user account. The version of the JVM I have installed is perhaps best described as antique: > system(paste("java -version"),show.output.on.console=T) java version "1.4.1" Java(TM) 2 Runtime Environment, Standard Edition (build 1.4.1-b21) Java HotSpot(TM) Client VM (build 1.4.1-b21, mixed mode) The HTML help search applet has worked for all versions of R up to and including 2.6.0, but in 2.6.1 the java applet is not initialised ("Applet SearchEngine notinited"). I have checked Appendix D of the admin & installat...

...email: c.fezzi at uea.ac.uk *************************************** Here is an example with 3 equations and 2 exogenous variables: ----- start code ------ N <- 1000 # number of observations library(MASS) ## parameters ## # eq. 1 b10 <- 7; b11 <- 2; b12 <- -1 # eq. 2 b20 <- 5; b21 <- -2; b22 <- 1 # eq.3 b30 <- 1; b31 <- 5; b32 <- 2 # exogenous variables x1 <- runif(min=-10,max=10,N) x2 <- runif(min=-5,max=5,N) # residual covariance matrix sigma <- matrix(c(2,1,0.7,1,1.5,0.5,0.7,0.5,2),3,3) # residuals r <- mvrnorm(N,mu=rep(0,3), Sigma=sigma)...

...s your domain name? If you have special > characters in the domain name > the workstations won't join properly. The domain > name should be all one > word. > > Kevin B. McCrory > Network Engineer - COPS > US Government Solutions > 13600 EDS Drive > Mail stop: A4S-B21 > Herndon, VA 20171 > * phone: +01-703-733-3255 > * mailto:kevin.mccrory@eds.com > * AKO mailto:kevin.mccrory@us.army.mil > > > > -----Original Message----- > From: > samba-bounces+kevin.mccrory=eds.com@lists.samba.org > [mailto:samba-bounces+kevin.mccrory=eds.co...

...tc/smbldap-tools/ directory, then edit the /etc/smbldap-tools/smbldap_conf.pm" instead of "CHange to the /opt/IDEALX/sbin/ directory then edit the /opt/IDEALX/sbin/smbldap_conf.pm" . Kevin B. McCrory Network Engineer - COPS US Government Solutions 13600 EDS Drive Mail stop: A4S-B21 Herndon, VA 20171 * phone: +01-703-733-3255 * mailto:kevin.mccrory@eds.com <mailto:kevin.mccrory@eds.com> * AKO mailto:kevin.mccrory@us.army.mil <mailto:kevin.mccrory@us.army.mil>

...h is converted using as.character, looses its last part: ``diagonal = 1e-12)'' Shorter formula is ok though. Best, H??vard ************ Browse[2]> formula.str y ~ -1 + b1 + b2 + b3 + b4 + b5 + b6 + b7 + b8 + b9 + b10 + b11 + b12 + b13 + b14 + b15 + b16 + b17 + b18 + b19 + b20 + b21 + b22 + b23 + b24 + b25 + b26 + b27 + b28 + b29 + b30 + b31 + b32 + b33 + b34 + b35 + b36 + b37 + b38 + b39 + b40 + b41 + b42 + b43 + b44 + b45 + b46 + b47 + b48 + b49 + elevation + f(idx, model = "sphere", sphere.dir = "global_temperature_80s", T.or...

...he input of file has been a problem for me. Take this case for example (http://code.google.com/codejam/contest/dashboard?c=agdjb2RlamFtchALEghjb250ZXN0cxjRzBQM), the files are usually of the form: A(number of lines for group 1) a11 a12 a13 a21 a22 a23 ... B(number of lines for group 2) b11 b12 b13 b21 b22 b23 ... I guess SAS may work pretty well in this kind of situation with data step. But I don't know how to handle them using R. Any suggestions? Thanks a lot. Best wishes, -- ??? Hesen Peng http://hesen.peng.googlepages.com/

...6*(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+X29*x2)^2+C19)A110*(X110*x1+X210*x2)+1*sqrt(B110*(X110*x1+X210*x2)^2+C110)" > cum2 [1] "A21*(X11*x1+X21*x2)+1*sqrt(B21*(X11*x1+X21*x2)^2+C21)A22*(X12*x1+X22*x2)+1*sqrt(B22*(X12*x1+X22*x2)^2+C22)A23*(X13*x1+X23*x2)+-1*sqrt(B23*(X13*x1+X23*x2)^2+C23)A24*(X14*x1+X24*x2)+-1*sqrt(B24*(X14*x1+X24*x2)^2+C24)A25*(X15*x1+X25*x2)+1*sqrt(B25*(X15*x1+X25*x2)^2+C25)A26*(X16*x1+X26*x2)+1*sqrt(B26*(X16*x1+X26*x2)^2+C26)A27*(X17*x...

Hi, Is there any package that I can use to simulate the Gaussian Mixture Model , which is a mixture modeling method that is widely used in statistical learning theory. I know there is a mclust, however, I think it is a little bit different from my problem. Thanks very much.. regards. -------------------------- Peng Jiang ?? Ph.D. Candidate Antai College of Economics &

...d rechecked the Administrator account and made sure the password is correct for the account. Not sure where to go from here... Some pointers on what to look at would be greatly appreciated... Thanks Kevin B. McCrory Network Engineer - COPS US Government Solutions 13600 EDS Drive Mail stop: A4S-B21 Herndon, VA 20171 * phone: +01-703-733-3255 * mailto:kevin.mccrory@eds.com <mailto:kevin.mccrory@eds.com> * AKO mailto:kevin.mccrory@us.army.mil <mailto:kevin.mccrory@us.army.mil>

...r. Most of the examples we cover in class are change models, i.e., we working with longitudinal data. Specific questions: in HLM 6.0, we build the following model; Y = P0 + P1*(CONFLICT) + P2*(TIMEYRS) + E Level-2 Model P0 = B00 + B01*(H0MCITOT) + R0 P1 = B10 + B11*(H0MCITOT) + R1 P2 = B20 + B21*(H0MCITOT) + R2 Can someone explain to me how to represent this in lmer syntax? I've tried e.g., lmer(MAT ~ 1 + CONFLICT + TIMEYRS + (1 + CONFLICT + + TIMEYRS | H0MCITOT)) But I don't get the same result. More generally: Should I be using the lme4 package, the nlme package, or somet...

...d> Blackdown Project was available, but it was written against the vulnerable JRE 1.4.2. I too gave up on this until I discovered Red Hat's "iced tea" project. The latest incantation for Fedora Core 8 loaded well on a RHEL5 box using rpm's: java-1.7.0-icedtea-1.7.0.0-0.19.b21.snapshot.fc8.x86_64.rpm java-1.7.0-icedtea-plugin-1.7.0.0-0.19.b21.snapshot.fc8.x86_64.rpm tzdata-java-2007h-1.fc8.noarch.rpm Then do: alternatives --config java and select the new 1.7 JRE. Usually FC8 stuff doesn't install well on RHEL5/Centos5/SAL5 boxes. However, these worked like a char...

..."category" "a11" "a12" "a13" "a21" [7] "a22" "a23" "a31" "a32" "b11" "b12" [13] "b13" "b21" "b31" "b32" "b33" "b41" [19] "b42" "c11" "c12" "c21" "c22" "c23" [25] "c31" "c32"...

...#################################################### # Part 1: calibration library(e1071) #calibration step calibrate<-read.table("calibration.txt", header=TRUE) calibrate$calibration<-as.factor(calibrate$calibration) calibrate.rf<-svm(calibration~B1+B14+B15+B16+B17+B18+B19+B20+B21+B24+B25+B26+B51+B52+B53+B54+B55+B56+B57+B58+B59+B60+B61+B62, data=calibrate, cost=6.8, gamma=0.08) #################################################################################################################### ##################################################################################...

...0) b10 <- c(0,0,0,0,0,0,0) b11 <- c(0,0,0,0,0,0,0) b12 <- c(0,0,0,0,0,0,0) b13 <- c(0,0,0,0,0,0,1) b14 <- c(0,0,0,0,0,0,1) b15 <- c(0,0,0,0,0,1,0) b16 <- c(0,0,0,0,1,0,0) b17 <- c(0,0,0,1,0,0,0) b18 <- c(0,0,1,0,0,0,0) b19 <- c(0,1,0,0,0,1,0) b20 <- c(0,1,0,0,0,1,0) b21 <- c(0,0,0,1,0,1,1) b22 <- c(1,0,0,1,0,1,1) b23 <- c(1,0,1,1,0,1,0) bb <- rbind (b1,b2,b3,b4,b5,b6,b7,b8,b9,b10,b11,b12,b13,b14,b15,b16,b17,b18,b19,b20, b21,b22,b23) a <- array(0, c(161,3), list(NULL,c("Subject","Item","Response"))) for (s in c(1:23)...

...- c(0,0,0,0,0,0,0) b11 <- c(0,0,0,0,0,0,0) b12 <- c(0,0,0,0,0,0,0) b13 <- c(0,0,0,0,0,0,1) b14 <- c(0,0,0,0,0,0,1) b15 <- c(0,0,0,0,0,1,0) b16 <- c(0,0,0,0,1,0,0) b17 <- c(0,0,0,1,0,0,0) b18 <- c(0,0,1,0,0,0,0) b19 <- c(0,1,0,0,0,1,0) # differs b20 <- c(0,1,0,0,0,1,0) b21 <- c(0,0,0,1,0,1,1) b22 <- c(1,0,0,1,0,1,1) b23 <- c(1,0,1,1,0,1,0) bb <- rbind(b1,b2,b3,b4,b5,b6,b7,b8,b9,b10,b11,b12,b13,b14,b15,b16,b17,b18,b19,b20,b21,b22,b23) a <- array(0, c(161,3), list(NULL, c("Subject","Item","Response"))) for (s in c(1:23))...

...0) b10 <- c(0,0,0,0,0,0,0) b11 <- c(0,0,0,0,0,0,0) b12 <- c(0,0,0,0,0,0,0) b13 <- c(0,0,0,0,0,0,1) b14 <- c(0,0,0,0,0,0,1) b15 <- c(0,0,0,0,0,1,0) b16 <- c(0,0,0,0,1,0,0) b17 <- c(0,0,0,1,0,0,0) b18 <- c(0,0,1,0,0,0,0) b19 <- c(0,1,0,0,0,1,0) b20 <- c(0,1,0,0,0,1,0) b21 <- c(0,0,0,1,0,1,1) b22 <- c(1,0,0,1,0,1,1) b23 <- c(1,0,1,1,0,1,0) bb <- rbind (b1,b2,b3,b4,b5,b6,b7,b8,b9,b10,b11,b12,b13,b14,b15,b16,b17,b18,b19,b20, b21,b22,b23) a <- array(0, c(161,3), list(NULL,c("Subject","Item","Response"))) for (s in c(1:23))...

...ants of the diagonal blocks. If a diagonal block is of dimension 1x1, then the determinant of the block is simply the value of the single element in the block. If the diagonal block is of dimension 2x2 then the determinant of the block may be computed according to the well-known formula b11*b22-b12*b21, where bij is the value in the i'th row and j'th column of the block. int i, q, info, *ipiv, one = 1; double *b, *A, *D, det; /* ** A and D are upper triangular matrices in packed storage ** A[] = a00, a01, a11, a02, a12, a22, a03, a13, a23, a33, ... ** use the following macr...

To whom it may concern, I am a student at Macaleste College, and next semester Macalester is going to offer a course for CellBio that is mainly statistically based. For the most part the students will be using R for analysis. The problem is there will be some simple differential equations for the students to solve. The committee that in charge of the classes corriculam would like only to