search for: eta2

Hi - I am trying to fit sem model with intercepts. Here is what I have in my model. Exogeneous vars: x1 (continous), x2 (ordinal), x3 (ordinal), x4(continuous) Endogeneous vars: y1 (continuous), y2 (ordinal), y3 (ordinal) SEM model: x1 -> eta1; x2 -> eta1; x3 -> eta2; x4 -> eta2; eta1 -> y1, eta1 -> y2, eta2 -> y2, eta2 -> y3 However, in these arrow models, I don't know how to add intercept onto it. I am trying to find an example code using sem package on how to incorporate intercept but cannot find any documents on the web. Or we can simpl...

...code. Many thanks in advance. Sincerely, Loki ################################################################# SAS Code: proc nlmixed data=repdat parms b0 -3 b1 -.135 a0 3 a1 4 sigsq 0.25; if vs.flag = 1 then do; eta1 = b0 + b1*lnbldT; llbin = anybc.cens.ind*eta1 - log(1+exp(eta1)); eta2 = a0 + a1*lnndsTs; llnorm = -1/(2*sigsq)*(lnbldT - eta2)**2 - .5*log(sigsq); ll = llbin + llnorm; end; else do; eta2 = a0 + a1*lnndsTs; eta1 = b0 + b1*eta2 + u; llbin = anybc.cens.ind*eta1 - log(1+exp(eta1)); ll = llbin; end; sigma2 = sigmasq*b1**2 /*variance of ra...

...the profile will sometimes work with this) # parameter assignments Nmin <- 0.8852 Nopt1 <- 16.78 gN <- 0.5511 E.nfert1 <- 0.3271 E.nfert2 <- 0.6132 Beta <- 0.8902 Dls <- 0.5378 eta1 <- 0.3791 eta2 <- 0.6332 PopStd <- 90468 beta <- Beta DIs <- Dls MnmN <- Nmin OptN <- Nopt1 # define model function Y.model <- function(gN, MnmN, OptN, DIs, beta, eta1, eta2, Popn, Dmax, AWC, SumEp, PotYield3, Nsupply)...

...der to compare this model with the same specification estimated with Ols Can you please help me or give me an hint? If it can be useful, i hereby put an estimation with stata's oprobit function (varnames are a little bit different) > . oprobit lz2_1 ly lnfiglimin lnfiglimagg lnadulti eta eta2 uomo elem laurea saluteno extra > Iteration 0: log likelihood = -7877.9258 > Iteration 1: log likelihood = -7633.4001 > Iteration 2: log likelihood = -7633.3455 > Ordered probit estimates Number of obs = 3747 >...

..., 4.7) ;CL (0, 38) ;V (0.0, .5) ; EMAX (0.0, 100) ; EC50 (0.0, 1, 10) ; N $OMEGA 0.04 0.04 0.04 0.04 0.04 0.04 $SIGMA 0.04 5 $ESTIMATION MAXEVAL=9999 NOABORT METHOD=0 POSTHOC PRINT=5 $COVARIANCE $TABLE ID TIME AGE BW GNDR TVKA TVV TVCL ETA1 ETA2 ETA3 EVID NOPRINT ONEHEADER FILE=TABLE.PAR $TABLE NOPRINT FILE=../3033.TAB ONEHEADER ID TIME AGE BW GNDR KA V CL TVKA TVCL TVV FLAG FLA2 EVID EFF CONC IPRED FLAG $TABLE NOPRINT ONEHEADER FILE=../3033par.TAB ID TIME AGE BW GNDR KA V CL TVKA TVCL TVV EMAX EC50 ETA1 ETA2 ETA3 EVID $TABLE ID TIME COM(...

Hello everybody After doing a MANOVA on a bunch of data, I want to be able to make some comment on the amount of variation in the data that is explained by the factor of interest. I want to say this in the following way: XX% of the data is explained by A. I can acheive something like what I want by doing the following: X <- structure(c(9, 6, 9, 3, 2, 7), .Dim = as.integer(c(3,

...two problems, related to the function sub, grep, regexpr and similia. The header of the file(s) I have to import is like this. c("y (m)", "BD (g/cm3)", "PR (Mpa)", "Ks (m/s)", "SP g./g.", "P (m3/m3)", "theta1 (g/g)", "theta2 (g/g)", "AWC (g/g)") To get rid of spaces and symbols in the names of the columns, I use read.table(... check.names=TRUE) and I get: str <- c("y..m.", "BD..g.cm3.", "PR..Mpa.", "Ks..m.s.", "SP.g..g.", "P..m3.m3.", &quot...

Hi, Having searched google '[R] aov effect size' without any results I wonder if I not completely miss something. Is there any R function that calculates the effect size of an AOV's main effect or interaction effect? It should be related to the F's and the degree of freedom of the error, but the multitude in numbers in aov() baffle me a bit. Thanks, ---david

...test.invisible[2] <- match("quali.sup", invisible) } else test.invisible <- rep(NA, 2) if (new.plot) dev.new() if (is.null(title)) title <- "Variables representation" coord.actif <- res.mca$var$eta2[, axes] if (!is.null(res.mca$quali.sup$eta2)) coord.illu <- res.mca$quali.sup$eta2[, axes, drop = FALSE] if (is.na(test.invisible[1])) { plot(coord.actif, xlab = lab.x, ylab = lab.y, xlim = c(0, 1), ylim = c(0, 1), pch = 20, col = col.var...

Hi, all, Does anybody got the most general expression of a unitary matrix? I found one in the book, four entries of the matrix are: (cos\theta) exp(j\alpha); -(sin\theta)exp(j(\alpha-\Omega)); (sin\theta)exp(j(\beta+\Omega)); (cos\theta) exp(j\beta); where "j" is for complex. However, since for any two unitary matrices, their product should also be a unitary matrix. When I

[This email is either empty or too large to be displayed at this time]