search for: delta2

a graph question. Thanks a lot in advance. I made two scatterplots on one graph (sigma vs. delta1, sigma vs. delta2) (20 observations of delta1, delta2 and corresponding sigma) the x-axis is sigma, the y-axis is either delta1 or delta2. I connected both scatterplots. To seperate them, one curves is a line with circles, the other curve is a line with squares on it. I want to make a notation either on the y-axis...

...and therefore the submatrices are of the form: Lambda Delta1 Delta1 ... Delta1 Delta1 Lambda Delta1 ... Delta1 v11=v22= ....... Delta1 ..... Lambda Delta2 Delta2 Delta2 ... Delta2 Delta2 Delta2 Delta2 ... Delta2 v12=v21= ....... Delta2 ..... Delta2 The elements of these submatrices depend only upon the three covariance parameters: the compound symmet...

...n. loess.aic <- function (x) { # extract values from loess object if (!(inherits(x,"loess"))) stop("Error: argument must be a loess object") span <- x$pars$span n <- x$n traceL <- x$trace.hat sigma2 <- sum( x$residuals^2 ) / (n-1) delta1 <- x$one.delta delta2 <- x$two.delta enp <- x$enp aicc <- log(sigma2) + 1 + 2* (2*(traceL+1)) / (n-traceL-2) aicc1<- n*log(sigma2) + n* ( (delta1/(delta2*(n+enp)))/(delta1^2/delta2)-2 ) gcv <- n*sigma2 / (n-traceL)^2 result <- list(span=span, aicc=aicc, aicc1=aicc1, gcv=gcv) return(result) }...

...Global Parameters ############ beeta=0.8 pq=10000 pf=10000 F=20 L=12600 theta=0.6 psale=0.6 mu=psale*(1-theta) alphah=0.15 Cg=6240 Cs=2820 A= 100 D=0.0001 greekp=0.43 K=100000 ##### Species Parameters########## b1=0.38 p1=16654 v1 = 0.28 N1=6000 g1=1 delta1=1 b2=0.4 p2=2797 v2 = 0.31 N2=10000 g2=1 delta2=1 ####################################### Objective function ############################################ eval_f = function (x) { return(-1 * ( ( (1-alphah) *log(F) ) + ( alphah* ( (x[1]*(((N1/A)*(g1^greekp)*(x[3]^b1))+((1-exp(-2*D*v1*N1))*x[4]))) + (x[2]*(((N2/A)*(g2^greekp)*(x[3]^b2))+...

Hello, I am running Dovecot 1.0.5 and seem to have lots of errors in my Dovecot logs. The worst errors are things like: Oct 30 16:11:42 delta2 dovecot: IMAP(x): Corrupted transaction log file /home/spamcop-net2/deputies/dovecot.index.log: end_offset (392) > current sync_offset (388) Oct 30 16:11:42 delta2 dovecot: IMAP(x): file mail-index-sync-update.c: line 841 (mail_index_sync_update_index): assertion failed: (view->hdr.messag...

Hi r-users;   I have this values: eign_val <- c(137.810447,3.538721,2.995161,1.685670) alp    <- 1.6549 ;  lamda <- eign_val lamda_m <- min(lamda)   First I calculated manually: delta0 <- 1 delta1 <- alp*delta0*(4-lamda_m*(1/lamda[1]+1/lamda[2]+1/lamda[3]+1/lamda[4]))  delta1 delta2 <- (alp/2)*(delta1*(delta1/alp) + delta0*((1-lamda_m/lamda[1])^2+ (1-lamda_m/lamda[2])^2+(1-lamda_m/lamda[3])^2+(1-lamda_m/lamda[4])^2)) delta2 delta3 <- (alp/3)*(delta2*(delta1/alp) + delta1*((1-lamda_m/lamda[1])^2+ (1-lamda_m/lamda[2])^2+(1-lamda_m/lamda[3])^2+(1-lamda_m/lamda[4])^2)  + del...

...v) ###### Global Parameters ############ beeta=0.8 pq=10000 L=12600 theta=0.6 psale=0.6 mu=psale*(1-theta) alphah=0.15 Cg=6240 Cs=2820 A= 100 D=0.0001 greekp=0.43 K=100000 ##### Species Parameters ########## b1=0.38 p1=16654 v1 = 0.28 N1=6000 g1=1 delta1=1 b2=0.4 p2=2797 v2 = 0.31 N2=10000 g2=1 delta2=1 ### Define functions with vector x = c(Lg, Ls, gamma1, gamma2, lamda) firstordercond <- function (x) { y=numeric(4) y[1]=(alphah/x[3])-(x[5]*((p1-(((theta+mu)*(((N1/A)*g1^greekp*x[1]^b1)+K))+((theta+mu)*(((1-exp(-2*D*v1*N1))*x[2])+K))))*(((N1/A)*g1^(greekp))*x[1]^b1+(2*v1*N1*D)*x[2])...

This is an OpenPGP/MIME signed message (RFC 2440 and 3156) --------------enig807B2312A20EAF60129FDDFA Content-Type: text/plain; charset=ISO-8859-1 Content-Transfer-Encoding: quoted-printable Hi, I am writing this email, because I am not sure if the issue I have discovered is a bug or not. For a few days I have been fiddling around with a small program that calculates the reflectance of

Hi all, I'm not familiar with R programming and I'm trying to reproduce a result from a paper. Basically, I have a dataset which I would like to model in terms of successive increments, i.e. (y denote empirical values of y) y_1 = y1, y_2 = y1 + delta1, y_3 = y1 + delta1 + delta2. ... y_m = y1 + sum_2^m delta j where delta_j donote successive increments in the y-values, i.e. delta j = y_j - y_(j-1). In order to estimate y-values, I'm assuming that delta j is approximately equal to kj**u, such that my regression model should be something like this: ^y_1 = a1 ^y_2...

...erience in maximization could tell me if a function like that could be maximized using R and if s/he could point me to a package which might be useful for that. Just to give you some more details on the loglikelihood function attached: * Data vectors:u, zr, s * Indicator variables: epsilon, delta1, delta2 * Matrix with indicator variables: alpha * Unknown parameters: mu, beta, sigma * Unknown parameter vector: omega (contains probabilities). I'd be very grateful for any helpful information. Klaus. -- ---------------------------------------------------- Klaus Langohr Departament d'Estad...

In order to implement a subtle memory access optimisation during post-RA scheduling, I want to be able to determine some properties about the memory access. If I have two registers referring to memory, how can I determine if they are derived from the same base-pointer? Often LLVM will optimise to use intermediate registers holding partial displacements, for example, when a 'struct'

...sm2.lf <- likfit.glsm(gdn.glsm2.prelf, cov.model = > "exponential", ini.phi=26, lambda=0) > -------------------------------------------------------------------- > likfit.glsm: likelihood maximisation using the function optim. > phi = 26 tausq.rel = 0 > Error in if (det(Delta2) != 0) { : missing value where TRUE/FALSE needed > In addition: Warning message: > cannot use argument lambda with the given objects in mcmc.obj in: > likfit.glsm(gdn.glsm2.prelf, cov.model = "exponential", ini.phi = 26, ------------- Regarding the warning: the argument 'lamb...

...as.integer(sum.drop.sqr), as.double(span*cell), as.character(surf.stat), fitted.values = double(N), parameter = integer(7), a = integer(max.kd), xi = double(max.kd), vert = double(2*D), vval = double((D+1)*max.kd), diagonal = double(N), trL = double(1), delta1 = double(1), delta2 = double(1), as.integer(surf.stat == "interpolate/exact")) if(j==1) { trace.hat.out <- z$trL one.delta <- z$delta1 two.delta <- z$delta2 } fitted.residuals <- y - z$fitted.values if(j < iterations) robust <- .Fortran(R_lowesw, as.doubl...

....character(surf.stat), fitted.values = double(N), parameter = integer(7), a = integer(max.kd), xi = double(max.kd), vert = double(2 * D), vval = double((D + 1) * max.kd), diagonal = double(N), trL = double(1), delta1 = double(1), delta2 = double(1), as.integer(surf.stat == "interpolate/exact")) After that, I'm kind of stuck in terms of tracking it down. Thanks for any input, Jason Law City of Portland, OR

...ath.ethz.ch Hi all, I'm not familiar with R programming and I'm trying to reproduce a result from a paper. Basically, I have a dataset which I would like to model in terms of successive increments, i.e. (y denote empirical values of y) y_1 = y1, y_2 = y1 + delta1, y_3 = y1 + delta1 + delta2. ... y_m = y1 + sum_2^m delta j where delta_j donote successive increments in the y-values, i.e. delta j = y_j - y_(j-1). In order to estimate y-values, I'm assuming that delta j is approximately equal to kj**u, such that my regression model should be something like this: ^y_1 = a1 ^y_2...

...he data frame with column names var1 var2 Unit time trt (trt represents the treatment) Now suppose I want to fit a joint model of the form for the *i* th unit var1jk = alpha1 + beta1*timejk + gamma1* trtjk + delta1* timejk:trtjk + error1jk var2 = alpha2 + beta2*timejk + gamma2* trtjk + delta2* timejk:trtjk + error2jk where j index time and k index the treatment received Using indicator variables I have been able to fit and run the code for a bivariate model using unstructured covariance matrix. However, I want to fit a model for a structured variance covariance matrix. The error struc...

Hi everybody, I have joined 3 of my 40 samba fileservers on my active directory domain. And today I have a weird error message while I am trying to join another fileserver. Error message is the following : # tail -f /var/log/samba/log.winbindd-idmap [2010/07/16 13:26:02, 1] winbindd/idmap.c:idmap_init(385) Initializing idmap domains [2010/07/16 13:27:47, 1]