search for: nu

Hi All, Here 's what I got using dgamma function : > nu<-.2 > nu*log(nu)-log(gamma(nu))+(nu-1)*log(1)-nu*(1) [1] -2.045951 > dgamma(1,nu,nu,1) [1] 0.0801333 > dgamma(1,nu,nu,0) [1] NaN Warning message: NaNs produced in: dgamma(x, shape, scale, log) Could anyone tell me what is wrong here ? I am using R-2.4.1 on windows XP. Thanks a lot.

...t this was incredibly easy due to Rcpp.package.skeleton() and the countless answers to quesions online and documentation! Secondly, I'm worried that my speedup has been so massive (over 500x !!!) that I think I've made a mistake, hence my post here. Here is all my code, if someone has a minute to point out anything wrong (or even if its correct and there is room for improvement, im pretty new to this) it would be much appreciated. I've had a simple look at the results and they look fine, but seriously, 500x faster?! function in R; library(compiler) Rrcomp <- cmpfun( fu...

...14,18,3,4,34, 32,4,9,9,8,36,44,0,9,27,28,55,72,12,1, 9,0,32,0,0,2,15,5,6,17, 63,61,9,15,15,0,2) #########goal is  found unique maximum for  5-parameter function with using  "maxLik" function #####the function in latex commands is as following:  \begin{eqnarray*} &&\ell(\xi,\omega,\nu,\lambda_1,\lambda_2)=n\log2-n\log\omega+n\log\Gamma(\frac{\nu+1}{2})-\frac{n}{2}\log(\nu\pi)\\ &-&n\log\Gamma(\frac{\nu}{2})-\frac{\nu+1}{2}\sum_{i=1}^{n}\log(1+\frac{(x_i-\xi)^2}{\omega^2\nu})+\sum_{i=1}^{n}\log\Phi(\lambda_1\frac{(x_i-\xi)}{\sqrt{\omega^2+\lambda_2(x_i-\xi)^2}}) \end{eqna...

Here is a copy of the last few lines in base-Ex.Rout.fail: > x <- seq(3,500);yl <- c(-.3, .2) > plot(x,x, ylim = yl, ylab="",type='n', main = "Bessel Functions Y_nu(x)") > for(nu in nus){xx <- x[x > .6*nu]; lines(xx,besselY(xx,nu=nu), col = nu+2)} > legend(300,-.08, leg=paste("nu=",nus), col = nus+2, lwd=1) > > x <- seq(10,50000,by=10);yl <- c(-.1, .1) > plot(x,x, ylim = yl, ylab="",type='n', mai...

Dear List: The square of the noncentral t-statistic with noncentrality parameter \delta is a noncentral F with noncentrality parameter \lambda=\delta^2. So, t^2_{\nu,\delta} = F_{1,\nu,\lambda=\delta^2}. Consequently, it should follow that t^2_{1-\alpha/2,\nu,\delta} = f_{1-alpha,1,\vu,\lambda=\delta^2}. However, this is not what is happening with the following code. The central distributions agree as they should but the noncentral distributions do not....

I have no experience with incorporating Fortran code and am probably doing something pretty stupid. I want to use the following Fortran subroutine (not written by me) in the file SSFcoef.f subroutine SSFcoef(nmax,nu,A,nrowA,ncolA) implicit double precision(a-h,o-z) implicit integer (i-n) integer l,i,nmax double precision nu,A(0:nmax,0:nmax) A(0,0) = 1D0 do l=1,nmax do i=1,l-1 A(l,i) = (-nu+i+l-1D0)*A(l-1,i)+A(l-1,i-1) end do A(l,0) = (-...

I have no experience with incorporating Fortran code and am probably doing something pretty stupid. I want to use the following Fortran subroutine (not written by me) in the file SSFcoef.f subroutine SSFcoef(nmax,nu,A,nrowA,ncolA) implicit double precision(a-h,o-z) implicit integer (i-n) integer l,i,nmax double precision nu,A(0:nmax,0:nmax) A(0,0) = 1D0 do l=1,nmax do i=1,l-1 A(l,i) = (-nu+i+l-1D0)*A(l-1,i)+A(l-1,i-1) end do A(l,0) = (-...

Hi there, I had a serious problem here . Consider the following Bayesian model(discretized variance gamma): #Likelihood J[i]<-lambda*G[i]+sigma*sqrt(G[i])*rnorm(0,1) G[i]<-rgamma(1/nu,1/nu) #Prior: nu<-rinvgamma(m,M) # Parameters lambda=-.04 ; sigam=.38; nu=6.48; m=10,M=10 ; T=5000 (length of data) An author claimed that he got posterior distribution of nu with standard deviation .0965 But I could only get posterior sd=.84 with the same setting. I have been...

I am looking for a method of dealing with the modified Bessel function K_\nu(x) for large \nu. The besselK function implementation of this allows for dealing with large values of x by allowing for exponential scaling, but there is no facility for dealing with large \nu. What would work for me would be an lbesselK function in the manner of lgamma which returned the log...

#bug 1: besselI() for nu<0 and expon.scaled=TRUE #tested with R-devel (2007-06-17 r41981) x <- 2.3 nu <- -0.4 print(paste(besselI(x, nu, TRUE), "=", exp(-x)*besselI(x, nu, FALSE))) #fix: #$ diff bessel_i_old.c bessel_i_new.c #57c57 #< bessel_k(x, -alpha, expo) * ((ize == 1)? 2. : 2.*exp(-x))/M_PI #-...

When running Dovecot (rc6 CVS dovecot-20060808) as a proxy, why would it log the following when it has not yet maxed out the allowed number of imap-login processes? Aug 10 16:23:34 director5 dovecot: imap-login: Disconnected: Connection queue full: I am seeing this logged several times per minute. From netstat, the server has 650 proxied IMAP connections and 450 proxied IMAPS connections, there are currently 10 imap-login proces...

Dear Folks, We had this assertion failure on Dovecot 2.1.1 on a particular mailbox. I upgraded to 2.1.5, and we have the same assertion failure, on the same mailbox: Apr 27 16:02:03 imap(nu-imaptest00258 at example.com): Panic: file mailbox-list-fs.c: line 156 (fs_list_get_path): assertion failed: (mailbox_list_is_valid_pattern(_list, name)) Apr 27 16:02:03 imap(nu-imaptest00258 at example.com): Error: Raw backtrace: /usr/lib64/dovecot/libdovecot.so.0(+0x4246a) [0x7f310c28146a] ->...

...d parser changed quite recently and that '#', '_' and '&' must not be escaped according to "Writing R extensions". However, the problematic .Rd portion seems (to me) OK \deqn{\theta(h) = \begin{cases} u_\beta \left(\mu + \frac{\nu(h)}{\Gamma(1 - \xi)} \right),&\xi \neq 0\\ \exp\left(\frac{\nu(h)}{\sigma}\right),&\xi = 0 \end{cases} }{u_beta (mu + nu(h) / Gamma(1 - \xi)), if \xi < 1, exp(nu(h)/sigma), otherwise} Removing the alignment tab character '&...

Hi, I'm not sure how use curve(dexGAUS(… None of the following four works: rt<- rexGAUS(100, mu=300, nu=100, sigma=35) m1<-gamlss(rt~1, family=exGAUS) curve(dexGAUS(rt=x, mu=300 ,sigma=35,nu=100), 100, 600, main = "The ex- GAUS density mu=300 ,sigma=35,nu=100") curve(dexGAUS(x=rt, mu=300 ,sigma=35,nu=100), 100, 600, main = "The ex- GAUS density mu=300 ,sigma=35,nu=100") cur...

...BiLinInit0= function(mCall,LHS,data) { model=as.character(mCall[[1]]) #I don't know what this does a <- tan(pi/4) # 45 degree angle (rising slope) b <- -tan(pi/4) # 45 degree angle (desending slope) c <- 0 # intercept nu <- 0.434 xy <- sortedXyData(mCall[["x"]], LHS, data) x <- xy[["x"]] x0 <- (min(x)+max(x))/2 + nu*((log(a, exp(1)) - log(b, exp(1))))/(a+b) #log defaults to base e value = c(nu,a,b,c,x0) names(value)=mCall[c("nu","a","b",&quot...

Hello, I am using SVM under e1071 package for nu-regression with 18 parameters. The variables are ordered factors, factors, date or numeric datatypes. I use the linear kernel. It gives the following error that I cannot solve. I tryed debug, browser and all that stuff, but no way. The error is: Error in get(ctr, mode = "function", en...

Hi, A. In a nutshell: The training error, obtained as "error (ret)", from the return value of a ksvm () call for a eps-svr model is (likely) being computed wrongly. "nu-svr" and "eps-bsvr" suffer from this as well. I am attaching three files: (1) ksvm.R from the the kernlab package,...

Can't anyone stop self-promotion and tell the poor guy what he needs. A T1/E1 X-over cable using an RJ-45 (8-cond.) is pinned out as follows: 1 - 4 2 - 5 3 - NU 4 - 1 5 - 2 6 - NU 7 - NU 8 - NU NU = Not Used I have not in my experience seen any problems with using a Good Quality Cat5 vs. Cat 3 (telco standard) cable for X-connects. YMMV, but you should be fine. As far as the shielding goes, I use UTP cables and Connectors all the time and some of my X-c...

...<- as.data.frame(cbind(runif(n=100),runif(n=100))) sample$Y <- ifelse(sample$V1>sample$V2,1,-1) ## Plot data; attach(sample) plot(V1, V2, type="p", pch=Y,main="Sample Data") ##Perceptron algorithm sample_m <- as.matrix(sample) w <- c(0,0); b <- 0; k <- 0; nu <- 1e-3 R <- max(sqrt(V1^2+V2^2)) repeat { for (i in 1:nrow(sample_m)){ if (sample_m[i,3]*(t(w)%*%sample_m[i,1:2] + b) <= 0) w <- w + nu*sample_m[i,3]*sample_m[i,1:2] b <- b +nu*sample_m[i,3]*R^2 k <- k+1 cat(k, w, b, "\n")...

...* pi * theta)), flag = FALSE) > [...] > > Similarly, the frailty.t() appears to use the joint negative log-density > of Student t random variables. HOWEVER, frailty.gamma() uses: > > frailty.gamma > [...] > list([...], penalty = -sum(coef) * > nu, flag = FALSE) > [...] > I would rather expect to see something like: > (1) penalty=sum(coef*nu-(nu-1)*log(coef)+lgamma(nu)-nu*log(nu)) > which is the joint negative log-density of gamma variables. Alternately, I > could also expect to see something like this: > (2) penalty=...