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

I'll second Bert's comments, also assuming this is not homework. In addition: Your use of "mid-point" is not a standard one (in my world), nor perhaps is that of "simulate". Let me attempt to re-state your problem: You wish to choose 10-tuples of integers 0 <= k <= 100 satisfying 1) 0 <= k_i <= 100 for i = 1:10; and 2) k_(i+1) - k_(i) >= 5 for i =

JRG: I don't think your specification is correct -- perhaps just a thinko. I think a 10-tuple of "reals" (scare quotes because of computer precision) with your specifications is what is wanted. Bert "An educated person is one who can entertain new ideas, entertain others, and entertain herself." On Sat, May 31, 2025 at 3:43?PM JRG <jrg at loesl.us> wrote: >

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

If this is a real problem and not homework, can you tell us the context? It is not at all clear (to me) what you mean by "simulate", i.e. what your target distribution is, which may depend on/be defined by the context. Bert "An educated person is one who can entertain new ideas, entertain others, and entertain herself." On Sat, May 31, 2025 at 11:52?AM Brian Smith

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