search for: k_1

Displaying 8 results from an estimated 8 matches for "k_1".

Did you mean: __1
2011 Sep 22
1
Error in as.vector(data) optim() / fkf()
...blem is exactly. I imagine that it has something to do with "tyield" being a matrix. Any help on explaining what's going on and how to solve this is much appreciated. Thank you, Kristian library(FKF) #loading Fast Kalman Filter package library(Matrix) # matrix exponential package K_1 = 0.1156 K_2 = 0.17 sigma_1 = 0.1896 sigma_2 = 0.2156 lambda_1 = 0 lambda_2 = -0.5316 theta_1 = 0.1513 theta_2 = 0.2055 #test data tyield <- matrix(data = rnorm(200), nrow =2, ncol =100) # defining dimensions m <- 2 # m is the number of state variables n <- 100 # is the length of the obs...
2011 Nov 12
1
State space model
...lt;- matrix(abs(rnorm(400)), nrow=10, ncol=40) m <- 2 # m is the number of state variables n <- ncol(x) # is the length of the observed sample d <- nrow(x) # is the number of observed variables. h <- 1/52 ## creating state space representation of 2-factor CIR model CIR2ss <- function(K_1, K_2, sigma_1, sigma_2, lambda_1, lambda_2, theta_1, theta_2, delta_0, delta_1, delta_2) { ## defining auxilary parameters, phi_11 <- sqrt((K_1+lambda_1)^2+2*sigma_1^2*delta_1) phi_21 <- sqrt((K_2+lambda_2)^2+2*sigma_2^2*delta_2) phi_12 <- K_1+lambda_1+phi_11 p...
2006 Jul 24
1
Referring on the value of a counter in a loop
...nd "class" (numeric). The latter runs from 1 to 78 and categorizes the data. I need to create a file "answer_1.txt" for the cases with class==1, and so on, up to "answer_78.txt" I have been able to do this for one value of "class" in the following way k_1<-answer[class==1] write.table(k_1, "answer_1.txt") Now, I would like to loop over "class", replacing the occurences of "1" in the code above successively with the integers from 1 to 78, so that I get my 78 files. (equivalent to typing: k_1<-answer[class==1]...
2011 Jul 26
1
Optimal Code for subset of Integer Sets
Hi, I wanted to know if there existed an good implementation in R of the following classical subset Integers : $P_{n,m} = {k_1, \ldots, k_m \in \mathbb{N} : k_1 + \ldots + k_m = n }$ for any integers $m < n$. There is an obvious not optimal code which would be to run through $m$ sums and put a boolean condition to report if the integers respect these condition but it would take $m^n$ steps... Thanks for your answers,...
2010 Oct 24
1
Optimize parameters of ODE Problem which is solved numeric
Hi, I have a data-matrix: > PID sato hrs fim health 214 3 4.376430 6.582958 5 193 6 4.361825 3.138525 6 8441 6 4.205771 3.835886 7 7525 6 4.284489 3.245139 6 6806 7 4.168926 2.821833 7 5682 7 1.788707 1.212653 7 5225 6 1.651463 1.436980 7 4845 6 1.692710 1.267359 4 4552 5 1.686448 1.220539 6
2000 May 03
0
Combinatorics for nonparametric tests
..."h") all.perm_function(n){ p_matrix(1,ncol=1) for(i in 2:n){p_pp_cbind(p,i) v_c(1:i,1:(i-1)) for(j in 2:i){v_v[-1] p_rbind(p,pp[,v[1:i]])}} p} Until now I hesitated to dwell on writing functions to generate all n!/(k!(n-k)!) k-combinations of an n-element set or more generally all n!/(k_1!k_2!...n_k!) choices of placing n objects into k boxes where box j has n_j objects. Not intending to re-invent the wheel I turn to the R community to ask if someone has already written these functions. Any help is thankfully appreciated. --- D.Trenkler --- ***********************************...
2004 Jul 12
6
proportions confidence intervals
Dear R users this may be a simple question - but i would appreciate any thoughts does anyone know how you would get one lower and one upper confidence interval for a set of data that consists of proportions. i.e. taking a usual confidence interval for normal data would result in the lower confidence interval being negative - which is not possible given the data (which is constrained between
2008 Jun 30
4
Rebuild of kernel 2.6.9-67.0.20.EL failure
Hello list. I'm trying to rebuild the 2.6.9.67.0.20.EL kernel, but it fails even without modifications. How did I try it? Created a (non-root) build environment (not a mock ) Installed the kernel.scr.rpm and did a rpmbuild -ba --target=`uname -m` kernel-2.6.spec 2> prep-err.log | tee prep-out.log The build failed at the end: Processing files: kernel-xenU-devel-2.6.9-67.0.20.EL Checking