search for: periodicity

Hi, I'm new here, I'd like to map a set of values in a data frame to another set of values in another frame (think a c++ map<>, a python dict or a hash lookup) Specifically, I have a data frame of information about trading periods in the NZ electricity market. For each period I have information about each bus (or node) in the network. I'd like add a "column" to the

Thankyou very much Marc for that nifty little script. When I use it on my real dataset though, the lines are fat in the middle and thinner towards the ends. I guess it's because "lines" draw one fitted line for each x, and if you have hundreds of x, this turns into a line that is thicker that it should be (due to rounding errors perhaps). I got the tip to use

Dear R Gurus, Thanks for any help in advance! Date.frame: Returns.names X id ticker date_ adjClose totret RankStk 258060 258060 13645T10 CP 2001-06-29 18.125 1877.758 My data frame is in the above format. I would like to filter by period, per id (every 125 days) each consisting of 250 days, I.e. 1-250, 126-375, etc. One important thing to note is that not all

...aks is highly variable. It's something like this: patient <- data.frame( day = c(1, 3, 5, 8, 10, 12, 15, 17, 19, 22, 24, 26), protein = c(5, 3, 10, 7, 2, 8, 25, 12, 7, 20, 10, 5) ) plot(patient$day, patient$protein, type="b") My goal is two-fold: firstly, I need to test for periodicity, and secondly, I need to try to predict the temporal location of future peaks. Of course, the peaks might be occurring on unmeasured days. I have been looking at this model: wave.form <- deriv3( ~ sin(2*pi*((day-offset)/period + 0.25)) * amplitude + mean, c("period", &qu...

Hello all, I am trying to understand how to specify nested factors when using coxph(), and if it is appropriate to nest these factors in my situation. In the simplest form, I am testing two different temperatures, with each temperature being performed twice in different experimental periods (e.g. Temp5 performed in Period A and C, Temp4 performed in Period B and D) I am trying to see if survival

...5<-c(0,0,0,0,0,10,0,0,0,0,10,0,0,0,0,0,0,10,0,0,0,0,10,0,0,0,0,0,0,10 ,0,0,0,0,10,0,0,0,0,0,0,10,0,0,0,0,10,0,0,0,0,0,0,10,0,0,0,0,10,0) par(mfrow=c(2,1)) spec.pgram(period.6, log="no") spec.pgram(period.5, log="no") The first series has period 6 and shows its periodicity in 1/6 and harmonics. In the second series I was expecting to see a signal occurring at frequency 1/5 but it does not. Can anybody explain me why? Thanks in advance, Nuno ______________________________________________ Centro de Oceanografia - IO-FCUL, Portugal Center for Quantitativ...

Dear R Users, I am working with gsub for the first time. I am trying to remove some characters from a string. I have hit the problem where the period is the shorthand for 'everything' in the R language when what I want to remove is the actual periods. In the example below, I simply want to remove the periods as I have removed the comma, but instead the complete string is wiped out. I

I am looking for help to analyze an unbalanced AB/BA cross-over design by requesting the type III SS ! # Example 3.1 from S. Senn (1993). Cross-over Trials in Clinical Research outcome<-c(310,310,370,410,250,380,330,270,260,300,390,210,350,365,370,310,380,290,260,90,385,400,410,320,340,220) subject<-as.factor(c(1,4,6,7,10,11,14,1,4,6,7,10,11,14,2,3,5,9,12,13,2,3,5,9,12,13))

Hi dear R community, I have up to 12 measures of a protein for each of 6 patients, taken every two or three days. The pattern of the protein looks periodic, but the height of the peaks is highly variable. I'm testing for periodicity using a Monte Carlo simulation envelope approach applied to a cumulative periodogram. Now I want to predict the location of the peaks in time. Of course, the peaks might be occurring on unmeasured days. Sadly, an NDA prohibits me from sharing the real data. The data look something like this: ##...

Hi its me again I don't mean to get on your nerves, but the use of R proofs to be a bit more complicated than envisaged. I would like to calculate the mean of a group of values, here "wage_accepted". The group is determined by the stage and period, so in the end there should be a column with the values of the wages in period 1 stage1, period 1 stage 2, period 2 stage1...

-----BEGIN PGP SIGNED MESSAGE----- ============================================================================= FreeBSD-SA-01:12 Security Advisory FreeBSD, Inc. Topic: periodic uses insecure temporary files Category: core Module: periodic Announced: 2001-01-29

I'm wondering if there's a way to set which periodic-announce file is played from my dialplan, much like setting the monitor-filename. Something like this: exten => s,n, Set(PERIODIC_ANNOUNCE=foo) This would be a great feature if it doesn't already exist. Or perhaps there's a better way to do this. Thanks for your time. -- Regards, Robert Broyles

R Help Community I'm trying to understand time series (TS) objects. Thought I understood but recently have run into a series of error messages that I'm not sure how to handle. I have 15 years of quarterly data and I typically create a TS object via something like... data.ts <- ts(mydata, start = 2002, frequency = 4) this create a matric as opposed to a vector object as I receive a

* OS name and version, Gentoo, updated * exact NUT version, 2.6.2 (also tried 2.6.0-r1) * NUT installation method: from source tarball, package or Subversion, portage * exact device name and related information (manufacturing date, web pointers, ?), Eaton 3S 550/700 * complete problem description, with any relevant traces, like system log excerpts, and driver debug output. You can obtain the

I'm attempting to run 250 permutations of a negative binomial GLM model for data on fish counts. Many of the models are fit appropriately, but some issue warnings such as "convergence not reached" or "step size truncated due to divergence." I've attempted using glm.nb from the MASS package and the negbin function from the AOD package, but both still cause some models to

This patch -Displays[1] current parameters for running domains ala xm list -Allow users to set[2] one or more parameters[3] without having to provide values for parameters they do not wish to change -Adds additional testing of sched-sedf via new xm-test testcases. With this patch applied, test 02_sedf_period_lower_neg.py exposes a bug[4]. I''ll follow up this email with a patch for

I have a large Excel file with SKU numbers (stock keeping units) and forecasts which can be mimicked with the following: Period <- c(1, 2, 3, 1, 2, 3, 4, 1, 2) SKU <- c("A1","A1","A1","X4","X4","X4","X4","K2","K2") Forecast <- c(99, 103, 128, 63, 69, 72, 75, 207, 201) PeriodSKUForecast <-

Dear R-list, I think it's really a newbie question but ... I try to model age-period-cohort models with polynoms in each effect. I have for each level of age and period some cases -k- and persons-years -py-. Models are A_xP_y: log(k_age,period/py_age,period)=f(polynom(x),age)+f(polynom(y),period). For exemple I try to adjust an A2P3 :

Hi, I am looking for some guidance on whether I can use the period(.) in function names and variable names. For example: my.function.name <- function(my.data.variable, my.radius, my.another.var, my.value = 10) { } Will this pose any problems regarding older and current version of R. -- SmartG [[alternative HTML version deleted]]

9d1fd58ff602 was bogous in not letting a new domain being created if its scheduling parameters --when running under the sedf scheduler-- were not fully specified, making creation fail like in this example here below: 2012-06-16 07:37:47 Z executing ssh ... root@10.80.248.105 xl create /etc/xen/debian.guest.osstest.cfg libxl: error: libxl.c:3619:sched_sedf_domain_set: setting domain sched sedf: