search for: kirsten

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.

...mp$Subject_ID[i]]) } I am interested to know... (1) Is there a way to incorporate a weighted median into this script, where the weights are the number of times each distance is recorded? (2) Can I transform my current matrix into one that gives each distance its own row? Help is much appreciated, Kirsten Beyer

...= 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 the vector (i.e., the same line repeated with the next yyyy-mm-dd and corresponding julian day). Any help would be much appreciated. Thank you, Kirsten Kirsten Barrett Mendenhall Postdoctoral Fellow, USGS Alaska Science Center 4210 University Drive Anchorage, AK [[alternative HTML version deleted]]

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, > 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))) &...

I have a list of ICD9 (disease) codes with various formats - 3 digit, 4 digit, 5 digit. The first three digits of these codes are what I am most interested in. I would like to either add zeros to the 3 and 4 digit codes to make them 5 digit codes or add decimal points to put them all in the format ###.##. I did not see a function that allows me to do this in the formatting command. This seems

...the warning message: In max(i) : no non-missing arguments to max; returning -Inf Best as I can tell the code is working properly and the output is as expected. I would like some help in understanding why he is getting this error message and what its implications are. I have his code. Sincerely, Kirsten Miles Support Specialist Research Computing Lab Charles L. Brown Science and Engineering Library kdm3d@virginia.edu [[alternative HTML version deleted]]

...I keep getting an error statement about an invalid ). I've tried several permutations to fix, but without luck. Also, can anyone help me bind columns together? I've tried: cbind[, c(15:25)] but get the error: object of type 'closure' is not subsettable Thank you in advance! Kirsten [[alternative HTML version deleted]]

...ols model with good fits. I am running into trouble combining the two (as outlined in the third step of the Fletcher et al. 2005 paper). Does anyone have any experience or any advice on doing this? How does one come up with an overall model for the data using this approach? Thanks for your help! Kirsten -- View this message in context: http://r.789695.n4.nabble.com/Conditional-model-in-R-tp4651188.html Sent from the R help mailing list archive at Nabble.com.

...ly 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 estimates. Does this make more sense? Thanks for your time! Kirsten On Sun, Aug 6, 2017 at 3:21 PM Ben Tupper <btupper at bigelow.org> wrote: > 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 th...

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&l...

Hello. I am looking for a function that will allow me to paste rows together without regard for the numbers of columns in the datasets to be joined. The only columns where it matters if they are aligned correctly are at the beginning - the rest of the columns represent differing numbers of ICD9 (disease) codes reported by each person(record) at a health visit. They are in no particular order.

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...

Hi all, Thank you for taking the time to read my message. I'm trying to make a figure that plots p-values by a range of different adjustment values. (Using the **logit** function in package **car**) My Statistical analyses were conducted on probability estimates ranging from 0% to 100%. As it's not ideal to run linear models on percentages that are bounded between 0 and 1, these

...ed to add to the key.title line (size, color, font, etc) seems to be making a bit of difference. key.title = title(main = "Asymptomatic\n Prevalence", cex = "0.75") I assume there's an override earlier in the code, but I have no idea what it is. Any suggestions? Thanks! Kirsten -- View this message in context: http://r.789695.n4.nabble.com/Legend-Truncated-Using-filled-contour-tp4644938.html Sent from the R help mailing list archive at Nabble.com.

...e other than the page linked above. I've searched all the help forums and archives I can find, and I'm coming up with nothing. Is this some common function that I've never heard of, being self taught? Does anyone know where I can find documentation on this function? Thanks so much! Kirsten -- Kirsten Simmons, MPH Polymath interested in productivity, how ideas spread, gardening, marketing, entrepreneurship and models of social networks and disease transmission [[alternative HTML version deleted]]

...opulation (400). From 101-200 take 100 samples that are between 5% and 15% of the population (800). From 201-300 take 100 samples that are between 5% and 15% of the population (300). I assume this would require a nested for loop. Does anyone have advice as to how to do this? Thanks for your time. Kirsten ## Mark-Recapture ## Estimate popoulation from repeated sampling ## Population size N <- 400 N ## Vector labeling each item in the population pop <- c(1:N) pop ## Lower and upper bounds of sample size lower.bound <- round(x = .05 * N, digits = 0) lower.bound ## Smallest possible sample...

...effectively another word for"function"... functions like cbind are called with argument lists delimited by parentheses, not brackets, and having a missing argument to the cbind function will be of no use. -- Sent from my phone. Please excuse my brevity. On July 24, 2017 5:23:57 AM PDT, Kirsten Morehouse <kmoreho1 at swarthmore.edu> wrote: >Hi everyone, > >I'm having some trouble with my ifelse statements. > >I'm trying to put 12 conditions within 3 groups. Here is the code I >have so >far: > >dat$cond <- ifelse(test = dat$cond == "cond1&qu...

...ength 1 219.457 219.457 606.4349 < 2.2e-16 *** Year 1 2.970 2.970 8.2073 0.004934 ** Length:Year 1 1.012 1.012 2.7952 0.097173 . Residuals 119 43.064 0.362 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Thanks very much for any help. Kirsten [[alternative HTML version deleted]]

I would like to code records in a dataset with a 1 if any of the columns 9-67 contain a particular code, and zero if they don't. I've been working with "subset" and it seems that something like subset(data, data[9:67]--"12345") would work, but I have been unsuccessful so far. It seems like a simple problem - any help is appreciated!

...s it identifies the boxes. My question is whether or not anyone has ever tried anything like this, or if you even think it can be done. Having spent a week on this, I'm starting to doubt my initial assumption that it should be an easy alteration to make. Thanks for any input you might have, Kirsten Simmons -- Kirsten Simmons, MPH Polymath interested in productivity, how ideas spread, gardening, marketing, entrepreneurship and models of social networks and disease transmission [[alternative HTML version deleted]]