similar to: Call for suggestions

Hello, I'm trying to do some differential equation modeling in R using the package 'deSolve.' Briefly, I'm trying to use the law of mass action (the details of which aren't really important) to structure a vector of rate equation which will be passed into an ODE function and solved with associated functions from 'deSolve.' Roughly what this entails is scanning

Hi everyone, I am trying to compare proportions among groups using the logistic regression approach as follows: 1) Fit the model log(p_i/(1-p_i)) = M + G_i, where p_i is the probability of success in group i and G_i is the effect of group i, i=1,..,I. 2) Test the hypotheses: Ho: G_1 = G_2 = ... = G_I (the probability of success is the same for all groups) versus Ha: at least two

Hi, is there a R function that order a matrix according to some criteria based on the rows(or cols) of that matrix? For example, let's say that my matrix S is composed by n rows S_1, S_2,.., S_n and that I compute some real value g_i=g(S_i) for each row. Then I want to order this set of g_i (from smaller to bigger) and order the correspondent row to the new position. Is it possible (apart

Here is a programming question. The code I am using is quite slow and I was wondering if there is a way to get rid of the for loop. I am dealing with "interaction" in 2x2 table, and am using Edwards's G_I (Likelihood, p. 194). I label the cells in the table as follows stim response "y" "n" total -------------------------------- y hit miss nsignal

Hi, I would appreciate if someone could help me on track with this problem. I want to compute some parameters from a system of equations given a number of sample observations. The system looks like this: sum_i( A+b_i>0 & A+b_i>C+d_i) = x sum_i( C+d_i>0 & C+d_i>A+b_i) = y sum_i( exp(E+f_i) * ( A+b_i>0 & A+b_i>C+d_i) = z A, C, E are free variables while the other

Greetings. For R gurus this may be a no brainer, but I could not find pointers to efficient computation of this beast in past help files. Background - I wish to implement a Cramer-von Mises type test statistic which involves double sums of max(X_i,Y_j) where X and Y are vectors of differing length. I am currently using ifelse pointwise in a vector, but have a nagging suspicion that there is a

Hi all, I want to do plot() in a loop to make 10 graphs, so I have some code like for (i in 1:10) { plot(... ... , xlab = expression(g[i]) ) } I expect g_1, g_2, and so on appear on x labels, but it simply prints g_i for each graph. Does anybody know how to get around this problem? Thanks. NL

I check the Mips backend for the following C code fragment compile result. It seems not correct. Is it my misunderstand or it's a bug. //ch8_3.cpp #include <stdarg.h> int sum_i(int amount, ...) { int i = 0; int val = 0; int sum = 0; va_list vl; va_start(vl, amount); for (i = 0; i < amount; i++) { val = va_arg(vl, int); sum += val; } va_end(vl);

Dear R users, I'm trying to fit a simple random intercept model with a fixed intercept. Suppose I want to assign a weight w_i to the i-th contribute to the log-likelihood, i.e. w_i * logLik_i where logLik_i is the log-likelihood for the i-th subject. I want to maximize the likelihood for N subjects Sum_i {w_i * logLik_i} Here is a simple example to reproduce

Which part of the generated code do you think is not correct? Could you be more specific? I compiled this program with clang and ran it on a mips board. It returns the expected result (21). On Tue, Feb 19, 2013 at 4:15 AM, Jonathan <gamma_chen at yahoo.com.tw> wrote: > I check the Mips backend for the following C code fragment compile result. > It seems not correct. Is it my

Hello :-) There was probably no good reason to do this build but the Debian 6.0 Squeeze repo version was 1.2.3, we needed 1.2.4 and I didn't think of using the package from unstable. Arch: x86_64 Platform: Linux 2.6 Debian 6.0 (Squeeze) Compiler: gcc version 4.4.5 (Debian 4.4.5-8) Version: 1.2.4 Status: no known problems Source: http://oligarchy.co.uk/xapian/1.2.4/xapian-core-1.2.4.tar.gz

On 28/10/11 16:30, xapian-discuss-request at lists.xapian.org wrote: > Date: Thu, 27 Oct 2011 16:29:26 +0100 > From: James Aylett <james-xapian at tartarus.org> > Subject: Re: [Xapian-discuss] Apache conf.d file for omega? > To: Charles <xapian at catcons.co.uk> > Cc: Xapian Discussion <xapian-discuss at lists.xapian.org> > Message-ID:

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

Hello., I configured Xapain in local machine (ubuntu 11.04). The configurations which i did are as follows Xapian Installation - 1) Xapian-core-1.2.7 - tar xvzf xapian-core-1.2.7.tar.gz - cd xapian-core-1.2.7 - ./configure --prefix=/usr/local - make - make install 2) Xapian-binding-1.2.7 - tar xvzf xapian-binding-1.2.7.tar.gz - cd xapian-binding-1.2.7 - ./configure --prefix=/usr/local

Dear R useRs, i try to compute the posterior mean for the parameters omega and beta for the following posterior density. I have simulated data where i know that the true values of omega=12 and beta=0.01. With the function postMeanOmega and postMeanBeta i wanted to compute the mean values of omega and beta by numerical integration, but instead of omega=12 and beta=0.01 i get omega=11.49574 and

Here's some code to make lissajous dance. I've attached a small sample GIF. Cheers, Rob Steele robsteele at yahoo dot com plot.lissajous = function(omega.x, omega.y, delta = 0, num.thetas = 200) { thetas = seq(0, 2 * pi, length = num.thetas) xs = sin(omega.x * thetas + delta) ys = cos(omega.y * thetas) plot(xs, ys, type = 'l', lwd = 3, ann = FALSE, axes = FALSE) }

NT4.0 SP4 gives me this error when trying to access a Samba 3.0.23c server. \\Omega300\files is not accessable The specific network password is not correct. omega300 is the samba server, running samaba 3.0.23.c on SCO Openserver 5.0.7. omg180 is a Windows NT4 SP4 box. omega is a Windows 98 box. Passwords are not encrypted. We are also using VisionFS (SCO's smb product, which they no longer

Hello, There is a build error by what looks to be a typo in xapian-applications/omega/omega.cc My checkout is about 15 minutes old, so it's possible it has been caught since then. omega.cc: In function ?int main(int, char**)?: omega.cc:204: error: ?pretty_tery? was not declared in this scope make[3]: *** [omega.o] Error 1 make[3]: Leaving directory

Dear all, I am a newbie to Xapian-Omega. I installed Xapian-Omega (1.2.14) a few days ago on my freebsd NAS server, and intend to use it to search my CIFS/SMB share. So far it works fine (I am able to search my share using the omega web interface), nevertheless I have few newbie questions. Hoping that somebody can answer my (hopefully not to silly) newbie questions. Indexing file names:

Cologne, 10.12.00 Dear Sir, Dear Madam, We are concerned about the possibility that the European Commission might introduce software patenting into the European Community because we think this will harm our profession. We make our living on data management and statistical analyses. Modern statistics crucially depends on algorithms [cf. Venables, W. N., & Ripley, B. D. (1999). Modern applied