similar to: Nested for loop

Hi Kirsten, I can run your example code but I can't quite follow your division of sampling. Can you restate the the task? Below is what I think you are asking for, but I have the feeling I may be off the mark. Set A: 400 samples, draw 100 in range of 5 to 15 Set B: 800 samples, draw 100 in range of 5 to 15 Set C: 300 samples, draw 100 in range of 5 to 15 Ben > On Aug 5, 2017, at

Hi Ben, That's exactly right! Except for each set it's the sample population that is 400, 800 or 300. I want to take 3 samples, each of 100, where only the population differs. I can do this separately, but I'm having trouble putting them all on the same graph. I'd like to have sample on the x axis (1-300) and estimate on the y axis. I want to show how population affects the

Hi! Thanks for taking the time to read this. The code below creates a graph that takes 100 samples that are between 5% and 15% of the population (400). What I'd like to do, however, is add two other sections to the graph. It would look something like this: from 1-100 samples take 100 samples that are between 5% and 15% of the population (400). From 101-200 take 100 samples that are between

Thanks Caitlin and Richard MH. Works great. Stephen From: Caitlin [mailto:bioprogrammer at gmail.com] Sent: Friday, November 10, 2017 12:33 PM To: Bond, Stephen Subject: Re: [R] update R version in windows install.packages("installr") updateR() rather... On Friday, November 10, 2017, Bond, Stephen <Stephen.Bond at cibc.com<mailto:Stephen.Bond at cibc.com>> wrote: Is

When I try to use stubdom-dm in lieu of qemu-dm, the guest domain crashes almost immediately if I specify an emulated NIC (''vif='') in the guest configuration file. If I don''t specify a vif, it works fine. If I use the same configuration with qemu-dm instead of sutbdom-dm, it works fine. Is there some reason QEMU emulated networking doesn''t work with stub domains

Hello R-Community, we are running aprogram to fit Non-linear differential equations to Aphid population Data and to estimate the birth and death parameters, here is the code: dat<-data.frame(Time=c(0:60),Cur=c(5,6.2,59,39,38,44,20.4,19.4,34.2,35.4,38.2,48.2,55.4,113.2, 97,112,115,126,136.6,140.6,147.2,151.6,157.8,170,202,210.4,221.2,224.4,248.2,266,

Hello, I have a problem with using replace() to convert a vector of dates from yyyy-mm-dd to julian date. For example, I type replace(x,2004-05-14,134) and I receive an error: Error in `[<-.data.frame`(`*tmp*`, list, value = 134) : new columns would leave holes after existing columns If I can successfully convert, I have a script that will convert all of the dates in

Hi Jim, Thanks for your help, I really appreciate it. Perhaps I'm misunderstanding, but does this formula run different ajustment values for this function? logit(p = doc$value, adjust = 0.025) I'm looking to plot the p-values of different adjustment values. Thanks so much, Kirsten On Wed, Jul 12, 2017 at 8:49 PM, Jim Lemon <drjimlemon at gmail.com> wrote: > Hi Kirsten,

Okay, I'm at the end of my amazing internet research skills, and now I'm having to offer this problem up to the mailing lists. In a fit of stupidity I've just upgraded my machine from 98 to XP (purely for research ;), and now I need to add it to the domain. My root user on the Linux box running samba is called guB~b1n5 and the Linux box is configured to use PAM and SSL, I've

Hi All, I'm trying to create two side-by-side contour plots with one legend by modifying the code found here: http://wiki.cbr.washington.edu/qerm/sites/qerm/images/b/bb/Example_4-panel_v1a.R I've been able to set up the other two functions called in the code, but I can't find reference to print.letterTrevor.R anywhere other than the page linked above. I've searched all the help

Not a reproducible example, so a bit of guessing here, but a) don't try to assign results to variables inside the ifelse. That is, remove all the single-equals signs and "test" variables. If you really need to conditionally assign variables then use "if"... but chances are good you don't need that. b) "closure" is effectively another word

Have also now changed the root username back to root to see if that helped (although it doesn't cause any problems with any other applications), but it does now mean that I get the same error message using the root user as I do using the user seri: The following error occurred attempting to join the domain "HOME": Access is denied. I have also added the line nt acl support = no to

Dear R experts, I want to plot a line chart with another secondary axis placed right to the standard secondary axis which one can access with the axis command, so that the data lines are seen in the same plot. Is there any way to do this in R? Many thanks, Kirsten.

Hello, I've installed the Becker and Wilks maps, mapdata, and mapproj packages so I can begin to try these out for some work I need to do on a map of Alaska but I don't know where to find a map of Alaska. Has anyone solved this already and could help? Thanks very much in advance, Caitlin

Hi all. I need to construct a plot showing words on the x-axis and how many times each word was given as a verbal response on the y-axis as solid bar (frequency). Is there a convenient function to do this in R? I considered hist(), but I'm not sure how to construct the text file. Example: apple, 2 pear, 14 house, 1 beach, 5 computer, 15 Thanks, ~Caitlin [[alternative HTML version

Hello, Looking over the swap_tail macro in pnorm.c, the comment inside the code indicates it is to swap the values of '*cum' and '*ccum'. According to the code below that is taken from the source file, the swap functionality is dependent on what values are provided for 'x' and 'lower'. #define swap_tail \ if (x > 0.) {/* swap ccum <--> cum */

Hi all. I need to plot 30 values and I'd like the x-axis to have 30 dates values, e.g., 5/1/2011, 5/2/2011, etc. so an observer could see how the values correspond as the month progresses. Is there a convenient function I could use? I considered creating two vectors, x and y, then simply plotting those. Thanks, ~Caitlin [[alternative HTML version deleted]]

Hi Kirsten, Perhaps this will help: set.seed(3) kmdf<-data.frame(group=rep(1:4,each=20), prop=c(runif(20,0.25,1),runif(20,0.2,0.92), runif(20,0.15,0.84),runif(20,0.1,0.77))) km.glm<-glm(prop~group,kmdf,family=quasibinomial(link="logit")) summary(km.glm) pval<-0.00845 padjs<-NA npadj<-1 # assume you have five comparisons in this family for(method in p.adjust.methods) {

sorry for the late reply. i dunno if anyone ever wrote you back. although perhaps not written for clarity, i think what appears in the code is correct. the routine only guarantees that one of cum or ccum has a valid return value depending on i_tail (the other might as well be garbage, but that isn't what's returned). j On Wed, 9 Jul 2003 Henrik.Seidel@schering.de wrote: >

Many thanks to Nick Ellis, Vadim Kutsyy, Charles Berry, and Bill Dunlap for providing thoughts and alternative solutions to the problem. I have included Nick and Charles' notes in full below and a summary of the others. Thanks also to Bill for telling me about an inconsistency in how the first argument of sample is interpreted. I had ignored this, resulting in a major bug. Each person