search for: qtr

...(pnorm (au_j,mean=b_j * a$R_m, sd= sqrt(sigma_j^2))- pnorm(al_j, mean=b_j * a$R_m, sd= sqrt(sigma_j^2) )), 1) ) ) ) ) } start.par = c(-0.01,0.01,0.1,1) #looping now runs=133/20+1 out <- matrix(NA, nrow = runs, ncol = 4, dimnames = list(paste("Qtr:", 1:runs , sep = ''), c("al_j", "au_j", "sigma_j", "b_j"))) ## Estimate parameters based on rows 0-20, 21-40, 41-60 of afull for (i in 1:runs) { index_start=20*(i-1)+1 index_end= 20*i a=afull[index_start:index_end,] out[i, ] <-...

...-rep(5,200) age<-c(age3,age4,age5) #combines ages into one vector unique<-1:1000 #gives each fish a unique id pop<-data.frame(unique,length,age) pop<-pop[sample(1:1000,size=1000,replace=FALSE),] #randomized the order of pop c.one<-pop[sample(1:1000,size=100,replace=TRUE),] a.one.qtr<-c.one[sample(1:100,size=25,replace=TRUE),] merge<-merge(c.one,a.one.qtr,by="unique",all=TRUE) What I would ultimately like to have is one row for all 100 in the sample and three columns (unique, length, age). And then some way to identify the 25 individual selected rows. Thank...

Hi All, I have a very simple question about a time series object: how to access values for a particular year and quarter (say)? Suppose, following http://www.stat.pitt.edu/stoffer/tsa2/R_time_series_quick_fix.htm I have read in data as a time series; here is how it looks. * Qtr1 Qtr2 Qtr3 Qtr4 1960 0.71 0.63 0.85 0.44 1961 0.61 0.69 0.92 0.55 . . . . . . . . . . 1979 14.04 12.96 14.85 9.99 1980 16.20 14.67 16.02 11.61* How do I access the value for 1961 quarter 3 (say)? Dipankar [[alternative HTML version del...

# Dear List, # I want to characterize a time series according to its Quarter components. # My data ("a.ts": http://docs.google.com/View?id=dfvvwzr2_478cr9k4cdb)? look like: #???????????????? Qtr1????????? Qtr2????????? Qtr3????????? Qtr4 #?? 1948 -0.0714961837? 0.0101747827? 0.0654816569 -0.0227830729 #?? 1949 -0.1175517556? 0.1151378692? 0.1015777858 -0.1971535900 #?? 1950? 0.0716002123? 0.2551020416? 0.0977574743 -0.0739337411 #??? ... #? The time series is 216 long #? The easiest way...

hi, i need help on placing legend/text in a time series plot. here is what i am doing (i am using rw1022 on windoze 2000): #read data file gdpn <- scan("jngdpsa.dat", list(year=0, qtr=0, gdp=0)); gdpr <- scan("jrgdpsa.dat", list(year=0, qtr=0, gdp=0)); #convert to time series object gdpn <- ts(gdpn$gdp, frequency=4, start=c(1955,2)); gdpr <- ts(gdpr$gdp, frequency=4, start=c(1955,2)); #plot two gdp time series in one graph postscript("fig2-1j.eps",...

The following data spans out to 2012 and when graphing out the data - R defaults dates on x axis only shows the years eg 2005 2007 2009 My question is how can I get the graph to show more dates on the axis - eg every qtr or 6 months or every month etc PRODUCT,DATE,AMOUNT_1,AMOUNT_2 A,3/1/2005,126647400,121942966.2 A,4/1/2005,121942966.2,122244045 A,5/1/2005,122244045,117227480.04 A,6/1/2005,117227480.04,117913720.32 A,7/1/2005,117913720.32,116767296. Regards Phil -- View this message in context: http://r.789695....

Folks, I want to make the following function more efficient, by vectorizing it: getCriterionDecisionDate <- function (quarter , year) { if (length(quarter) != length(year)) stop ("Quarter and year vectors of unequal length!"); ret <- character(0); for (i in 1:length(quarter)) { currQuarter <- quarter[i]; currYear <- year[i]; if ((currQuarter < 1) |

...t problem I have is that 4th quarter 2010 is missing--not even NA, there is no record for Year=2010 and Quarter =4, so instead of 260 rows, I only have 259. To solve this, I created a temporary DF with "Year" and "Quarter" that was complete Yr.temp = rep(1947:2011, rep(4,65)) Qtr.temp = rep(1:4, 65) Temp.df = data.frame(cbind(Yr.temp,Qtr.temp)) and merged the two, so now I have the NA's. so, my DF is gdpdata with the above four columns, 260 rows. My first question: what is the difference between gdp.ts <- ts(gdpdata$GDP, start=1947, end=2011, fr=4) and gdp2.ts &lt...

Hi everyone,I have a data daily data (x) for 10 years starting from 04-01-1995 to 03-31-2005. I was able to get the yearly sum for the ten years using aggregate(x, years, sum). But this gave me the yearly sum for 1995 (Apr- Dec); 1996 (Jan-Dec) ---------2005 (Jan-Mar). But I want to get the aggregates for Apr-1995 to Mar 1996, Apr 1996- mar 1997 and so on. your help will be higly appreciated.

Has anyone in the UK got ANI to work on an inbound call ? Using asterisk 1.4 trunk and zaptel 1.4 trunk, with a Euro-ISDN 30 Julian ______________________________________________________________________ This email has been scanned by the MessageLabs Email Security System. For more information please visit http://www.messagelabs.com/email

Hi all- I have a dataset (year, month, hour, co2(ppm), num1,num2) [49,] 2006 11 0 383.3709 28 28 [50,] 2006 11 1 383.3709 28 28 [51,] 2006 11 2 383.3709 28 28 [52,] 2006 11 3 383.3709 28 28 [53,] 2006 11 4 383.3709 28 28 [54,] 2006 11 5 383.3709 28 28 [55,] 2006 11 6 383.3709 28 28 [56,] 2006 11 7 383.3709 28 28 [57,] 2006

Hi, I have a ts object with a frequency of 4, i.e., quarterly data, and I would like to calculate the mean for each quarter. So for example: > ts.data=ts(1:20,start=c(1984,2),frequency=4) > ts.data Qtr1 Qtr2 Qtr3 Qtr4 1984 1 2 3 1985 4 5 6 7 1986 8 9 10 11 1987 12 13 14 15 1988 16 17 18 19 1989 20 If I do this manually, the mean for the 1st quarter would be mean(c(4,8,12,16,20)), which is 12. But I am wondering if there is a R function that c...

...t auth_verbose = yes Password file: # squirrelmail imap users user1:xe9l8sIMpcTsg:1008:1008::::chroot:Maildir:/home/cmp/Maildir/user1 user2:m7VKCgWVmZYMY:1008:1008::::chroot:Maildir:/home/cmp/Maildir/user2 user3:ELh9mthLVjEG2:1008:1008::::chroot:Maildir:/home/cmp/Maildir/user3 # normal users xyz1:QtR.iN1DbL/UM:1004:1004::::chroot:mbox:/var/mail/xyz1 xyz2:y38.qLIgrvYic:1002:1002::::~/Mail:mbox:/var/mail/xyz2:INDEX=/home/xyz2 ... Problem shown in log: Dec 19 16:45:30 localhost dovecot-auth: PAM: pam_authenticate(xyz2) failed: User not known to the underlying authentication module Dec 19 16:45:3...

I know this isn't the "Cisco list", but enough people here are wired into the VoIP world that perhaps someone has heard if Cisco has released a SIP image for the 7920 yet... JT

Dear List, I ran into some problems with time-series-Data. Imagine a data-structure where observations (x) of test attendants (i) are made a four times (q) a year (y). The data is orderd the following way: I y q x 1 2006 1 1 1 2006 3 1 1 2006 4 1 1 2007 1 1 1 2007 2 1 1 2007 3 1 1 2007 4 1 2 2006 1 1 3 2007 1 1 3 2007 2 1 I am looking for a way to count the attendants that at least have

Having problems with the autocomplete on my linux box and firefox. When I use form tags it breaks, without form tags everything works fine. Exact copy of the example provided, but I added form tags.. The initial search and select works but then the box sort of "locks up" and I can''t change it at all, If I do try to type something nothing changes but a second later the whole

I have two data frames. One (arr) has all arrivals to an airport for a year, and the other (gw) has the dates and quarter hour of the day when the weather is good. arr has a Date and quarter hour column. >names(arr) [1] "Date" "weekday" "hour" "month" "minute" [6] "quarter" "ICAO"

As preliminary steps in splitting virt-p2v to an own repository, continue making p2v more independent within libguestfs. This is accomplished by the following changes: - have only the authors in the about dialog, and read them from a local AUTHORS file - few more cleanups This is still not complete, although I believe most of the work needed is done, and it still makes p2v usable within

If the connection is not local, paths of disks will refer to the remote host, which were mistakenly handled as local paths (in the best case failing to open a non-existing disk, and in the worst case opening a different disk!). In case the disks are remote resources like ssh or ceph, nothing guarantees that the hostname can be reached from the local machine, or even that it is actually the same on