similar to: curve fit

how to realize the following SQL command in R? select distinct A, B, count(C) from TABLE group by A, B ; quit; Best Regards --------------------------------- Sponsored Link Get a free Motorola Razr! Today Only! Choose Cingular, Sprint, Verizon, Alltel, or T-Mobile. [[alternative HTML version deleted]]

Is there anyway to change any y[i] value (i=2,...6) to make following NLS workable? x <- c(0,5,10,15,20,25,30) y <- c(1.00000,0.82000,0.68000,0.64000,0.66667,0.68667,0.64000) lm(1/y ~~ x) nls(1/y ~~ a+b*x^c, start=list(a=1.16122,b=0.01565,c=1), trace=TRUE) #0.0920573 : 1.16122 0.01565 1.00000 #Error in numericDeriv(form[[3]], names(ind), env) : # Missing value or

x<-c(-1.873....,-0.121) # 23 numerics; x.arma12 <- armaFit(x ~ arma(1,2)) #estimates y[t]= -0.11465 - 0.23767 y[t-1] - 0.14230 e[t-1] -0.85770 e[t-2] + e[t]; # ? how to predict 46 steps ahead based on 23 data points? # the following doesn't work since n is in armaSim rather than armaFit; predict(x.arima12, n.ahead=46) # Thanks ---------------------------------

I typed the following commands but it still use maxiter=50 after the 2nd command: nls.control(maxiter = 1000) nls(......) Thanks! --------------------------------- [[alternative HTML version deleted]]

> y [1] 1 11 42 64 108 173 214 > t [1] 1 2 3 4 5 6 7 > nls(1/y ~ c*exp(-a*b*t)+1/b, start=list(a=0.001,b=250,c=5), trace=TRUE) 29.93322 : 0.001 250.000 5.000 Error in numericDeriv(form[[3]], names(ind), env) : Missing value or an infinity produced when evaluating the model # the start value for b is almost close to final estimates, # a is usually

I simulate a uniform data and then transformed into Burr(1,3,1) data, which is of pdf: f(x)=[3*(x^2)] / [(1+x^3)^2], x>0 How can I perform a goodness-of-fit test (k-s, anderson-darling,chisq,cramer-von mises,...) on it (should highly accept) to get test-statistics & p-values? Thanks! Sincerely, Shelton Jin -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-

Hi, May I ask how to download and install R software freely from internet? Thanks a lot! Best Regards, Shelton -------------------------------------------------- Industrial Statistician - Quality Engineering Alcatel Fiber Plant 2512 Penny Road Claremont, NC 28610-0039 Tel: 828-459-3092 Fax: 828-459-8840 xiaodong.jin at cable.alcatel.com

I'm using R 2.1.1 and just successfully installed packages tseries, fseries. I try to run example http://www.maths.lth.se/help/R/.R/library/tseries/html/garch.html But it shows > x.arch <- garch(x, order = c(0,2)) # Fit ARCH(2) Error: couldn't find function "garch" Then I run command > help.search("garch") it shows the R information.

May I ask how to include following procedures into R script (such as first.r) such that it will do an automatical call next time after I open R? Packages -> Load Packages -> splines Thanks, Shelton --------------------------------- [[alternative HTML version deleted]]

I have a problem moving from multistratum aov analysis to lmer. My dataset has observations of ampl at 4 levels of gapf and 2 levels of bl on 6 subjects levels VP, with 2 replicates wg each, and is balanced. Here is the summary of this set with aov: >> summary(aov(ampl~gapf*bl+Error(VP/(bl*gapf)),hframe2)) > >Error: VP > Df Sum Sq Mean Sq F value Pr(>F) >Residuals

Dear Listers, I am drawing a plot of two density curves, for male and female incomes. I would like to shade/hatch/color (whatever) the areas under the curves which are distinctive for each gender. This is the code I have tried so far: m <- density(topmal.d$y, bw = "sj") f <- density(topfem.d$y, bw = "sj") par(mfrow = c(1,1)) plot(x = c(0,400), y = c(0,0.02), type =

Always keep the mailing list in cc. The code runs for each row in the data. However I get the feeling that there is a mismatch between what you think that is in the data and the actual data. ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance Kliniekstraat 25 1070

Hi, I am wondering if there is a simple way to fix the problem I am having. For unknown reason, I could not get the full name of the factors to be printed in the summary. I have tried to used summary.lm as well but the problem still persists. SJ$Weekday <-

I have two series of events over time and I want to construct a graph of the relative frequency/density of these events that allows their distributions to be sensibly compared. The events are the milestones items in my project on milestones in the history of data visualization [1], and I want to compare trends in Europe vs. North America. I decided to use a graph of two overlaid density

Hi all, I'm trying to resolve a weird issue with SIP routing. I have a number of SIP trunks, from a selection of providers, all of which are registered in sip.conf: [general] context=default allowguest=no allowoverlap=no udpbindaddr=0.0.0.0 tcpenable=yes tcpbindaddr=0.0.0.0 transport=udp bindport=15060 srvlookup=yes allowsubscribe=yes

On 26 September 2010 18:56, Nathan Jeffords <blunted2night at gmail.com> wrote: > The syntax for the invoke instruction is a little misleading. %x is a value > that is being generated by the instruction, not passed to is. It is no > different in that regard as to say '%x = call @eh.exception ...'. Since you > don't specify the type in that type of assignment, I chose

I did a CentOS 4.4 fresh server install. Using it to run Squid. Named is also running but I get this error. [root at sj ~]# /etc/init.d/named status rndc: connection to remote host closed This may indicate that the remote server is using an older version of the command protocol, this host is not authorized to connect, or the key is invalid. [root at sj ~]# /etc/init.d/named restart Stopping

For a real data column X value ranged between (56.4521,32317.9) with missing values, I need to overlay 2 plots: histogram & fitted gamma density. I use following to generate histogram. xbk_seq(50,33000,by=100) hist(x,breaks=xbk) But I don't know how to get "fitted gamma density"? In SAS proc capability, I got Shape=2.59, Scale=3481). But when I do plot(dgamma(x,

How to make it correct? [1] Error in parse(file, n = -1, NULL, "?") : syntax error at 1: {\ #source('c:/Documents and Settings/xxxxxxx/My Documents/Archives/R/test.r') #source("c:\\Documents and Settings\\xxxxxxx\\My Documents\\Archives\\R\\test.r") [2] Error in file(file, "r", encoding = encoding) : unable to open connection In

Dear R-ians, I'm looking for a computational simplified formula to calculate a measure for heterogeneity (let's say H ): H = sqrt [ (Si (Sj (Xi - Xj)?? ) ) /n ] where: sqrt = square root Si = summation over i (= 0 to n) Sj = summation over j (= 0 to n) Xi = element of X with index i Xj = element of X with index j I can simplify the formula to: H = sqrt [ ( 2 * n * Si (Xi) - 2 Si (Sj