similar to: Vector comparison to matrix

I'm more like a caveman when it comes to programming tools. So, with that in mind, is there a way to use readBin in a batch format to read in pieces of a large binary file? Thank you for the consideration of my question. Todd Remund

I am looking for a Kalman filter that can handle a control input. I thought that l.SS was suitable however, I can't get it to work, and wonder if I am not using the right function. What I want is a Kalman filter that accepts exogenous inputs where the input is found using the algebraic Ricatti equation solution to a penalty function. If K is the gain matrix then the exogenous input

I have obtained from transfer functions, the state space matrices for the following state space model. x* = Ax + Bu y = Cx + Du I have A, B, C, and D, now I would like to take the exogenous inputs and simulate the data using the state space model. I know there is a simulate function in the package dse1, but I am unsure as to what type of TSmodel to create to put into it. Could anyone give me

I have had a strange problem with rgl. I am currently switching from a desktop computer to a workstation. When I get a very large data set and run spectro3D, the computations run then the window with the 3d plot appears it doesn't show the data inside and immediately shuts down R. The computers have OpenGL on them, but a colleague suspects the configuration of rgl with OpenGL may be the

I’m ordering a new computer to increase my ability to handle large data sets. I’ve tried the dual core type and also the dual processor with dual cores each, and have not been satisfied. This seems to agree with all the other postings on the help list. I don’t want to do any simulations, I just want to deal with a large data set and do things such as resampling time series data from 50,000

I'm trying to create a 3D plot using wireframe with certain parts removed. I would like to get rid of the part of the outer cube that crosses over the plot leaving the back two walls and the axes. It would also be useful to put lines in the plot separate from the wireframe call. I've looked through the documentation and have not been able to find such a request. Here is my code

Dear R help I was hoping you might be able to show me how to write a loop function take would ccomplish this task. library(prob) { a <- sample ( 1:20, 100, replace=T) b<-sample(5:24,100,replace=T) } dd <- data.frame(a,b) dd # code piece I am looking for if(subset(dd,c(1,23,ordered=F))is found))( print subset) else( continue evaluating subsets) subset(dd,isin(dd,c(1,23), ordered =

Hello I am trying to compare a vector with a Matrix's rows.The vector has the same length as the number of columns of the matrix, and I would like to find the row numbers where the matrix's row us the same as the given vector. What I am doing at the moment is using apply as follows: apply(Matrix,1,function(x)all(x%in%LHS)) but this isn't too fast actually. I would like to know if

I have looked through the documentation and have not been able to find a way of using an xlim, ylim, or zlim type option on rgl.surface. I know that persp3d has the option, but seems to only be able to expand the axes not reduce them. Is there anyone who has an idea of how to do this? Thank you for your time. [[alternative HTML version deleted]]

I know that there is a function, (spectro3D), that produces the Power Spectrogram. Are there R functions that produce the Frequency Spectrogram and the Phase Spectrogram? Thank you for your time. [[alternative HTML version deleted]]

I would like to ensure that all the text in my plot is in the arial font, like in excel and word. I have been saving my graphics to a file using the png() function, and have tried the postscript() function. In the documentation I have not been able to find a way of changing the family argument to an arial font. Does it exist? How would I go about creating a png or ps graphic using Arial?

I have tried the following code to plot and save many box plots in files. The code works by itself, but does not run in a loop. The error indicates that there is no device. The error is printed below the code. Code: for(i in 1:len) { A = as.numeric(delta[Delta[i,1]]) B = as.numeric(delta[Delta[i,2]]) C = as.numeric(delta[Delta[i,3]]) D = as.numeric(delta[Delta[i,4]]) data

I have fit a model with a more complex covariance structure, but the fit reports a false convergence. I have read from past posts that this can be an indication of over-specification. I went ahead and fit a model with a simpler covariance structure. It doesn't seem like I can compare the two likelihoods or the AIC or BIC to compare the two model since the one model had false convergence.

Hello, I have a problem with a LDAP backed based Samba PDC. Last week, due an hardware problem, I lost my primary LDAP server and PDC. I reinstalled the LDAP server and populated it with the old data, I also reinstalled Samba. The problem is that I can't log in to samba as root (cn=root,dc=info,dc=uvt,dc=ro). All others user accounts ar working except root. Eg.:

Dear All, I have encountered a strange bug that occurs (I think) with Unicode related file names in samba. Let me illustrate the example: I have a Win2k box with IP 192.168.0.2. This box has an NTFS folder called "Music" on it, which is shared for only one user. I have a Debian Unstable box with IP 192.168.0.1 which coincidentally acts as the router for the Win2k box. Both

Hello, I think I've found a bug in the C function do_grep located in src/main/grep.c. It seems to affect both the latest revisions of R-2-13-branch and trunk when compiling R without optimizations and with it's own version of pcre located in src/extra, at least on ubuntu 10.04. According to the pcre_exec API (I presume the later versions), the ovecsize argument must be a multiple of 3 ,

Hi R users, I have a question. How can I use for loop to do pair comparisons. For example, > x<-c(1,2,3) > result<-matrix(data=NA, nrow=choose(3,2), ncol=1) > for(i in 1: length(x)) + { result[i,]<-ifelse(x[i] > x[i+1], yes="Big", no="Small") + result} > result [,1] [1,] "Small" [2,] "Small" [3,] NA

Hello list, I'm doing a bootstrap analysis where some models occasionally fail to converge. I'd like to automate the process of restricting AIC to the models that do converge. A contrived example of what I'd like to do is below: resp <- c(1,1,2) pred <- c(1,2,3) m1 <- lm(resp~pred) m2 <- lm(resp~poly(pred,2)) m3 <- lm(resp~poly(pred,3)) # Fails, obviously ## Some

the man page for relational operators (see, e.g., ?'<') says: " Binary operators which allow the comparison of values in atomic vectors. Arguments: x, y: atomic vectors, symbols, calls, or other objects for which methods have been written. " it is somewhat surprizing that the following works: '<'(1) # logical(0) '<'() #

the man page for relational operators (see, e.g., ?'<') says: " Binary operators which allow the comparison of values in atomic vectors. Arguments: x, y: atomic vectors, symbols, calls, or other objects for which methods have been written. " it is somewhat surprizing that the following works: '<'(1) # logical(0) '<'() #