similar to: Summation in R

Hi there, I have a multiple for loop over a list of data frames for ( i in 1:(N-1) ) { for ( j in (i+1):N ) { for ( p in 1:M ) { v_i[p] = alist[[p]][i,"v"] v_j[p] = alist[[p]][j,"v"] } rho_s = cor(v_i, v_j, method = "spearman") rho_p = cor(v_i, v_j, method = "pearson"

Hi, I've got an array, say with i,jth entry = A_ij, and a vector, say with jth entry= v_j. I would like to multiply each column of the array by the corresponding vector component, i,e. find the array with i,jth entry A_ij * v_j This seems so basic but I can't figure out how to do it without a loop. Any suggestions? Michal.

Hello all, Is there a way in R to compute the multivariate normal density of every pair of entries in a vector efficiently instead of using for loop? For example Suppose I have a vector a=c(v_1,...,v_p)=c(0.5343909, -0.7784353, -0.0568370, 1.8772838, -1.3183407, 0.8227418,...) I want to compute density(v_i, v_j) for every pair of entries (i,j) (i!=j) in a. The joint bivariate distribution

On Wed, Jan 28, 2015 at 7:50 AM, Hal Finkel <hfinkel at anl.gov> wrote: > Hi Vijender, > > Thanks for posting this, there is wide support here for improving our support for reductions of various kinds, both in flavor and robustness. I've cc'd some others who have previously discussed this. > > James has advocated in the past for an intrinsic for horizontal reductions,

HI, When we run wine <program.exe>, it looks like wine by default tries to pick up files from drive_c\windows\system32. What do we do if the executable files are installed in a different location?. Wine cannot load it and fails. Is there a PATH variable or something we can set for WINE, just like we set in a OS? Cheers

Please see the following code for the runtime comparison between split() and mysplit.data.frame() (they do the same thing semantically). mysplit.data.frame() is a fix of split() in term of performance. Could somebody include this fix (with possible checking for corner cases) in future version of R and let me know the inclusion of the fix? m=300000 n=6 k=30000 set.seed(0) x=replicate(n,rnorm(m))

Please see the following code for the runtime comparison between split() and mysplit.data.frame() (they do the same thing semantically). mysplit.data.frame() is a fix of split() in term of performance. Could somebody include this fix (with possible checking for corner cases) in future version of R and let me know the inclusion of the fix? m=300000 n=6 k=30000 set.seed(0) x=replicate(n,rnorm(m))

On Thu, 4 Feb 2010, Conrad Parker wrote: >> Would there be any interest in hosting this alongside the current >> (definately write-only) java library? > > of course, the more better implementations the better! what vcs are you > using? Currently the vcs isn't anywhere public. I've dropped a snapshot at http://gagravarr.org/misc/vj-20100204-1.tar.gz in case

M P?opklpnlbyj?nvnmm M. ?plppkbkv??k? knbnnnnnn??????llp?????tx hikkkhgxxx vj jul?l Sent from my iPhonejukuj?b?mjl jnmnmmm Sorry for keeping things short Gustaf Granath (phd) Plant Ecology Uppsala University

Hello, can anybody here explain what a Jonckheere-Terpstra test is and whether it is implemented in R? I just know it's a non-parametric test, otherwise I've no clue about it ;-( . Are there alternatives to this test? thanks for help, Arne

Hi, I need a simple tool for benchmarking my NT client against my Ubuntu Samba server. I get 20MB/s during a file copy but I think my limit is the client's hard disk. Dose anyone know a simple exe tool which simulate file I/O agains the Samba server. I like to avoid using my lokal disk during the test. Have looked around but can't find a simpe tool which dose the job without the need of

Hi list, i wonder if this is possible with some dovecot configuration. Selected mailboxusers will only mark their mails as deleted, when a mail is being deleted through any client. No automatic expunge. This is default behaviour in Outlook, and i think in Thunderbird too, but unfortunately not in Open-Xchange for example. So i was asking myself if it would be possible to simply turn this

Full_Name: Dina Rego Version: R Cocoa GUI 1.14 OS: Intel Mac 10.4.4 Submission from: (NULL) (24.199.98.228) I recently bought the new intel imac and every program I had on my old computer (shareware and all) works on my new computer except R. When you double click it, it seems like it's going to launch (it does that box rushing to the screen thing) but nothing happens. I need this program

Dear all Version: samba 3.0.21a Server: Windows 2000 mixed (R2K-PDC) steps -- net rpc join domain -U administrator%xxxxx , (success joined) -- start winbindd => OK -- wbinfo -u => failed, error message kerberos_kinit_password RO120$@R2K-PDC.COM failed: Cannot resolve network address for KDC in requested realm ads_connect for domain R2K-PDC failed: Cannot resolve network

Hi all, Anyone experienced in the LaTeX format? I'm trying to use the xtable package to create nice anova tables, but how do I do to produce a pdf from the resulting LaTeX table? I've tried WinShell and MiKTeX, but I couldn't get any of them working... Here's an example of the output in R: % latex table generated in R 2.9.2 by xtable 1.5-6 package % Tue Nov 24

Hello, I was looking at the following test case which is very relevant in imaging applications. int sad(unsigned char *pix1, unsigned char *pix2) { int sum = 0; for( int x = 0; x < 16; x++ ) { sum += abs( pix1[x] - pix2[x] ); } return sum; } The llvm IR generated after all the IR

Dear R experts I need to write a function that incorporates double summation, the problem being that the upper limit of the second summation is the index of the first summation, i.e: sum_{j=0}^{x} sum_{i=0}^{j} choose(i+j, i) where x variable or constant, doesn't matter. The following code obviously doesn't work: f=function(x) {j=0:x; i=0:j; sum( choose(i+j,i) ) } Can you help? Thanks

Dear R users... I made the R-code for this double summation computation http://www.nabble.com/file/p19213599/doublesum.jpg ------------------------------------------------- Here is my code.. sum(sapply(1:m, function(k){sum(sapply(1:m, function(j){x[k]*x[j]*dnorm((mu[j]+mu[k])/sqrt(sig[k]+sig[j]))/sqrt(sig[k]+sig[j])}))})) ------------------------------------------------- In fact, this is

Dear R users... I made the R-code for this triple summation computation http://www.nabble.com/file/p20517134/a.jpg ------------------------------------------------- Here is my code.. x=seq(.1,1,.1); l=10 y=seq(1,10); m=10 z=seq(.1,1,.1); n=10 sum(sapply(1:l, function(i) {sum(sapply(1:m, function(j) {sum(sapply(1:n, function(k){exp(x[i]*y[j]*z[k] )/gamma(y[j]+1)}))^(1.5) }))}))

Dear R users, I try to compute this summation, http://www.nabble.com/file/p25054272/dd.jpg where f(y|x) = Negative Binomial(y, mu=exp(x' beta), size=1/alp) http://www.nabble.com/file/p25054272/aa.jpg http://www.nabble.com/file/p25054272/cc.jpg In fact, I tried to use "do.call" function to compute each u(y,x) before the summation, but I got an error, "Error in X[i, ]