similar to: R sings

I am having a problem with the package SSPIR. The code below illustrates it. I keep getting the message: "Error in y - f : non-conformable arrays." I tried to tweak the code below in many different ways, for example, substituting rbind for cbind, and sometimes I get a different error message, but I could not find a variation of this code that would work. Any help will be greatly

I'm currently running into a little trouble with the kfilter method, and would love some clarification if you are able to offer it. When trying to run kfilter, I've been running into errors that seem to result from having mismatched dimensions. Specifically, the dimension of my observations is 2, while the dimension of the state space is 4. In the filterstep function (file sspir_kfs.R),

Hello, I apologize for such a basic question, but I have been trying to do this in multiple packages without much success. I am trying to set up a state space model for Kalman filtering. I am using package dlm. The DLM is specified by: observation: y(t) = F(t)*theta(t) + v(t) state: theta(t) = G(t)*theta(t-1) + w(t) I have no problem setting up a simple example where F is constant. I am

Hi all, please cc me, i'm not on the list. Second: All google findable information about problems setting up ntlm_auth for squid with winbind are read and checked more than three times. After breaking a running setup under debian squeeze, i go back to debian lenny to circumvent the actual MIT kerberos problem[1]. [1] http://bugs.debian.org/cgi-bin/bugreport.cgi?bug=566977#57 Now i

Hi all, I am working for nest box occupancy data for birds and would like to construct a Markov transition matrix, to derive transition probabilities for ALL years of the study (not separate sets of transition probabilities for each time step). The actual dataset I'm working with is 125 boxes over 14 years that can be occupied by 7 different species, though I have provided a slimmed down

Dear masters, I have data that looks like this: #val Freq1 Freq2 0.000 178 202 0.001 4611 5300 0.002 99 112 0.003 26 30 0.004 17 20 0.005 15 20 0.006 11 14 0.007 11 13 0.008 13 13 ...many more lines.. Full data can be found here: http://dpaste.com/173536/plain/ What I intend to do is to have a cumulative graph with "val" as x-axis with "Freq1" & "Freq2" as

Hello, I recently have started to learn about R, and am interested in using it to generate data that would then get used to produce audio. Basically, I am "envisioning" using R to create the data, based on probability functions (and fractal data), which then gets output as a text data file, and then processed by CSound to produce the finished audio. My interests are based on some of

Hi, I have loaded a dataset in R : data = label freq1 freq2 news 54 35 fun 37 21 milk 19 7 food 3 3 .... etc And I have a vector flist<-c("fun","food") Now I want to use the vector 'flist' for selecting these values from 'data' so that I get the following dataset : label freq1 freq2 fun 37 21 food

Dea'R' helpers I have following data - prob = c(0.1, 0.2, 0.3, 0.4, 0.5) frequency = c(100, 75, 45, 30, 25) no_trials = c(10, 8, 6, 4, 2) freq1 = rbinom(frequency[1], no_trials[1], prob[1]) freq2 = rbinom(frequency[2], no_trials[2], prob[2]) freq3 = rbinom(frequency[3], no_trials[3], prob[3]) freq4 = rbinom(frequency[4], no_trials[4], prob[4]) freq5 = rbinom(frequency[5],

I am getting the following error from MSVC when calling ov_open() from my app: Unhandled exception in Progname.exe (KERNEL32.DLL): 0xC0000005: Access Violation The call stack shows: KERNEL32! bff7b9a6() VORBISFILE! 00471d1d() I am running Win98 and v1.0 off the vorbis libs. The app is single- threaded. The following code works just fine in a test app I built, but when I add it to another

Dear R-Users, I'm using SSPIR package for a spatio-temporal application. Is it possible to modify the structure of the involved matrixes (Fmat, Gmat, Vmat,Wmat)? I want to create a model like this #y(t)=k*theta(t)+epsilon(t) #theta(t)=h*theta(t-1)+eta(t) #epsilon(t) N(0,V) V=sigma2*I #eta(t) N(0,W) W=sigma2_eta where the state variable theta has dimension 1(p=1) and at

Hi all, it's Advent time and I took the opportunity to program a song for you. The sound, i.e. speech synthesis and reverb, is rendered by the LLVM-JIT and the synthesizer is controlled via ALSA-MIDI by the Haskell live sequencer. http://www.youtube.com/watch?v=0EQCgi5qa3E "Alta trinita beata" You find additional information in the video description. Best, Henning

Dear list, I am trying to re-analyse something. I do have two time series, one of which (ts.mar) might help explaining the other (ts.anr). In the original analysis, no-one seems to have cared about the data being time-series and they just did OLS. This yielded a strong positive correlation. I want to know if this correlation is still as strong when the autocorrelations are taken into account.

Hi, I managed to make R core dump (linux and Mac OSX versions), but I think I've figured out why. First, here's the message I get on core dump (on linux - no message on Mac): R: ../inst/include/tnt/fmat.h:529: TNT::Vector<T> TNT::matmult(const TNT::Fortran_Matrix<T>&, const TNT::Vector<T>&) [with T = double]: Assertion `A.num_cols() == x.dim()'

Greetings, I sent the message below to the developer of the contributed R package "sspir", but have yet to receive any response. I would be very grateful for any advice people have on the matter. Thanks, Mark -------- Original Message -------- Subject: possible bug in sspir? Date: Tue, 19 May 2009 16:08:41 -0700 From: Mark Scheuerell <mark.scheuerell at noaa.gov> To:

Thank you for the suggestion, Bill. The result is not quite what I would like. Here's sample code for you or anyone else who may be interested: Al1 = c('A','C','C','C') Al2 = c('G','G','G','T') Freq1 = c(0.0078,0.0567,0.9434,0.9908) MAF = c(0.0078,0.0567,0.0566,0.0092) m1 = data.frame(Al1=Al1,

Hello R useRs, I have a function which returns a list of functions : freq1 <- function(x) { lev <- unique(x[!is.na(x)]) nlev <- length(lev) args <- alist(x=) if (nlev == 1) { body <- c("{", "sum(!is.na(x))", "}") f <- function() {} formals(f) <- as.pairlist(args) body(f) <- parse(text = body) namef <-

Hi, I have a grouped data set and would like to calculate weighted proportions for a large number of factor variables within each group member. Rather than using dplyr::count() on each of these factors individually, the idea would be to do it for all factors at once. Does anyone know how this would work? Here is a reproducible example: ############################################################

Optimize opus decode (float only) use case using ARM NE10. Mainly effects opus_ifft and ctl_mdct_backward and related functions. Work based on previous Encode optimization using ARM NE10 library. TBD: Add commit id of upstream Encode NE10 optimization patch so that users have reference of how to enable this optimization Signed-off-by: Viswanath Puttagunta <viswanath.puttagunta at

Optimize opus decode (float only) use case using ARM NE10. Mainly effects opus_ifft and ctl_mdct_backward and related functions. Work based on previous Encode optimization using ARM NE10 library. TBD: Add commit id of upstream Encode NE10 optimization patch so that users have reference of how to enable this optimization Signed-off-by: Viswanath Puttagunta <viswanath.puttagunta at