similar to: pattern matching

Hello Listers, I am trying to sample a vector to create a new one of sample length, witha sum equal to the sum of the initial vector: initial = 10, 30, 10 (sum=50) sample example = 5, 35, 10 (sum=50) or 25, 15, 10 (sum=50), etc ... My problem is to control the sum, so it stays constant. Any suggestions would be very helpful ! Thank you in advance, Eric

Hi Folks, I''m trying to setup a test envoriment which is composed by a puppetmaster running in my laptop (Macos Snow Leopard) and a puppet client running on a EC2 instance at Amazon. In order to allow the client to talk with the master I''m using an SSH Remote port forward, ie I login into from the laptop into the EC2 instance witha "-R 8140:localhost:8140" flag. This

Hi, I am simulating an event that has 15 possible outcomes and I have a vector 'pout' that gives me the probability of each outcome - different outcomes have different probabilities. Does anyone know a simple way of simulating the outcome of my event? If my event had only two possible outcomes (0/1) I would pick a uniform random number between 0 and 1 and use it to choose between the two

First off, Marc Schwartz posted this link earlier today, read it. http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-are-p_002dvalues-not-di splayed-when-using-lmer_0028_0029_003f Second, your email is not really descriptive enough. I have no idea what OR is, so I have no reaction. Third, you're comparing estimates from different methods of estimation. lmer will give standard errors that

dear all, I am writing a sweave documentation for my analysis, and I am plotting huge scatter plot data for microarray. unlucly this take a lot of resource to my pc because of the quality of the image which is to high (I see the PC get stuck for each single spot). how can I overcome this problem? is there a way to make lighter image? john [[alternative HTML version deleted]]

Hello All, I've decided to go ahead and plug my fax machine into one half of a Sipura SPA-2000. I've Google the list archive for examples of default fax contexts but not come up witha clearly working result. Can someone here please send me a config fragment that theyt know to work? Thanks, Michael -- Michael Graves mgraves@pixelpower.com Sr. Product

Hello, Does anybody know if Samba is compadable with Vista Yet? If so, is there a way for a group pr user policy to allow samba users to things in windows vista, like install programs w/o admin access. I know Vista has policies that can do this, where it cna basicly grant the user access, witha prompt or w/o a prompt. Can Samba support this too?

Hi, Another weird thing I'm seeing on my logs is this message about registering witha rootdse. ldb: ldb_tdb: Unable to register control with rootdse! ldb_wrap open of secrets.ldb My server works "normally", altough I'd like to know what this means. Thanks. --- David Gonzalez H. DGHVoIP - OPEN SOURCE TELEPHONY SOLUTIONS Phone Bogot?: +(57-1)289-1168 Phone Medellin:

dear all, we am trying to improve the performance of my R code, with the implentation of some function with custom C code. we found difficult to import and export/import data structure such us matrices or data.frame into the external C functions. we already tried the solution from "Writing R Extensions" form the R webpage. do you have any other solution or more advanced documentation

I am running some simulations which involve fitting a glm with Binomial error ( in fact Bernoulli) and logistic link. Running R 0.63 I have no problems. However I need more speed so migrated to our UNIX boxes which run Version 0.90.1 (December 15, 1999) Here I get an error Error in (if (is.empty.model(mt)) glm.fit.null else glm.fit)(x = X, y = Y, : y must be univariate unless binomial

I have a number of files containing anywhere from a few dozen to a few thousand integers, one per record. The statement "refdata18 = read.csv("K:\\MerchantData\\RiskModel\\Capture.Week.18.csv", header = TRUE,na.strings="")" works fine, and if I type refdata18, I get the integers displayed, one value per record (along with a record number). However, when I try "

I am considering replacing my old checkpoint and watchguard firewalls witha single Linux box using iptables and shorewall. I have two ISP''s (with separate routing tables), two DMZ''s, at least one VPN to a remote office, and a local trusted network. The configuration will look like: +----------------+ | | net0 ----------+ eth1

Hi, Can anyone shed some light on this: My PC's are on a separate power circuit. Each PC has its own UPS. I have an old Honda generator. When the power goes out, I fire up the generator. However, the UPS's still run on battery. If I connect the PC directly to the socket and bypass the UPS, the PC works. I was told by a UPS supplier, that the power output from the generator was

On Tue, 2004-08-03 at 14:05, Vikram Adve wrote: > I have been at Microsoft the last couple of days and so couldn't > jointhe discussion earlier. No worries. I knew you'd chime in sooner or later :) > Here's my view of the name issue, and (thereason this is long), a > little about how I think we want users to viewthis tool: > > First, I think the name should

I am using R to make two-way ANOVA on a number of variables using g <- aov(var ~ fact1*fact2) where var is a matrix containing the variables. However the outcome seem to be dependent on the order of fact1 and fact2 (i.e. fact2*fact1) gives a slightly (factor of 1.5) different result. Any ideas why this is? Thanks for any help Anders

dear all, I have a very simple question but I could not figure out. I need to make plots with grid in the background. something like I old retrive like this a=runif(100)*10 b=runif(100)*10 plot(a,b, pch=20, xlim=c(0, round(max(a))), ylim=c(0, round(max(b)))) vs=seq(0, max(a), 0.5) for(i in 1:length(vs)){ abline(v=vs[i], col="lightgrey") } hs=seq(0, max(b), 0.5) for(i in

Why the error is coming? even though the length of outcome.new$compkey and outcome.new$armkey were exactly same. Can anyone help? setwd("D:/AZ") library("RODBC") cdb_cnct <- odbcConnectExcel("AZIF_DC_GVK_NSCLC_MSALL_287papers_02072012_141450_v1_4.xls") outcomes <- sqlFetch(cdb_cnct, "Outcomes_info") odbcClose(cdb_cnct) rm(cdb_cnct)

Dear R users, I get the following (I think puzzling) result when doing the following: > a <- list(3,4,5) > a[[2]] <- NULL > a [[1]] [1] 3 [[2]] [1] 5 I would have expected the result to be: [[1]] [1] 3 [[2]] NULL [[3]] [1] 4 as in the outcome of: > list(3, NULL, 4) Is this a desired effect? If so, could it be built in a 'help(NULL)' file? If you think it is

Dear R users, I have a linear model of the kind outcome ~ treatment + covariate where 'treatment' is a factor with three levels ("0", "1", and "2"), and the covariate is continuous. Treatments "1" and "2" both have regression coefficients significantly different from 0 when using treatment contrasts with treatment "0" as the

Hi, Could someone help with coding this in R? I need to select one row per patient i in clinic j. The data is organized similar to that shown below. Two columns - patient i in column j identify each unique patient. There are two columns on outcome. Some patients have multiple rows with each row representing one visit, coded for in the column, visit. Some patients have just one row indicating