search for: nugget

dear list I am student M.S. statistics in department statistics . I am working in the function "nls" in the [R 2.3.1] with 246 data and want to fit the "exp" model to vectors( v and u ) but I have a problem to use it u 5.000000e-13 2.179057e+03 6.537171e+03 1.089529e+04 1.525340e+04 1.961151e+04 2.396963e+04 2.832774e+04 3.268586e+04 3.704397e+04 4.140209e+04

...Asym ~ Type * Treatment, lrc + c0 ~ 1), start = c(32.412, 0, 0, 0, -4.5603, 49.344), weights=varConstPower(fixed=list(const=0.1, power=1)), verbose=T) CO2.nlme.CAR<-update(CO2.nlme.var, corr=corCAR1()) CO2.nlme.gauss<-update(CO2.nlme.var, correlation=corGaus(form=~as.numeric(conc)|Plant,nugget=F), data=CO2) CO2.nlme.exp<-update(CO2.nlme.var, correlation=corExp(form=~as.numeric(conc)|Plant,nugget=F), data=CO2) But, if i change each of the 'conc' numbers slightly so that they are no longer identical between subjects i can only get the corCAR1 correlation to work while R cra...

...sults. Even though are not accurate, still might be a good indicator how much improvement there can be. A.Case. The classic: for 1:100 for (i in c(1:dimz)){ print(sprintf('Creating the %d map',i)); Shadowlist[,,i]<- GaussRF(x=x, y=y, model=model, grid=TRUE,param=c(mean,variance,nugget,scale,Whit.alpha)) } user system elapsed 1825.699 303.100 1063.352 -------------------------------------------------------------------------- B.Case. Same as above but with lapply instead of for Shadowlist<-lapply(1:dimz, function(i) { print(sprintf('Creating the %d ma...

...Measurements have been taken on 470 plots in March each year and georeferenced. Average density is supposed to be different each year. In a first approach, I would like to use a general model of this type, taking spatial correlation into account: gls(IKAfox~an,correlation=corExp(2071,form=~x+y,nugget=1.22),data=renliev) but I get > gls(IKAfox~an,correlation=corExp(2071,form=~x+y,nugget=1.22),data=renliev) Error in getCovariate.corSpatial(object, data = data) : Cannot have zero distances in "corSpatial" I understand that the 470 geographical coordinates are repeated thr...

...estimates of fish biomass (lbloat) in 1,000 meter long intervals (elementary sampling units) from multiple transects (each 20-30 km long, tranf) in two different lakes and three different years. bloat.gamm1 <- gamm(lbloat ~ s(depth), correlation=corSpher(c(30000, 0.01),form = ~ x+y|tranf, nugget=TRUE), data=fish3) bloat.gamm2 <- gamm(lbloat ~ lakef + s(depth), correlation=corSpher(c(30000, 0.01),form = ~ x+y|tranf, nugget=TRUE), data=fish3) bloat.gamm3 <- gamm(lbloat ~ lakef + s(depth, by=lakef), correlation=corSpher(c(30000, 0.01),form = ~ x+y|tranf, nugget=TRUE), data=fish3...

...lp, Can anyone explain the meaning of the warning, Singular precision matrix in level -1, block 1 ? Or how to track down where it comes from? More precisely, using the nlme package, I'm issued with the warning itt2 <- lme(lrna~rx.nrti+lbrna, random=~1|patid, cor=corExp(form=~days|patid,nugget=T), weights=varPower( form=~lbrna),data=rna3) Warning messages: 1: Singular precision matrix in level -1, block 1 2: Singular precision matrix in level -1, block 1 the output is: Linear mixed-effects model fit by REML Data: rna3 Log-restricted-likelihood: -4990.142 Fixed: lrna ~ rx.nrti + l...

...thenticate ok, but then I get this error. The Finder cannot complete the operation because some data in "smb:// <theserver/theshare>" could not be read or written. (Error code -36). I check the logs and found these two entries every time I've tried to connect. May 3 12:43:49 Nugget kernel[0]: smb_maperr32: no direct map for 32 bit server error (0xc00000bd) May 3 12:43:49 Nugget kernel[0]: smb_maperror: Unmapped DOS error 1:52 I also did a tcpdump while trying to connect and can see that successful connections are occurring between me and the server on 137/ UDP, 139/TCP, and...

...effects), and differ only in that the second model includes a corExp structure for spatial autocorrelation. (more detailed explanation of the models at end). Specific questions: 1. The second model estimates two additional parameters in the process of fitting the corSpatial object - the range and nugget of the spatial autocorrelation. Based on this, I would expect the second model to have two fewer residual degrees of freedom. However, the summary function reports that both models have the same number of residual degrees of freedom. Why is this? (Interestingly, the difference in AIC between the t...

Dear Members, I am trying to fit a variogram model using fit.variogram function from the gstat package. The figure showing my experimental variogram can be seen here: https://i.stack.imgur.com/UZXw4.png My code line for this operation is: > c2.vgm.fit<-fit.variogram(c2.vgm.exp,vgm(nugget=0, psill=400,model="Exp",range =40000),fit.method = 7) The system throws following error message: Error in switch(model, exponential = fit.exponential(v.object, c0 = nugget, : EXPR must be a length 1 vector Any help to know the cause of the error and avoiding it would be highly app...

I was at VON in Boston today and saw on the Digium stand a Cisco 7960 with a picture of Tux and the Asterisk log on its display. I WANT IT!!!!! Anyone know where I can download this file please? -- Mark, G7LTT/KC2ENI Randolph, NJ http://www.g7ltt.com

...le. I thought it would be a simple step to loop through the file and insert the values in to the db. In addition to the generated admin_controller.rb file containing list, update etc i put def csv_preview require ''faster_csv'' @csv_trans_action = FasterCSV.read("c:/nugget/statement.csv") end def csv_load require ''faster_csv'' @csv_trans_actions = FasterCSV.read("c:/nugget/statement.csv") for trans_action in @trans_actions for csv_trans_action in @csv_trans_actions trans_action.trans_date_at = csv_tran...

Hi, I'm trying to introduce a (spatial) exponential correlation structure (with range=200 and nugget.effet of 0.3) in a lme model of this form: lme(ARBUS~YEAR, random=~1|IDSOUS). The structure of the data is "IDSOUS" "XMIN" "YMAX" "YEAR" "ARBUS" with 2 years of data and 5600 points for each year. I do: corstr<-corExp(value=200, form=~...

...VM I wanted to have started automatically. Since updating to 4.6, this doesn't work.? Period. I'm having a hard time finding what needs to change in order to get VMs to auto start on reboot of a server.? Can someone please point me to documentation that I have missed or share the gold nugget? -- Craig Thompson, President Caldwell Global Communications, Inc. -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.centos.org/pipermail/centos-virt/attachments/20160815/97459d4c/attachment-0002.html>

...;> library(RandomFields) >> library(MASS) >> >> x <- c(0, 200.0, 0.025) >> y <- c(0, 1450.0, 0.025) >> >> model <- "stable" >> scale <- 0.1 >> alpha <- 2 >> mean <- 0 >> variance <- 1.0 >> nugget <- 0 >> >> f <- GaussRF(x=x, y=y, model=model, grid=TRUE, gridtriple=TRUE, > + param=c(mean, variance, nugget, scale,alpha)) > Error: cannot allocate vector of size 3625515 Kb > In addition: Warning messages: > 1: Reached total allocation of 2560Mb: see h...

...n the hope that one of you folks might be able to correct me. Thanks in advance for your time. The issue is in 2 parts. (A) I've used gls to fit a model with two fixed effects and a corExp object. By my count, this fitting process estimates 5 parameters: (Intercept), l10area, newx, range, and nugget. With 118 total df, there should be 118 - 5 = 113 residual df. However, the output from summary.gls reports 115 residual degrees of freedom. Is this an error in summary or gls, or is there an error in my count? (B) Summary.gls reports logLik=-273.6. Using my count of 5 estimated parameters, the AI...

...orm=~STX1+STY1|tran)) Variogram(test2, form=~STX2+STY2|tran) # Varigoram has distance of 0! plot(Variogram(test2, form=~STX2+STY2|tran)) mod1<-lme(log(yabund)~factor(log)+STX1+STY1, random=~1|tran, data=dat, corr=corSpatial(form=~STY1+STX1|tran, type="exponential", nugget=T), method="REML") mod2<-lme(log(yabund)~factor(log)+STX2+STY2, random=~1|tran, data=dat, corr=corSpatial(form=~STY2+STX2, type="exponential", nugget=T), method="REML") # Error in getCovariate.corSpatial(object, data = data) : # Cannot have zero di...

I've got soft phone that allows me to dial SIP URI's. I'd like to route these calls through a provider to be completed, because I'm beind a NAT box and doing it directly doesn't work. Right now I've got an extension defined like this: Dial(IAX2/${FWDUSERID}:${FWDPASSWD}@${FWDSERVER}/**356<username>) This will connect a call to FWD and call a user at FWD. It works

...lack of knowledge pertaining to random fields in general and not the RandomFields package specifically, in which case someone can set me straight. I like the versatility of the extended definition form that can be used in RandomFields. For example, specifying an isotropic exponential model with a nugget of 1 as: modelA=list("+", list("$", var=2, scale=2, list("exp")), list("$", var=1, list("nugget")) ) Ultimately I would like to work with anisotropic models (and data), and...

...an NFS mount is showing 64M in size, but when copying the file to a local drive, it shows 2.5MB in size. My NFS server is hardware Raided with a volume stripe size of 128K were the volume size is 20TB. My NFS clients are the same distro as the server being Centos. Is this due to my stripe size? Nuggets are appreciated. - aurf

[hoping to redeem myself for my last spurious bug report] From ?corRatio: Letting d denote the range and n denote the nugget effect, the correlation between two observations a distance r apart is (r/d)^2/(1+(r/d)^2) when no nugget effect is present and (1-n)*(r/d)^2/(1+(r/d)^2) when a nugget effect is assumed. This disagrees with the C code (corStruct.c) /* Rational class */ static double ratio_corr(...