search for: fitt

..."Warning message: 'newdata' had 1 rows but variable(s) found have 24 rows " here is the code i am using for locfit......... xt<-read.table("x.dat",header=FALSE) yt<-read.table("y.dat",header=FALSE) tst<-read.table("test.dat",header=FALSE) fitt<-locfit.raw(sapply(xt,"[",1:24),sapply(yt,"[",1:24), deg=1,maxk=700) fr<-fitted(fitt) result<-fr hltt<-predict(fitt,newdata=tst) hltt<-predict.locfit(fitt,newdata=sapply(tst,"[",1:1)) res_tt<-hltt cat(file="prd_model.dat",result, sep=&q...

Hello, I am analyzing some climate time series data using the Mann Kendall package and was wondering if there was a way to calculate the trend using Sen's nonparametric estimator slope in R? Thank you in advance, Barry _________________________ Barry Baker, Ph.D. Global Climate Change Initiative The Nature Conservancy 2424 Spruce St., Suite 100 Boulder, CO 80302 Tel: (303)-541-0322 Fax:

I am fittting a linner regression with constrained parameters, saying, all parameters are non-negative and sum up to 1. I have searched historical R-help and found that this can be done by solve.QP from the quadprog package. I need to assess the significance of the coefficient estimates, but there is no sta...

...ount of transformed inputdata to find the best fit and then step backwards down to e.g. rank 5 in order to get a smooth curve, i could start right from the beginning with only very few but therefore very flexible funktions (litle similar to box-plot) which adopt themselve automaticaly to an optimal fitt and r-squered. can anyone tell me if such a flexible adoptive function also exists in r+? (i could not find anything on the r+ pages with these search words). waiting desperately for help joerg

Hello I?m need fitt growth curve with data length-age. I want to evaluate which is the function that best predicts my data, to do so I compare the Akaikes of different models. I'm now need to evaluate if changing the initial values changes the parameters and which do not allow to estimate the model. To do this I u...

Dear colleagues, ?arima? returns directly the 1-step ahead errors but I am interested in obtaining other h-step ahead errors for several ARIMA models I have fitted. Is there any way I can obtain this with R? Any help would be appreciated. Sincerely, Nuno Prista _________________________ CO - FCUL, Lisboa, Portugal CQFE - ODU, Norfolk, USA

...lt;- step(model) out <- data.frame(FID=covar$FID, IID=covar$IID, st$model) # formatting for plink write.table(out, outfiles[i], col.names=TRUE, row.names=FALSE, quote=FALSE, sep="\t") } and i end up having 500,000 files for each SNP 1- Is this way is approperiate in assessing the fittness of logistic model ? 2- how can i now use all of these file in a way that i dont need to write the glm formula for each file ? 3- Is there is a better way to save the results so i will not have to deal with 500,000 files ? Please need your help Alajmi [[alternative HTML version deleted]]

...forums, but I wanted post a mail about it here to be sure that it gets noticed by you compiz developers. The problem is that compiz seems to make the screen edges insensitive to clicks, meaning that e.g. a panel item isn't clickable when the pointer is directly at the screen edge. This breaks Fitt's Law and should therefore be fixed for usability reasons. Here is the forum thread about this issue, in which somebody seems to have already tracked down the cause of the problem: http://forum.compiz-fusion.org/showthread.php?t=6885 I also reported a bug about it on bugzilla.freedesktop.org:...

...8 EDT .Commissioners/ drwx------ 5 vmail vmail 512B 2014-07-08 13:41:19 EDT .Commissioners.Cordrey/ drwx------ 5 vmail vmail 512B 2014-07-03 17:01:31 EDT .Dover/ drwx------ 5 vmail vmail 512B 2014-07-19 08:32:05 EDT .Drafts/ drwx------ 5 vmail vmail 512B 2014-07-18 19:07:15 EDT .Fitts/ drwx------ 5 vmail vmail 512B 2014-07-08 13:41:19 EDT .Gering/ drwx------ 5 vmail vmail 512B 2014-07-19 08:32:05 EDT .Green/ drwx------ 5 vmail vmail 512B 2014-07-18 18:59:06 EDT .Haywood/ drwx------ 5 vmail vmail 512B 2014-07-03 17:01:31 EDT .Hornik/ drwx------ 5 vmail vma...

...Accept Sess-Refresh : uas Sess-Expires : 1800 secs Min-Sess : 90 secs RTP Engine : asterisk Parkinglot : parkinglot_test Use Reason : No Encryption : No File "debug.txt" shows that when a call is parked, it is NOT sending over to parkinglot_test, but parkinglot_fitts. Debug.txt CODE: SELECT ALL == Using UDPTL CoS mark 5 == Using SIP RTP CoS mark 5 -- Called 101-test == Extension Changed 101[local-extensions-test] new state Ringing for Notify User 102-test -- SIP/101-test-00000035 is ringing -- SIP/101-test-00000035 is ringing == Extension...

...an),di = S-1, mi = S-1, do = S-1, fr = S-1, b1 = S-1, b2 = S-1, b3= S-1, b4= S-1, dum = S-1, alpha0 = S^2, alpha = S, beta = S) upperB = c(a = 10*abs(Mean), di = 1-S, mi = 1-S, do = 1-S, fr =1-S, b1 = 1-S, b2 = 1-S,b3 = 1-S, b4 = 1-S, dum = 1-S, alpha0 = 100*Var, alpha = 1-S, beta = 1-S) fitt<-maxLik(start=param, logLik=garch2,method="BHHH", x=dat2$r_csi,Di=dat2$Di,Mi=dat2$Mi,Do=dat2$Do,Fr=dat2$Fr,y=dat2$r_t,z=dat2$r_sp,d=dat2$f) Note that optim always breaks down: nlminb and the BFGS and BHHH algorithmus from the maxLik-package work fine. The estimated coefficients are...

?tryCatch -- Sent from my phone. Please excuse my brevity. On October 20, 2017 7:37:12 AM PDT, Evangelina Viotto <evangelinaviotto at gmail.com> wrote: >Hello I?m need fitt growth curve with data length-age. I want to >evaluate >which is the function that best predicts my data, to do so I compare >the >Akaikes of different models. I'm now need to evaluate if changing the >initial values changes the parameters and which do not allow to >estimate &...