search for: nile

1. The function "KalmanLike" seems to change its inputs AND PREVIOUSLY MADE copies of the inputs. Consider the following (using R 2.1.0 patched under Windows XP): > Fig2.1 <- StructTS(x=Nile, type="level") > unlist(Fig2.1$model0[2:3]) a P 1120 286379470 > tst2 <- tst <- Fig2.1$model0 > tst23 <- tst[2:3] > tst23u <- unlist(tst23) > nile.KL <- KalmanLike(nile, tst2) > unlist(tst[2:3]) a P 798.3682...

When doing ?Nile, the url for the data source is dead. It says http://www.= ssfpack.com/dkbook/ but this has changed to=20 http://www.ssfpack.com/DKbook.html Version: platform =3D i386-redhat-linux-gnu arch =3D i386 os =3D linux-gnu system =3D i386, linux-gnu status =3D major =3D 2 minor =3D 10.0 year =3D...

Dear all, I'm trying to understand how the strucchange package is working and I have been looking at the examples given for the Fstats() function. The first example (Nile), shows one peak in the F-stats and one breakpoint is estimated, that can be plotted using the following code ## Nile data with one breakpoint: the annual flows drop in 1898 ## because the first Ashwan dam was built data("Nile") plot(Nile) ## test the null hypothesis that the annual flow...

Is it possible to use some form of robust regression with the breakpoints routine so that it is less sensitive to outliers? --Rich Richard Kittler Advanced Micro Devices, Inc. Sunnyvale, CA

People, It appears that there is no way of getting Boxplots to plot using Mean, SD, Max & Min - is there something else that would do what I want? I couldn't find it . . Thanks, Phil. -- Philip Rhoades GPO Box 3411 Sydney NSW 2001 Australia E-mail: phil at pricom.com.au

UNSUBSCRIBE AT THE BOTTOM ------------------------------------------------------- Dear Subscriber/Member, You don't need to spend hundreds of dollars on Norton SystemWorks 2003! You can get the exact same thing from us! Norton SystemWorks 2003 is the most technologically advanced method of AntiVirus and Computer Stability set of programs ever to be offered at a public level. Get it today

Hi R community, I'm running R 2.7.2 on Windows XP SP2. I'm trying to (1) plot loess lines for each of my groupings using the same color for each group; (2) plot loess predicted values. The first part is easy: data1 <- data.frame(Names=c(rep("Jon",9),rep("Karl",9)),Measurements=c(2,4,16,25,36,49,64,81,100,1,2,5,12,17,21,45,54,67),PlotAt=c(1:9,1:9)) data2 <-

Dear all, I'd like to discuss about a possible bug in function StructTS of stats package. It seems that the function returns wrong value of the log-likelihood, as the added constant to the relevant part of the log-likelihood is misspecified. Here is an simple example: > data(Nile) > fit <- StructTS(Nile, type = "level") > fit$loglik [1] -367.5194 When computing the log-likelihood with other packages such as KFAS and FKF, the loglikelihood value is around -645. For the local level model, the likelihood is defined by -0.5*n*log(2*pi) - 0.5*sum(log(F_t) +...

Dear all, I am doing some time series analysis with R now. The problem is, when I create a time series object using "ts" function, then after fitting the model, the predicted values cannot be plotted with "ts" object together using "lines". ie. ts.series<-ts(x,start, end) plot(ts.series) lines(predict(fit.model)) This doesn't work. Does anyone know about

Here I am expecting to be configured following scenario: User calls : it will play a sound will ask for input DTMF, then call will be given to particular extension for any DTMF entered. But its not working as expected. I have attached the dial plan file. -------------- next part -------------- An HTML attachment was scrubbed... URL:

I am attempting to get fax detection to work. I am using a NETjet-s card under ISDN4Linux. Asterisk does not seem to be detecting the fax tone. I have tried following as a test: [MainMenu] exten => s,1,Answer exten => s,2,DigitTimeout(3) exten => s,3,ResponseTimeout(5) exten => s,4,Background(Welcome) exten => s,5,Background(MainMenu) exten => fax,1,Dial(Zap/1,,d) [FaxTest]

Dear R Users, I have 12 data frames, each of 12 rows and 2 columns. e.g. FeketeJAN MEAN SUM_ AMAZON 144.4997874 68348.4 NILE 5.4701955 1394.9 CONGO 71.3670036 21196.0 MISSISSIPPI 18.9273250 6511.0 AMUR 1.8426874 466.2 PARANA 58.3835497 13486.6 YENISEI 1.4668313 592.6 OB 1.4239179 559.6 LENA 0.9342164 387.7 NIGER 4.7245709 826.8 ZAMBEZI 76....

> niles@atheos.net wrote: > > Christopher, > > > > Just as a total guess, check to be sure the PoE portion > >of his cat5 cable are not either grounded, or touching > > each other. Also, be sure he's not connected to a PoE > >capeable switch since cisco is reverse...

I am trying to finalize the use of rsync for updatiung a new nfs server before we take the old one offline. I keep getting the following warning during the rsync process: Warning: unexpected rad size of 0 in map_ptr Any ideas where this comes from and how to make it go away? I am using rsync 2.5.2 on Solaris 8 to pull data from rsync 2.5.2 on Solaris 7. Bob roconnor@vectorpartners.com

I would like to get the result of acf min of lag 2 and max of lag 50. When I use time series ( acf, lag.max = 50, type="covariance"), I got lag 0 to lag 50. How do I get lag 2 to lag 50? Sincerely, Stephen

Does anyone know where I can get R code for plotting the Brown , Durbin and Evans cumsum procedure ( 1975 ) ? I wrote my own code but I am a little worried that my confiodence bands may not be correct ( I find the formula in the original paper confusing and S+Finmetrics has a formula but that formula implies that there should be 4 lines as far as I can tell ) so I would like to see someone

...g to fix for someone with a clue. Unfortuately, I lack clue wrt to wine development. The program is available from http://www.aopa.org/flight_planner/ However, you need an AOPA account to login. I'm happy to run any tests for anyone that's interested in exploring this issue. Thanks, Rick Niles.

I have a question, I am unclear on how to park a call. I know that you are supposed to be able to press "#" and then transfer the call to extension 700. However, * doesn't seem to be graping the dtmf. I am using dtmfmode=inband. Asterisk is in the media path as well. Thanks in advance Justin

Hi all, Using Beamer, in order to highlight a piece of R code I do something like this - note the "\structure" and "\alert" commands: \begin{semiverbatim} > mleOut <- \structure{dlmMLE}(Nile, + parm = c(0.2, 120), # initial values for optimizer + lower = c(1e-7, 0)) \alert<2>{# V must be positive} > mleOut$convergence \alert<3>{# always check this!!!} [1] 0 \end{semiverbatim} How can I get a similar effect using Sweave? Thank you in a...

Dear all, I have a matrix called combine86 which looks as follows: > combine86 Sim Mean Obs Mean Sim Sum Obs Sum AMAZON 1172.0424 1394.44604 553204 659573 NILE 262.4440 164.23921 67973 41881 CONGO 682.8007 722.63971 205523 214624 MISSISSIPPI 363.0758 142.59883 124535 49054 AMUR 143.5857 89.30434 36040 22594 PARANA 702.3793 388.03030 162952 89635 YENISEI 208.1396 174.52722 83464 70509 OB...