similar to: many weighted means: is there a simpler way?

Hi R users, I would like to add year information to each point in a scatter plot. The following example shows what i mean: df=data.frame(year=c(2001,2002,2003),a=c(8,9,7),b=c(100,90,95)) df plot(b~a,data=df) Now, I would like to see which point belongs to 2001, 2002 and 2003 in the above graphic. Thank you very much, ozan [[alternative HTML version deleted]]

Dear R Users, I have the following data frame: v1 <- c(rep(10,3),rep(11,2)) v2 <- sample(5:10, 5, replace = T) v3 <- c(0,1,2,0,2) df <- data.frame(v1,v2,v3) > df v1 v2 v3 1 10 9 0 2 10 5 1 3 10 6 2 4 11 7 0 5 11 5 2 I want to add a new column v4 such that its values are equal to the value of v2 conditional on v3=0 for each subgroup of v1. In the above example, the

Hi, Has there been any development in the OVF support for Xen? It was in the roadmap of Xen 3.4 (http://www.xen.org/download/roadmap.html) /Ozan _______________________________________________ Xen-devel mailing list Xen-devel@lists.xensource.com http://lists.xensource.com/xen-devel

Hi, Has there been any development in the OVF support for Xen? It was in the roadmap of Xen 3.4 (http://www.xen.org/download/roadmap.html) /Ozan _______________________________________________ Xen-devel mailing list Xen-devel@lists.xensource.com http://lists.xensource.com/xen-devel

Hi As you may''ve seen from my last posts I''m struggling to get the setup on my subject line to work without flaws. I''ve tried compiling the Kernel from Source, I''ve also tried Fedora 14, which has a flaky pv_ops Dom0 kernel which DOESN''T work at all, but none of these seem to work. I''m a mid-time XEN user, I just love it and want to stay

Hi As you may''ve seen from my last posts I''m struggling to get the setup on my subject line to work without flaws. I''ve tried compiling the Kernel from Source, I''ve also tried Fedora 14, which has a flaky pv_ops Dom0 kernel which DOESN''T work at all, but none of these seem to work. I''m a mid-time XEN user, I just love it and want to stay

Hi, As far as I am aware, the model.matrix function does not return perfect metadata on what each column of the model matrix "means". The columns are named (e.g. age:genderM), but encoding the metadata as strings can result in ambiguity. For example, the dummy variables created when the factors var0 = 0 and var = 00 both are named var00. Additionally, if a level of a factor variable

Dear All, I am economist and working on poverty / income inequality. I need descriptive statitics like the ratio of education expentitures between different income quintiles where each household has a different weight. After a bit of google search I found 'Hmisc' and 'quantreg' libraries for weighted quantiles. The problem is that these packages give me only weighted quintiles;

Hello, I am looking at the job satisfaction data below, from a problem in Agresti's book, and I am not sure where the degrees of freedom come from. The way I am fitting a binomial model, I have 168 observations, so in my understanding that should also be the number of fitted parameters in the saturated model. Since I have one intercept parameter, I was thinking to get 167 df for the Null

Hi, can anyone tell me how to nest two fixed factors using glmer in lme4? I have a split-plot design with two fixed factors - A (whole plot factor) and B (subplot factor), both with two levels. I want to do GLMM as I also want to include different plots as a random factor. But I am interested on the effect of A a B and their interaction on the response variable. I tried

Hi, I was wondering if there is an easy way to accomplish the following in R: Often I want to apply a function, e.g. weighted.quantile from the Hmisc package to grouped subsets of a data.frame (grouping variable) but then I also need to hand over the weights which seems not possible with summaryBy or aggregate or the like. Is there a function to do this? Currently I do this with loops but it

I've been running a glmer logit on a very large data set (600k obs). Running on a 10% subset works correctly, but for the complete data set, R completes apparently without error, but does not display the results. Given these jobs take about 200 hours, it's very hard to make progress by trial and error. I append the code and the sample and complete output. As is apparent, I upgraded R

for s10u2, documentation recommends 3 to 9 devices in raidz. what is the basis for this recommendation? i assume it is performance and not failure resilience, but i am just guessing... [i know, recommendation was intended for people who know their raid cold, so it needed no further explanation] thanks... oz -- ozan s. yigit | oz at somanetworks.com | 416 977 1414 x 1540 I have a hard time

Dear forum members, How can I force a negative slope in a linear regression even though the slope might be positive? I will need it for the purpose of determining the trend due reasons other than biological because the biological (genetic) trend is not positive for these data. Thanks. Julia Example of the data: [1] 1.254 1.235 1.261 0.952 1.202 1.152 0.801 0.424 0.330 0.251 0.229

I have a binary file with data sequence in the order [age,weight][age,weight] .... I know the length of the data and I want to load it into a data.frame. of course a way to do this is to read age and weight seperately and then use cbin(age,weight) to combine them into a dataframe, but is there a better solution? Thanks in advance!

Hi. I''m building a setup where I have users and groups, and a user can have an optional title in a group (eg. "owner"). I''m doing this using "has_many :through" associations, using a Role model object as the join hub: class User < ActiveRecord::Base has_many :roles, :dependent => :destroy has_many :groups, :through => :roles end class

Dear R users, I am looking for a simple way to define a contrast in a linear model. I have a data frame with two factors (f1, f2) and a dependent y. x <- data.frame(y =rnorm(30), f1=gl(2, 15), f2=gl(3, 5, 30)) Now I want to specify the following contrast: "f1= 1 or 2 and f2=1" vs. "f1= 1 or 2 and f2=3" The best I can come up with is the following: x$new[x$f1==1 &

When talking to Joel and Brian a while ago, I realized that it is a very common occurance to create an instruction, then immediately insert it into a basic block. At the time, this required keeping iterators around (updating them as neccesary) to determine where to insert the instruction, and it also required calling getInstList().insert(....) or something similar after creating the instruction

> 1. copyRegToReg() won't be able to use register classes to pick a copy opcode. For instance, an XMM register will no longer be copied by MOVSS or MOVSD. Given just the physical register, MOVAPS will be used. Is that a problem? I haven't had time to really look into it, but have been playing around with the idea that instead of two register classes copyRegToReg and some of the load

Signed-off-by: Micha? Miros?aw <mirq-linux at rere.qmqm.pl> --- drivers/gpu/drm/cirrus/cirrus_drv.c | 23 +---------------------- 1 file changed, 1 insertion(+), 22 deletions(-) diff --git a/drivers/gpu/drm/cirrus/cirrus_drv.c b/drivers/gpu/drm/cirrus/cirrus_drv.c index 69c4e352dd78..85ed8657c862 100644 --- a/drivers/gpu/drm/cirrus/cirrus_drv.c +++ b/drivers/gpu/drm/cirrus/cirrus_drv.c @@