search for: ngrid

...llowing code has been cut and pasted directly from the web resource - the only change I have made is to fill in my values for T, n, low/high se/sp and the alpha beta for the distributions) prevalence.bayes<-function(theta,T,n,lowse=0.5,highse=1.0, sea=1,seb=1,lowsp=0.5,highsp=1.0,spa=1,spb=1,ngrid=20,coverage=0.95) { ibeta<-function(x,a,b) { pbeta(x,a,b)*beta(a,b) } ntheta<-length(theta) bin.width<-(theta[ntheta]-theta[1])/(ntheta-1) theta<-theta[1]+bin.width*(0:(ntheta-1)) integrand<-array(0,c(ntheta,ngrid,ngrid)) h1<-(highse-lowse)/ngrid h2...

Hi, I am doing kernel density plots, and am trying to make the lines thicker. I comparing three groups, in sm.density.compare. I tried changing lwd to make the line sthicker right on the density compare call, but was not able to do it. There is not an option in sm.options to specify line thickness, as well as cex.ylab or cex.xlab- I tried it and it does not change the thickness of the lines.

...s(t(res))), log(abs(t(synthesised))))) Error in loessFit(log(abs(t(res))), log(abs(t(synthesised)))) : ? binary operation on non-conformable arrays Timing stopped at: 0.57 0.06 0.64 NA NA > system.time(fit <- plotPrintTipLoess( log(abs(t(res))), log(abs(t(synthesised))))) Error in 1:layout$ngrid.c : NA/NaN argument In addition: Warning messages: 1: In object$printer : ? $ operator is invalid for atomic vectors, returning NULL 2: In object$M : $ operator is invalid for atomic vectors, returning NULL 3: In object$A : $ operator is invalid for atomic vectors, returning NULL 4: In object$M : $...

Dear all, Does anybody know how to get more evaluation points in performing Nonparametric analysis of repeated measurements data with "sm" library. The following command gives the estimation for 50 points, by I would like to increase to 100 points But I do not know how to do that. library(sm) provide.data(citrate, options=list(describe=FALSE)) provide.data(dogs,

...ion=na.omit) All is fine untill I bacame... MCMC scan 1000 of 5000 (CPU time: 32.125 s) MCMC scan 2000 of 5000 (CPU time: 61.797 s) MCMC scan 3000 of 5000 (CPU time: 91.281 s) MCMC scan 4000 of 5000 (CPU time: 120.750 s) MCMC scan 5000 of 5000 (CPU time: 150.047 s) Error in matrix(foo$f, nrow = ngrid, ncol = ngrid) : Versuch ein Attribut von NULL zu setzen.... Some idea? Best regards, rodrigo -- Besuchen Sie meine Web Seite http://rodrigoherrera.blogspot.com/ Viel Spass beim Surfen ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ Habe nun, ach! Philosophie, Jurist...

...~STR, nest=T, data=data) model <- svyglm(fmla, design=design) > #works fine as evidenced by summary (not shown) termplot(model, > data=model.frame(design), partial.resid=T, se=T, > smooth=make.panel.svysmooth(design)) Waiting to confirm page change... Error in rowsum.default(c(rep(0, ngrid), w), c(1:ngrid, findInterval(mm[, : incorrect length for 'group' > make.panel.svysmooth(design) function (x, y, span = 2/3, col.smooth = "red", col = par("col"), bg = NA, pch = par("pch"), cex = 1, ...) { if (is.null(bandwidth)) bandwi...

...ifferent values. My second question is regarding the sm library: d<-sm.regression(lprice, XX[,i], h=sd(lprice), display="none") will only generate 50 point estimates while NROW(XX[,i]) = 3897 and when I do d<-sm.regression(lprice, XX[,i], h=sd(lprice), display="none", ngrid=NROW(lprice)) I get the right dimension of the estimated points but again they are not estimated at the points in lprice. Any help is greatly appreciated. P.S. I have Bowman and Azzalini book but unfortunately it does not clarify the procedures in sm.regression() Jean,

Hi Everyone, I have been trying to read some Arrayvision files( 2 channel cDNA) and am having some problem. My code is : setwd('C:/work/data/limma/ndd1'); files <- c('ndd1_1.txt','ndd1_2.txt','ndd1_3.txt'); RG=read.maimages(files,"arrayvision",sep="\t"); #Normalisation MA=normalizeWithinArrays(RG); #plotPrintTipLoess(MA); #Fit Linear

Hi Samba Team, We are currently running on Samba 2.2.3a version. As part of the technology upgrade, management is willing to upgrade the current samba versions of in production machines. Can you please tell me how long I can use the 2.2.3a version? And if you suggest me to upgrade to a higher version, which version will best suite for us ?. Thanks in Advance.....Hisham Unless expressly stated

...library: > > > > d<-sm.regression(lprice, XX[,i], h=sd(lprice), display="none") > > will only generate 50 point estimates while NROW(XX[,i]) = 3897 > > and when I do > > d<-sm.regression(lprice, XX[,i], h=sd(lprice), display="none", > > ngrid=NROW(lprice)) > > > > I get the right dimension of the estimated points but again they are not > > estimated at the points in lprice. > > > > > > Any help is greatly appreciated. > > > > P.S. I have Bowman and Azzalini book but unfortunately it does n...