similar to: extending code to handle more variables

This is a continuation from a previous posting of mine: The following algorithm below is what I want to accomplish: Z(xk) = Average(Yi, i belongs to Ik), where Ik contains all i such that for each j, |Xi,j - xkj?? 2. Here, j = 1, 2 and i corresponds to the elements in each X and/or xk >data x1 x2 y 1 1 2 2 2 6 3 3 12 Now, consider a second data frame or matrix (xk):

Consider a data frame (x) with 2 variables, x1 and x2, having equal values. It looks like: x1 x2 1 1 2 2 3 3 Now, consider a second data frame (xk): xk1 xk2 0.5 0.5 1.0 0.5 1.5 0.5 2.0 0.5 0.5 1 1.0 1 1.5 1 2.0 1 0.5 1.5 1.0 1.5 1.5 1.5 2.0 1.5 0.5 2 1.0 2 1.5 2 2.0 2 I have written code to calculate some differences between these

I keep getting the 'memory.size' error message when I run a program I have been writing. It always it cannot allocate a vector of a certain size. I believe the error comes in the code fragement below where I have multiple arrays that could be taking up space. Does anyone know a good way around this? w1 <- outer(xk$xk1, data[,x1], function(y,z) abs(z-y)) w2 <- outer(xk$xk2,

useR's, Consider some variables and a data frame of points: x1 <- c(1,2,3) x2 <- c(3,4,5) xk1 <- seq(min(x1)-.5, max(x1)+.5,.5) xk2 <- seq(min(x2)-.5, max(x2)+.5,.5) expand.grid(xk1=xk1,xk2=xk2) xk1 xk2 1 0.5 2.5 2 1.0 2.5 3 1.5 2.5 4 2.0 2.5 5 2.5 2.5 6 3.0 2.5 7 3.5 2.5 ... 46 2.0 5.5 47 2.5 5.5 48 3.0 5.5 49 3.5 5.5 I want to restrict the data frame to only contain

Unfortunately I was unable to build s390 binaries since multilib was dropped in Fedora 24 on s390x. Going from the source of the 'file' command it seems as if it prints "32-bit" (the architecture is really 31 bit). --- generator/actions.ml | 12 ++++++++++++ src/filearch.c | 10 ++++++++++ test-data/Makefile.am | 2 ++

--==M2004020923123623794 Content-Type: text/plain; charset=us-ascii Content-Transfer-Encoding: 7bit --- The message cannot be delivered to the following address. --- anna@metso.com Permanent error involving remote host. 554 Transaction failed (too many hops) --==M2004020923123623794 Content-Type: message/delivery-status Reporting-MTA: mailgate01@metso.com Final-Recipient:

Here is a drop-in replacement for the R incomplete beta function. src/math/pbeta.c It is a slightly modified version of the cephes library one from Netlib. In the few cases I tried it seems to give at least 14 digit agreement with the one in S-PLUS (its hard to get more). I'm not sure what performance is like. I'd like to know if it helps with some of the problems which have been

Dear R list, Can anyone help with a plotting question? I'm trying to display some data on a plot and I've almost got the format I need (see code below), but 2 things I can't get: 1. How to get "Jan","Feb","Mar" on the x=axis instead of 1:3? 2. How to get "T"s on the end of my error bars like you have in standard scientific plots? Any comments

useR's, I am trying trying to find out if there is a faster way to do a certain computation. I have successfully used FOR loops and the apply function to do this, but it can take some time to fully compute, but I was wondering if anyone may know of a different function or way to do this: > x [1] 1 2 3 4 5 > xk [1] 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 I want to do:

--xk+LtGr+JJ Content-Type: text/plain; charset=us-ascii Content-Transfer-Encoding: 7bit Attached is a conversion script for converting current R doc files to the proposed new format ... I'm not sure if it handles everything correctly, as I made just some random checks. The main purpose for posting the thing to the list is such that those interested can get a look at the new format. Of

Hi everybody, I've been looking for a function that combines all variables from a dataset because I need to do multivariate regression. If we have linear regression with an expression like f(x) = a0 + sum(ai*xi) what I want to do is something like f(x) = a0 + sum(ai*xi) + sum(sum(bij * xi * xj)) + sum(sum(sum(cijk*xi*xj*xk))) + ... So I need a function that combines all the values from

useR's, I would like to know if there is a way to avoid using FOR loops to perform the below calculation. Consider the following data: > x [,1] [,2] [,3] [1,] 4 11 1 [2,] 1 9 2 [3,] 7 3 3 [4,] 3 6 4 [5,] 6 8 5 > xk Var1 Var2 Var3 1 -0.25 1.75 0.5 2 0.75 1.75 0.5 3 1.75 1.75 0.5 4 2.75 1.75 0.5 5 3.75 1.75

I have the following problem. Consider a `classical' test which works for k .ge. 2 samples. Possible interfaces are e.g. xxx.test(x, g) x ... all data, g ... corresponding groups xxx.test(x1, ..., xk) xxx.test(list(x1, ..., xk)) etc etc. Clearly, the first and the second one are nice, but cannot be combined without making `g' (i.e., `group') a named argument. Hence, in

Dear WizaRds, I am struggling with the use of twophase in package survey. My goal is to compute a simple example in two phase sampling: phase 1: I sample n1=1000 circuit boards and find 80 non functional phase 2: Given the n1=1000 sample I sample n2=100 and find 15 non functional. Let's say, phase 2 shows this result together with phase 1: ...................phase1........

Hi, I'm using samba 3.0.6. I see computer_name folders created in homedir. This never happened before I switch from Domain to ADS. Could someone please give me a few pointers? Regards, Norman Zhang [global] workgroup = ARKONDOMAIN realm = HQ.ARKONNETWORKS.COM server string = Samba Server %v security = ADS obey pam restrictions = Yes

Dear readers of the list, I have a problem a comparison of two data from a vector. The comparison yields FALSE but should be TRUE. I have checked for mode(), length() and attributes(). See the following code (R2.10.0): ----------------------------------------------- # data vector of 66 double data X =

Hi folks I am using the lmer function (in the lme4 library) to analyse some data where individuals are clustered into sets (using the SetID variable) with a single fixed effect (cc - 0 or 1). The lmer model and output is shown below. Whilst the fixed effects are consistent with stata (using xtmixed, see below), the std dev of the random effect for SetID is very very small (3.5803e-05)compared to

Dear R experts, I'm trying to find a code to calculate the p-value from the score test for the logistic regression. My fit is like this: logit=beta0+beta1*x1+beta2*x2 +....+ betak* xk. And my H0 is beta1=beta2=...=betak =0. Any help will be highly appreciated. Thank you! Ying

Not directly related to Asterisk but I'm sure one or two of you will get hot and bothered over this. :-) http://deancollinsblog.blogspot.com/2008/10/thecus-n7700.html Regards, Dean Collins Cognation Inc dean at cognation.net <mailto:dean at cognation.net> +1-212-203-4357 New York +61-2-9016-5642 (Sydney in-dial). +44-20-3129-6001 (London in-dial).

Hello folks With a small delay, the September release of xk-c is there. Have fun! Cheers Sergey ========================================= Akshay S Dinesh (2): typofix tel-sarala key layout Whitelist Indian keyboard layouts. #101532 Baurzhan Muftakhidinov (2): Enable level3 ralt_switch by default in Kazakh (ext) layout. Extend more keys in Kazakh (ext) layout. Jonas