similar to: How to run regressions over increasing time series

I am trying to create basic descriptive statistics for a spatial dataset using moving windows of variable sizes. I found the filter(dataset, rep(1/a,1)) very useful for obtaining averages for a given moving window size "a" (a is the window size within a column of data in a matrix - the filter is applied to all columns automatically) I could not find anything analogous for

Hello, I am currently writing my Master's Thesis on a topic regarding the analysis of memory safety and termination of LLVM programs. This includes alignments in LLVM IR, but I am not sure if I understand their semantics correctly. I have written a program (see attachment) which uses the instruction store i32 1, i32* %7, align 4 to store an integer at an address that I forced to be

Hi all, This could be naivety/stupidity on my part rather than a problem with model output, but here goes.... I have fitted a fairly simple model m1<-glm(count~siteall+yrs+yrs:district,family=quasipoisson,weights=weight,data=m[x[[i]],]) I want to know if yrs (a continuous variable) has a significant unique effect in the model, so I fit a simplified model with the main effect ommitted...

Hi all: Using the example below, is there a way to add Y axis titles to each graphic instead of sharing the same title? library(ggplot2) RT = matrix(c(814, 500, 424, 394, 967, 574, 472, 446),4,2) colnames(RT) = c('repetition','alternation') rownames(RT) = c('7-yrs','11-yrs','15-yrs','21-yrs') rt <- melt(RT) names(rt) <- c("age",

hi all, i've taken over the maintainership of dovecot's pkgsrc on NetBSD - please send me your ideas, questions and comments related to dovecot on NetBSD. FreeBSD's port has currently some more 'features' than the pkgsrc, i'm going to add some of them - any feedback is appreciated. bye, TOM -------------- next part -------------- A non-text attachment was scrubbed...

Hello, I have a data frame containing time (e.g. GMT), and I would like to create/add a new variable that would be the computation of the elapsed time since the first observation. Does anyone have a suggestion for an easy way to do this? I am having trouble creating a new variable that would contain just the first time observation (then I could take difference between actual time and

Dear R-helpers, I have a data.frame (bcpe.lat.m) containing 13 countries, ages 0-50yrs per month, and the corresponding mu&sigma (see below). * I would like to limit the age range to include all 12 months for the 1st 5 years and only whole years for all ages thereafter for each of the countries present in the data frame. * I would like to create separate data.frames according

Folks, This question is somewhat related to a previous posting of mine. I just can't seem to create a generic solution. Here is a function that I found searching around the internet: splitIt <- function(x, n) {split(x, sort(rank(x) %% n))} I use it like so: > splitIt(1:12, 2) $`0` [1] 1 2 3 4 5 6 $`1` [1] 7 8 9 10 11 12 Or > splitIt(1:12, 4) $`0` [1] 1 2 3 $`1` [1] 4 5 6

I've helped a lot of people on the mailing lists and on IRC #asterisk. and wanted to let people know that I will be in Europe between May 19 and June 21. Stockholm (VON 2005), Brussels (holiday/vacation), Amsterdam (holiday/vacation), and Madrid (Astricon). There are several weeks during my trip that I have no current plans for and may add other cities to my itinerary. I'm looking

Hello, I am trying to match lat/lon from one dataset with the lat/lon from a second dataset and use that rows data for calculations. I am using match, but this is finding the first match and not comparing the pair, how can I determine if the lat/lon are the same? See example below. Is there a better way to determine to a matching pair of lat/lon values? Example Datasets: > data2

Folks, I call the function calcAmorts like so: calcAmorts(prevAm, amort, myDates) Note that I use the package lubridate. The last line where do.call is called to first divide all the rows by the first row and then rbind gives the following error: Error in do.call("rbind", apply(amortsByYears, 1, "/", amortsByYears[, : second argument must be a list By contrast if

Hi R-users, I have the following problem with CrossTable function within ?gmodels? package: the output of the function (format ?spss? and asresid=T) can not be stored within another object. For example: >library(gmodels) >data(infert, package = "datasets") > CrossTable(infert$education, infert$induced)->aa # the function prints everything ok on the screen > aa # works

I have been timing a particular model fit using lmer on several different computers and came up with a peculiar result - the model fit is considerably slower on a dual-core Athlon 64 using Goto's multithreaded BLAS than on a single-core processor. Here is the timing on a single-core Athlon 64 3000+ running under today's R-devel with version 0.995-5 of the Matrix package. >

Dear R Users,i am very new to R. I want your help on an issue regarding distance matrix and cluster analysis i had discharge data of 4 rivers(a,b,c,d) in 4 vectors each having 364 values > dput(qmu)structure(list(a = c(0.26, 0.25, 0.25, 0.25, 0.24, 0.23, 0.22, 0.21, 0.21, 0.21, 0.2, 0.19, 0.19, 0.19, 0.19, 0.18, 0.18, 0.18, 0.17, 0.17, 0.17, 0.17, 0.17, 0.17, 0.17, 0.17, 0.17, 0.17, 0.17, 0.17,

hi, I have been trying to do this in R (have implemented it in Excel) but I have been using a very inefficent way (loops etc.). I have matrix A (columns are years and ages are rows) and matrix B (columns are birth yrs and rows are ages) I would like to first turn matrix A into matrix B And then I would like to convert matrix B back again to the original matrix A. (I have left out details of

I have the following piece of code that combines lists comprised of components of varying length into a list with components of constant length. I have found 2 ways to do it, and the faster of the two is posted below along with sample results. Do you have any suggestions on how to decrease the calculation time by modifying the code? > ####Function########### >

greetings to one and all. a new iso burn w/ centos 6.8 failed to boot on a 686 32 bit mid tower box, but will boot a 32 bit laptop, so i dropped back to 6.7 which did boot. after install with centos 6.7, attempt to upgrade fails with 'repolist 0'. all desire repo files have enable = 1. searched thru past 4 yrs of personal archives, nothing found related to 'repolist 0'. ran

Hi all, I want to know how to perform the test Breslow-Day test for homogeneity of odds ratios (OR) stratified for svytable. This test is obtained with the following code: epi.2by2 (dat = daty, method = "case.control" conf.level = 0.95, units = 100, homogeneity = "breslow.day", verbose = TRUE) where "daty" is the object type table svytable consider it, but

we have got asterisk 1.0 (over 1 yrs old) version and very old zaptel version. That code is working only with 8 or less ports (accumulative) on digium fxs/fxo cards (2 cards with 4 ports each). the questoin is, what if we want 12 ports?..well, really, i don't understand the limitations? is it simply zaptel driver code fix? or kernel fix? or technology limitation? donno any tips would

hi list, first of all, many thanks to the developers of dovecot for this great piece of software! i'm running dovecot 0.99.9.1 on netbsd-current - runs very fine until i try to use sylpheed from my client to connect via ssl: --snip-- Jun 13 08:48:32 devbox imap-login: SSL_accept() failed: error:1408F455:SSL routines:SSL3_GET_RECORD:decryption failed or bad record mac [192.168.100.88] Jun 13