search for: xi

...you reply. Thanks! Chris JAGS model ... model { for(i in 1:N){ ########################## # measurement equation model ########################## for(j in 1:P){ y[i,j]~dnorm(mu[i,j],psi[j]) } ########################## # alp[i] corresponds to the intercepts ########################## mu[i,1]<-xi[i,1]+alp[1] ## Ext1 mu[i,2]<-lam[1]*xi[i,1]+alp[2] mu[i,3]<-lam[2]*xi[i,1]+alp[3] mu[i,4]<-xi[i,2]+alp[4] ## Soc1 mu[i,5]<-lam[3]*xi[i,2]+alp[5] mu[i,6]<-lam[4]*xi[i,2]+alp[6] mu[i,7]<-xi[i,3]+alp[7] ## Aca1 mu[i,8]<-lam[5]*xi[i,3]+alp[8] mu[i,9]<-lam[6]*xi[i,3]+alp[9] mu[i...

...e PASEI 6: In read.spss("gss-pspp.sav", to.data.frame = T) : gss-pspp.sav: File contains duplicate label for value 99.9 for variable MASEI 7: In read.spss("gss-pspp.sav", to.data.frame = T) : gss-pspp.sav: File contains duplicate label for value 99.9 for variable SPSEI 8: In xi >= z[1L] | xi <= z[2L] | xi[xi == z[3L]] : longer object length is not a multiple of shorter object length 9: In xi >= z[1L] | xi <= z[2L] | xi[xi == z[3L]] : longer object length is not a multiple of shorter object length 10: In xi >= z[1L] | xi <= z[2L] | xi[xi == z[3L]] :...

...groups of variables that are all in one dataset. This is the code: #Stand counter st<-1 #Collections stands<-numeric(67) slopes<-numeric(67) intercepts<-numeric(67) mses<-numeric(67) rsquares<-numeric(67) pValues<-numeric(67) #Start lists for X and Y values within each stand xi<-numeric(0) yi<-numeric(0) #Set the first element to the starting X and Y values xi[1]=X[1] yi[1]=Y[1] #Start looping working through your data, record by record for (i in 2:length(X)) { #If you are in the same stand as on the last record, continue to #collect X and Y values if(Stand[i]...

Hi, I am trying to fit a Johnson SB distribution using fitdist function in fitdistrplus Library. I have defined the Johnson SB distribution from ( http://www.ntrand.com/johnson-sb-distribution/) . But it gives me the follwing errors. Any help would be appreciated #xi = xi #lambda =l #delta =d #gamma = g djohn = function(x,xi,l,d,g) (d/(l*sqrt(2*pi)*((x-xi)/l)*(1-((x-xi)/l))))*exp[-0.5*(g + d*log(((x-xi)/l)/(1-((x-xi)/l))))^2] pjohn = function(x,xi,l,d,g) pnorm(g + d*log(((x-xi)/l)/(1-((x-xi)/l)))) qjohn = function(p,xi,l,d,g) xi + (l*exp((qnorm(p) - g)/d...

...ted, right? The real problems come out when I try to fit a General Pareto distribution. Using the same method, I calculated the negative log likelihood function for the General Pareto distribution: > # negative log-likelihood function for general pareto distribution > nll_gpd <- function(xi, log_sigma, mu) { # shape, scale, location + n <- length(interarrival_times) + x <- interarrival_times + cat("xi:",xi,"; log_sigma:",log_sigma,"; mu:",mu,"\n") + n*log_sigma + (xi+1)*sum(log(1+xi*(x-mu)/exp(log_sigma)))/x...

...else if (R_IsNA(REAL(x)[i]) && R_IsNA(REAL(y)[j])) return 1; else if (R_IsNaN(REAL(x)[i]) && R_IsNaN(REAL(y)[j])) return 1; else return 0; } /* Between 1.34x and 1.37x faster on my 64-bit Ubuntu laptop */ static int requal2(SEXP x, int i, SEXP y, int j) { double xi, yj; if (i < 0 || j < 0) return 0; xi = REAL(x)[i]; yj = REAL(y)[j]; if (!ISNAN(xi) && !ISNAN(yj)) return xi == yj; if (R_IsNA(xi) && R_IsNA(yj)) return 1; if (R_IsNaN(xi) && R_IsNaN(yj)) return 1; return 0; } /* Another extra 1.18x...

I am assuming the variable out is the output parameter. However, I don't understand what is out[1+xi*xx>0]? Can someone explain this to me? Thanks in advance, Chen dGEV <- function(x, xi, mu = 0, sigma = 1, logvalue=FALSE) { xx <- (x-mu)/sigma #use the new dGumbel which passes mu and sigma: #if (xi==0) out <- dGumbel(xx,logvalue=TRUE)-log(sigma) if (xi==0) {...

Dear helpers, I am using R version 2.5.1 to estimate a multinomial logit model using my own maximum likelihood function (I work with share data and the default function of R cannot deal with that). However, the computer (I have an Athlon XP 3200+ with 512 GB ram) takes quite a while to estimate the model. With 3 categories, 5 explanatory variables and roughly 5000 observations it takes 2-3 m...

...= integer())) nms <- names(allargs) if (is.null(nms)) nms <- character(length(allargs)) cl <- NULL perm <- rows <- rlabs <- vector("list", n) nrow <- 0 value <- clabs <- NULL all.levs <- list() for (i in 1:n) { xi <- allargs[[i]] nmi <- nms[i] if (is.matrix(xi)) allargs[[i]] <- xi <- as.data.frame(xi) if (inherits(xi, "data.frame")) { if (is.null(cl)) cl <- oldClass(xi) ri <- attr(xi, "row.names&quo...

A post last week contained the design document for a shadow page table mechanism. This post contains more information, and a pointer to the code. As a recap, the "XI Shadow Mechanism" is a design for shadow page table code for fully virtualized HVM domains running on a 64-bit Xen hypervisor. This work was undertaken to address a number of goals. These are enumerated in the document posted last week and include: - ability to run fully virtualized 32, 32P...

Hi all, I have a list of 20000 data, and the list look like below. I wonder what is the simplest way to extract 'kappa' value (or 'xi' or 'alpha' for the matter) from each of the data. How can I simply code it without having to change the list to a dataframe first? Many thanks! $X19997 xi alpha kappa 784.7718640 165.4065141 -0.2709599 $X19998 xi alpha kappa 860.7832604 21...

...3 x 4 d 4 y 5 a 5 x 6 e 6 z > ## All is fine. ## Third data frame has a different third column variable > fifi <- data.frame(v1 = c('a', 'e'), v2 = c(9, 10), v4 = ordered(c('y', 'z'))) > rbind(foo, bar, fifi) Error in match.names(clabs, names(xi)) : names do not match previous names: v4 > ## Output looks fine. ## Third data frame is missing a column. > gaga <- data.frame(v1 = c('a', 'e'), v2 = c(7, 8)) > rbind(foo, bar, gaga) Error in paste(nmi[nii == 0], collapse = ", ") : object "nii" n...

Hi, The automatic package builder Xi has moved to system within the Dovecot infrastructure. This also means that the repositories are available at a new host name: xi.dovecot.fi. The old host name still works and is a CNAME pointing to the new venue. I expect the old domain to be active for years to come. So, there is no need to u...

Hi, The automatic package builder Xi has moved to system within the Dovecot infrastructure. This also means that the repositories are available at a new host name: xi.dovecot.fi. The old host name still works and is a CNAME pointing to the new venue. I expect the old domain to be active for years to come. So, there is no need to u...

...length of x until x1 is produced in the first of the loop. Thanks. Steven On 2/28/2024 6:22 PM, Eik Vettorazzi wrote: > Hi Steven, > It's not entirely clear what you actually want to achieve in the end. > > As soon as you "know" x1, and assuming that the different "xi" do not > differ in length in the real application, you know the length of the > target vector. > Instead of the loop, you can use 'Reduce' without having to initialize > a starting vector. > > # generate sample vectors, put them in a list > > xi<-lapply(1...

Dear all, I am stuck on the following problem with integrate(). I have been out of luck using RSiteSearch().. My function is g2<-function(b,theta,xi,yi,sigma2){ xi<-cbind(1,xi) eta<-drop(xi%*%theta) num<-exp((eta + rep(b,length(eta)))*yi) den<- 1 + exp(eta + rep(b,length(eta))) result=(num/den)*exp((-b^2)/sigma2)/sqrt(2*pi*sigma2) r=prod(result) return(r) } And I am interested...

Hi, For EAs data structure in inode/block are little different from them in bucket. These two patches try to make them same for the most part. The first patch set xh_free_start and xh_name_value_len when EAs in inode/block. xh_free_start is useful to keep the minimum offset of the xattr name/value. But xh_name_value_len is not very useful because we don't have "hole" when EAs in

Hi, I'm trying to customize a window with 2 graphs. I'm able to do the first one with something like this general example par(mfrow=c(1,2),cex.axis=0.85,cex.lab=0.80,mai=c(1.3,1,0.5,0),las=3) bplot<-barplot(bar.values,names.arg=cf.names,width=0.5,ylab="% Area held") abline(h=0.3,lty=3,col="red") abline(h=0.1,lty=3,col="blue") points(bplot,target$target/xi[nrow(xi),],pch=20) text(bplot,target$target/xi[nrow(xi)...

Hi, I have fixed the problems in version 1 patches. These two patches based on the latest main line kernel. Thanks, tiger > For EAs data structure in inode/block are little different from them in > bucket. These two patches try to make them same for the most part. > > The first patch set xh_free_start and xh_name_value_len when EAs in > inode/block. xh_free_start is useful to

Dear All, I am starting to use the tikzDevice package, but I am experiencing some (minor problems). If I try to run the script at the end of the email, I get the following error Error in tools::texi2dvi("normal.tex", pdf = T) : Running 'texi2dvi' on 'normal.tex' failed. Messages: sh: /usr/bin/texi2dvi: not found How do I install the missing texi2dvi? I found some reference to it at http://bit.ly/bSYbxk but I have not understood well whether it is should already be...