similar to: State space model

Dear R users, When running the program below I receive the following error message: fit <- optim(parm, objective, yt = tyield, hessian = TRUE) Error in as.vector(data) : no method for coercing this S4 class to a vector I can't figure out what the problem is exactly. I imagine that it has something to do with "tyield" being a matrix. Any help on explaining what's going on

Dear all The following is a very basic and beginner's question on loops. Suppose you have data (say, 1000 cases) with two variables named "answer" (string) and "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

Hi, Can anyone tell me how to construct a simple VAR(1) time series with two variables using the dse package? I would like to end up with two time series y_1t = \phi_11 y_1,t-1 + \phi_12 y_2,t-1 + e_1t y_2t = \phi_21 y_1,t-1 + \phi_22 y_2,t-1 + e_2t Best regards, Sam.

Hallo everybody, I'm new in r and I"ll appreciate some help! I have a matrix of nrow=30 and ncoll=54,and I would like to generate 50 simulations with tha same size of the matrix!!!That is to say that I want to generate 50 matrices -for my 50 simulations - with the same dimensions! I took my 1st matrix according to the formula that I want to implement: D<-mean_m + U_i*mat_DELTA

Assume that we collect below data : - subjects = 20 males + 20 females, every single individual is independence, and difference events = 1, 2, 3... n covariates = 4 blood types A, B, AB, O http://r.789695.n4.nabble.com/file/n4245397/CodeCogsEqn.jpeg ?m = hazards rates for male ?n = hazards rates for female Wm = Wn x ?, frailty for males, where ? is the edge ratio of male compare to female Wn =

The officially sanctioned way to put the expression "lambda_1 = x" in a title is something like this: title(substitute(lambda[1] == lamb, list(lamb = x))) but suppose I have two lambdas and would like something like "lambda_1 = x_1 , lambda_2 = x_2" to appear. What then? Undoubtedly I'm missing something blindingly obvious with lists, but having tried several

Hi all; I recently have been used 'maxLik' function for maximizing G2StNV178 function with gradient function gradlik; for receiving this goal, I write the following program; but I have been seen an error  in calling gradient  function; The maxLik function can't enter gradlik function (definition of gradient function); I guess my mistake is in line ******** ,that the vector  ‘h’ is

For some small sample sizes I would like to exemplify the combinatorics underlying certain nonparametric tests like Mann-Whitney-Wilcoxon, Kruskal-Wallis and Spearman's rank correlation. I have written a function all.perm which delivers all permutations of 1,2,...,n and which works reasonably well. This can be used to generate P(R=r) of Spearman's rank correlation:

Hello R-experts, I wonder whether any of the R-packages cover the Hilbert-Huang Transform methodology (HHT)? Regards, Torbjorn -- Torbj?rn Lorentzen | torbjorn.lorentzen at bjerknes.uib.no |torbjorn.lorentzen at uni.no | http://www.bjerknes.uib.no/ Phone: +47 55 58 25 05 | Cellphone: +47 906 972 36 | Bjerknes Centre for Climate Research | Geophysical Institute | University of Bergen |

Dear R users, I am trying to carry out MLE of the time-varying CAPM using the FKF package. My approach so far has been to try and adapt the example given in the help file found using ?fkf which demonstrates the MLE of an ARMA(2,1) model. When I attempt to run my R code (given below) I get the following error: Error in fkf(a0 = sp$a0, P0 = sp$P0, dt = sp$dt, ct = sp$ct, Tt = sp$Tt, : Some of

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$

List members, Via searches I've seen similar discussion of this topic but have not seen resolution of the particular issue I am experiencing. If my search on this topic failed, I apologize for the redundancy. I am attempting to generate random covariance matrices but would like the corresponding correlations to be uniformly distributed between -1 and 1. The approach I have been using is:

Hi all, what would be the best way to install windows xp pro over XEN. I guess its only possible with NFS,HHT OR FTP shares. Does anyone have idea about it? How to have these share so that i could have Windows xp as a guest OS over xen(FC6)? Thanks -- View this message in context: http://www.nabble.com/installing-xen-over-XEN-tf2794004.html#a7794684 Sent from the Xen - User mailing list archive

Rainer M Krug writes: > Thorsten Wiegand used in his paper Wiegand T., and K. A. Moloney 2004. > Rings, circles and null-models for point pattern analysis in ecology. > Oikos 104: 209-229 a statistic he called O-Ring statistic which is > similar to Ripley's K, only that it uses rings instead of circles. > > http://www.oesa.ufz.de/towi/towi_programita.html#ring

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

Dear R-help, Could someone please try to explain this paradox to me? What is more likely to show up first in a string of coin tosses, "Heads then Tails", or "Heads then Heads"? ##generate 2500 strings of random coin flips ht <- replicate(2500, paste(sample(c("H", "T"), 100, replace = TRUE), collapse =

Hi, If I set df=2 in my smooth.spline function, is that equivalent to running a linear regression through my data? It appears that df=# of data points gives the interpolating spline and that df = 2 gives the linear regression, but I just want to confirm this. Thank you, Steven

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

Hi, I would like to know whether there exist algorithms to compute the coefficients or, at least, the degree of the minimal polynomial of a square matrix A (over the field of complex numbers)? I don't know whether this would require symbolic computation. If not, has any of the algorithms been implemented in R? Thanks very much, Ravi. P.S. Just for the sake of completeness, a

Hi Everybody, please can you help me how to call BayesX in R in order to see the graph already exist in BayesX Thanks ---------- Forwarded message ---------- From: <r-help-request@r-project.org> Date: Sat, Jan 22, 2011 at 5:00 AM Subject: R-help Digest, Vol 95, Issue 22 To: r-help@r-project.org Send R-help mailing list submissions to r-help@r-project.org To subscribe or