similar to: Referring on the value of a counter in a loop

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

Hi, I'm trying to estimate the parameters of a state space model of the following form measurement eq: z_t = a + b*y_t + eps_t transition eq y_t+h = (I -exp(-hL))theta + exp(-hL)y_t+ eta_{t+h}. The problem is that the distribution of the innovations of the transition equation depend on the previous value of the state variable. To be exact: y_t|y_{t-1} ~N(mu, Q_t) where Q is a diagonal

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:

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$

I am a beginner when it comes to using R, though fortunately I already know something about statistics. I think factor analysis should be used sparingly, but I occasionally use it. It doesn't seem to me that factanal() provides Kaiser's Measure of Sampling Adequacy, which should be computed for factor problems based on a small number of subjects, though perhaps it is elsewhere. Does anyone

Hello all, at the moment we are using samba 3.6.6 as pure fileserver on Debian Wheezy with external company-wide OpenLDAP as "passdb backend" (ldapsam). We want to use Debian Jessie, which have samba4 in stock, unfortunately out setup doesn't work with Samba4 anymore (NT_STATUS_UNSUCCESSFUL) I've read somewhere that Samba4 only works with internal LDAP server, is it correct?

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 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

I am sharing some printers attached to Win95 boxes, and would like to be able to do more than connect from Linux with smbclient and print via translate. I would like to be able to specify a printcap entry for these so as to be able to print from pine or lynx or Netscape (when I get X running here at work). I haven't found anything in Samba docs or LDP docs yet to tell me how to write such

Dear R-list members, I've had a hard time trying to solve a non-linear system (nls) of equations which structure for the equation i, i=1,...,4, is as follows: f_i(d_1,d_2,d_3,d_4)-k_i(l,m,s) = 0 (1) In the expression above, both f_i and k_i are known functions and l, m and s are known constants. I would like to estimate the vector d=(d_1,d_2,d_3,d_4) which is solution

Hello. I have searched the archives but have not found anything. I need to solve a constrained optimisation problem for a nonlinear function (“maximum entropy formalism”). Specifically, Optimise: -1*SUM(p_ilog(p_i)) for a vector p_i of probabilities, conditional on a series of constraints of the form: SUM(T_i*p_i)=k_i for given values of T_i and k_i (these are constraints on

Hi All, I''ve downloaded a Puppet module and I''m trying to add a parameter to it by editing lib/puppet/type/<resource>.rb. I simply added: newproperty(:pcfree) do desc "My description here" end In the corresponding file in lib/puppet/provider, I do something with :pcfree. Whenever I call the resource from Puppet with my new parameter, it keeps

Dear Community, does anybody have an idea on how to estimate a system of three seemingly unrelated regressions, two of which being TOBIT and one OLS? Background: I am currently estimating a translog cost function and two corresponding cost share equations using systemfit and the seemingly unrelated regression ("SUR") specification. However, I consider it more appropriate to estimate

Do you think I'd break a lot of people's filters if I removed the prefix? :) Anyone strongly for/against removing it? It seems kind of annoying to me whenever I happen to think about it. -------------- next part -------------- A non-text attachment was scrubbed... Name: signature.asc Type: application/pgp-signature Size: 198 bytes Desc: This is a digitally signed message part URL:

I am trying to see if this is even possible. We have a network in my office of mostly Windows machines, all connected to the internet via a firewall. We have a co-located Linux server in another state. I want to centralize my authentication setup, so that there is but a single username/password combincation used accross all systems. I have configured SSHD on the co-located Linux server to

Does anyone know how can I reset my printers when a job is removed? The problem is that when Win95 users remove a job from their print queue viewer, the print job gets successfully removed from the server, but the printer is leaved in a waiting state for more input. The only way is to turn off the printer, then turn it on again. Here goes my print command from smb.conf (the -h switch in the lpr

As most of you already know from a previous announcement, the R user conference useR! 2008 is scheduled for August 12-14, 2008, and will take place at the University of Dortmund. PRE-CONFERENCE TUTORIALS Before the official program, half-day tutorials will be offered on Monday, August 11th. We invite R users to submit proposals for three hour tutorials on special topics on R. The

As most of you already know from a previous announcement, the R user conference useR! 2008 is scheduled for August 12-14, 2008, and will take place at the University of Dortmund. PRE-CONFERENCE TUTORIALS Before the official program, half-day tutorials will be offered on Monday, August 11th. We invite R users to submit proposals for three hour tutorials on special topics on R. The

We are happy to inform you that the online abstract submission and registration for `useR! 2008' is now available online from http://www.R-project.org/useR-2008/ This meeting of the R user community will take place at the Technische Universit?t Dortmund, Dortmund, Germany, August 12 to 14, 2008. The conference schedule comprises invited lectures and user-contributed sessions as well as

We are happy to inform you that the online abstract submission and registration for `useR! 2008' is now available online from http://www.R-project.org/useR-2008/ This meeting of the R user community will take place at the Technische Universit?t Dortmund, Dortmund, Germany, August 12 to 14, 2008. The conference schedule comprises invited lectures and user-contributed sessions as well as