Re: [R] predict.drm not generating confidence intervals

, at 5:19 PM, Peter Ehlers wrote: Brant, See below. On 2010-11-28 12:25, David Winsemius wrote: Puzzled. Why are the data you offer to predict() for the independent variable, conc, all NA's? Is there something reversed or inverted about how drc functions handle formulas. -- David. On Nov 28, 2010

Re: [R] Array help

,,] and StatTemps[7:1123:34,,] but that didn't work. Any ideas? Thanks so much! It may not be an array (since time values don't play very well with that data structure.) If there is time as an index, it may be a more complex object such as a time-series or zoo. Check with str(). -- David

Re: [R] Plot data inside matrix

do this with few commands?? Possibly depending on the correctness of my assumptions: xyplot.zoo -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] selecting only corresponding categories from a confusion matrix

perfect for me. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self

Re: [R] Troubles in plotting to a postscript file (not to png)

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org

Re: [R] data.frame and formula classes of aggregate

values) if there are any NAs. Wouldn't one expect that the 2 forms (data frame vs formula) of aggregate would give the same result? thanks very much david freedman, atlanta -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing

Re: [R] FW: how to use by() ?

/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

Re: [R] reference text variables as column name to plot

/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

[R] And about sum() ...Re: Minor warning about seq

, but was actually adding the fractional risk to each of the q[1:4]'s before summation occurred. (Should have been ... sum(q[1:4] , frac*q[5]) -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r

Re: [R] warning creating an as.array method in a package

on my system that it may be Matrix. showMethods(as.array) Function: as.array (package base) x=ANY x=Matrix -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] Can't Destroy Dim Names

10 ... - attr(*, dimnames)=List of 3 ..$ dim1: chr [1:3] a b c ..$ dim2: chr [1:3] d e f ..$ dim3: chr [1:3] g h i - Christian Langkamp christian.langkamp-at-gmxpro.de -- View this message in context: http://r.789695.n4.nabble.com/Can-t-Destroy-Dim-Names-tp876633p3066413.html David

Re: [R] Pass an operator to function

is a bit non-standard: ?do.call do_this - function(a, fn=, b) {do.call(fn, list(a , b))} do_this(a=1, b=4) [1] FALSE do_this(a=1, b=0) [1] TRUE -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch

Re: [R] How to mean, min lists and numbers

cumulative mean of a vector, i.e, c( mean(vec[1]), mean(vec[1:2]), ,,, mean(vec) ): vec - sample(1:20) sapply(1:length(vec), function(x) mean(vec[1:x]) So, how to mean/min a list of lists and numbers to return a list? Not a list and not working on a list of lists. A vector. -- David Winsemius

Re: [R] How to mean, min lists and numbers

, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self

Re: [R] in continuation with the earlier R puzzle

. -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained

Re: [R] How to select the column header with \Sexpr{}

-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

Re: [R] Comparison of two very large strings

commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] Comparison of two very large strings

On Jul 12, 2010, at 6:46 PM, David Winsemius wrote: On Jul 12, 2010, at 6:03 PM, harsh yadav wrote: Hi, I have a function in R that compares two very large strings for about 1 million records. The strings are very large URLs like:- http://query.nytimes.com/gst

Re: [R] Continuing on with a loop when there's a failure

this is not a homework assignment! -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented

Re: [R] Xyplot or Tin-R problem?

) #??? -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained

Re: [R] Continuing on with a loop when there's a failure

: David Winsemius dwinsem...@comcast.net To: Josh B josh...@yahoo.com Cc: R Help r-help@r-project.org Sent: Mon, July 12, 2010 8:09:03 PM Subject: Re: [R] Continuing on with a loop when there's a failure On Jul 12, 2010, at 6:18 PM, Josh B wrote: Hi R sages, Here is my latest problem. Consider

Re: [R] Continuing on with a loop when there's a failure

On Jul 13, 2010, at 9:04 AM, David Winsemius wrote: On Jul 13, 2010, at 8:47 AM, Josh B wrote: Thanks again, David. ...but, alas, I still can't get it work! (BTW, it did work.) Here's what I'm trying now: for (i in 1:2) { mod.poly3 - try(lrm(x[,i] ~ pol(x1, 3) + pol(x2, 3), data=x

Re: [R] Continuing on with a loop when there's a failure

0.3245 x2^2 -0.94888 0.82067 -1.16 0.2476 x2^3 0.06362 0.05098 1.25 0.2121 ...so what exactly would I query to design my if/then statement? From: David Winsemius dwinsem...@comcast.net To: David Winsemius dwinsem...@comcast.net Cc: Josh B josh...@yahoo.com; R Help r-help@r

Re: [R] Continuing on with a loop when there's a failure

On Jul 13, 2010, at 10:26 AM, David Winsemius wrote: On Jul 13, 2010, at 9:24 AM, Josh B wrote: In my opinion the try and tryCatch commands are written and documented rather poorly. Thus I am not sure what to program exactly. Didn't you see the silent parameter? Its seems

Re: [R] how to extract information from anova results

commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] Continuing on with a loop when there's a failure

if/then statement? From: David Winsemius dwinsem...@comcast.net To: David Winsemius dwinsem...@comcast.net Sent: Tue, July 13, 2010 9:09:04 AM Subject: Re: [R] Continuing on with a loop when there's a failure On Jul 13, 2010, at 9:04 AM, David Winsemius wrote: On Jul

Re: [R] sort file names in numerical order

://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] loop troubles

the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] loop troubles

a name worth remembering. --Theodore Roosevelt -Original Message- From: David Winsemius [mailto:dwinsem...@comcast.net] Sent: Sunday, July 18, 2010 5:36 AM To: Joe P King Cc: r-help@r-project.org Subject: Re: [R] loop troubles On Jul 17, 2010, at 9:09 PM, Joe P King wrote: Hi all, I

Re: [R] Identify points (was Plot error)

of a row i.e in my table: name value_1 value_2 bill 14 ben 2 2 jane 3 1 I click on the data point at 2,2 and it would read out ben Check out HWidentify or HTKidentify in pkg:TeachingDemos. -Peter Ehlers David Winsemius, MD West Hartford, CT

Re: [R] NA preserved in logical call - I don't understand this behavior because NA is not equal to 0

the logical vector which seems entirely wasteful or using subset(). David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org

Re: [R] Continuing on with a loop when there's a failure

) + pol(x2, 3), data=x) anova(mod.poly3)[1,3] [1] 0.6976063 Any ideas? From: David Winsemius dwinsem...@comcast.net To: Josh B josh...@yahoo.com Cc: Peter Konings peter.l.e.koni...@gmail.com; R Help r-help@r-project.org Sent: Sun, July 18, 2010 3:33:07 PM Subject: Re: [R] Continuing on with a loop

Re: [R] NA preserved in logical call - I don't understand this behavior because NA is not equal to 0

, the flaw is that R doesn't apply this philosophy entirely consistently: I think factor and table should default to exclude = NULL. Hadley -- Assistant Professor / Dobelman Family Junior Chair Department of Statistics / Rice University http://had.co.nz/ David Winsemius, MD West Hartford, CT

Re: [R] loop troubles

] sigma[i] = 1/sqrt(pn[i]) mutotals[i,]-mu.out[i] } Joe King 206-913-2912 j...@joepking.com Never throughout history has a man who lived a life of ease left a name worth remembering. --Theodore Roosevelt -Original Message- From: David Winsemius [mailto:dwinsem...@comcast.net] Sent

Re: [R] R: package plotrix

. It specifies what can be attached and TIFF files are NOT in the list of acceptable types. Uwe might have gotten a copy but the rest of us did not. -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch

Re: [R] specifying column names in a vector of characters and the use?

On Jul 18, 2010, at 10:09 PM, Seth wrote: Hi, What I would like to do is have a data.frame with column names and have these column names stored as strings in another vector. Then I would like to be able to access the data.fram columns via referencing the vector of names. The code

Re: [R] Grouping and stacking bar plot for categorical variables

/index.html __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West

Re: [R] invalid type error

-1.7 And the reason the ,drop=FALSE is needed is that otherwise the vector and the matrix would have been non-conformable after the column vector became a row vector. -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing

Re: [R] Indexing by logical vectors

__ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT

Re: [R] ACCTGMX to 1223400 in R?

in teh zoo and an strapply function in the gsubfn package that might be even more powerful, but I am insufficiently talented to give you a one-liner using them. Please advise. Thank you. -- -- David Winsemius, MD West Hartford, CT __ R-help@r

Re: [R] Constrain density to 0 at 0?

On Jul 19, 2010, at 9:57 PM, Farley, Robert wrote: I'm plotting some trip length frequencies using the following code: plot( density(zTestData$Distance, weights=zTestData$Actual), xlim=c(0,10), main=Test TLFD, xlab=Distance, col=6 ) lines(density(zTestData$Distance,

Re: [R] trouble getting table of coeffs with quantreg with fixed effects

are the fixed effects -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal

Re: [R] Indexing by logical vectors

. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained

Re: [R] Help with time in R

(or possibly used hte data.frame function) , you are probably experiencing the effects of the default stringAsFactors=TRUE behavior. Try setting stringsAsFactors to FALSE. Then you can use the standard Date functions: ?Dates -- David Winsemius, MD West Hartford, CT

Re: [R] Help with time in R

On Jul 20, 2010, at 8:41 AM, David Winsemius wrote: On Jul 20, 2010, at 7:33 AM, Sarah Chisholm wrote: Hi, I have a problem with the time formatting in R. I have entered time in the format MM:SS.xyz and R has automatically classified this as a factor, but I need it numerically. However

Re: [R] Means from selected columns in a data frame

was rather ambiguous) that function is also available. -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting

Re: [R] Means from selected columns in a data frame

On Jul 20, 2010, at 9:25 AM, David Winsemius wrote: On Jul 20, 2010, at 9:13 AM, Marcus Drescher wrote: Hi all, I have a dataframe with survey data. Now I want to calculate means from several but not all columns (e.g. a1, a2, a3) and save them in a new separate column (e.g. a). Well

Re: [R] Convert Unix (Epoch) timestamp to DD/MM/YY and HH:MM:SS

. You can find the format codes with: ?strptime print(Sys.time(), format=%m/%d/%Y %H:%M) # your request is ambiguous w.r.t month and day order [1] 07/20/2010 09:39 EDT -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list

Re: [R] Nesting functions in loops that result in error messages breaking the loop

) # or Design for (i in 1:3) { mod.poly3 - try(lrm(x[,i] ~ pol(x1, 3) + pol(x2, 3), data=x), silent=TRUE) if(class(mod.poly3)[1] != 'try-error') {results[1,i] - anova(mod.poly3)[1,3]} David Winsemius, MD West Hartford, CT __ R-help@r-project.org

Re: [R] Servreg $loglik

__ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT

Re: [R] If help

On Jul 20, 2010, at 1:14 PM, Heiman, Thomas J. wrote: Hi Y'all, I have some data in a table with 2 columns. There are two values: Reduction and No Reduction. I am trying to make a new variable change which recode the combinations from column 1 and 2 into a single number. Here is a

Re: [R] Registered / trademark signs

Works on OSX 10.5.8 + R2.11.1 -- David. On Jul 20, 2010, at 1:36 PM, Henrique Dallazuanna wrote: I don't know if this works in OSX, but in XP: plot(0) title('\u00ae - \u2122') On Tue, Jul 20, 2010 at 2:22 PM, Dennis Fisher fis...@plessthan.com wrote: Colleagues, What is the easiest

Re: [R] square brackets in expression in plots

__ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT

Re: [R] Limited output

What does this return? options()$max.print # [1] 9 # which is the default And what does this return? str(your.data.frame) -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] Constrain density to 0 at 0?

the histogram and then overlay a smooth fitted distribution curve. ?hist require(MASS) ?fitdistr -- David. -Original Message- From: David Winsemius [mailto:dwinsem...@comcast.net] Sent: Monday, July 19, 2010 19:31 To: Farley, Robert Cc: r-help@r-project.org Subject: Re: [R] Constrain

Re: [R] blank pdf output when called in a loop

://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] best way to apply a list of functions to a dataset ?

. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT

Re: [R] date from weeknumber

On Jul 20, 2010, at 6:37 PM, H Rao wrote: Hi, Is there a function to get the last(or first) day of the week, given the week number of the year? For eg, week number for 7/20 is 29 as obtained by format(Sys.Date(),%U), is there a function which returns 7/25 - the last day of week # 29

Re: [R] Obtaining the unmerged cases from one of the two data set

] NAME SALARY B.SALARY dim(common)[1] [1] 13 David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] Chi-square distribution probability density function:

. See also: ?Chisquare ?gamma -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented

Re: [R] Contingency Table Analysis - specific cell to specific cell comparisons

. Consult your instructor for further help. Say, to compare pairs of cells: along a column: 3 vs 31, 9 vs 31, 3 vs 9 along a row: 36 vs 9 or even across column and row: 36 vs 31, and 36 vs 3 David Winsemius, MD West Hartford, CT __ R-help@r-project.org

Re: [R] Error message: attempt to set rownames on object with no dimensions

don't use the as.factor I have no problems but I don't think this is the correct way to proceed... The other variables I use can have any value. -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch

Re: [R] Interactions in GAMMs

not changes since the book was published.) -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] [SPAM](Aktien) Re: Chi-square distribution probability density function:

On Jul 21, 2010, at 12:36 PM, Knut Krueger wrote: David Winsemius schrieb: On Jul 21, 2010, at 6:52 AM, Knut Krueger wrote: Hi to all I found an formular of an ** ***p-Value Calculator for the Chi-Square test* *http://www.danielsoper.com/statcalc/calc11.aspx* *with the formula* http

Re: [R] Chi-square distribution probability density function:

On Jul 21, 2010, at 1:18 PM, Knut Krueger wrote: David Winsemius schrieb: And exactly why did you think I offered ?Chisquare I was completely on the wrong way, and tried to find a solution with the formula instead to substitute the formula. So I tried to implement pchisq

Re: [R] Command that is conditional upon file retrieval: is it possible?

or tryCatch which might be more specific. Andrew -- -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] post hoc test for lme using glht ?

in contrMat(table(mf[[nm]]), type = types[pm]) : less than two groups Thanks, Rob -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

Re: [R] xtable

*** David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained

Re: [R] Sweave special characters problem

, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self

Re: [R] lme4 on Mac OS X

, Rhelp08, Rhelp02, functions ), matchesPerPage = 100, ...) RSiteSearch(string=string, restrict = restrict, matchesPerPage = matchesPerPage, ...) David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing

Re: [R] t.test in for loop

in the vars- vector (running t.test(data1) by hand works). How can I remove the quotes before performing the test? a) don't quote them b) put them in a list c) print() the results. -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing

Re: [R] Bar Plot Bars Bleed off Plotting Area

-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help

Re: [R] , how to express bar(zeta) in main title in boxplot

, nor is it in the argument list for boxplot.default or its bxp function. The documentation for the ... argument does not suggest, to me anyway, that main would passed on to other graphical functions. -- David Winsemius, MD West Hartford, CT __ R

Re: [R] , how to express bar(zeta) in main title in boxplot

On Jul 22, 2010, at 5:01 PM, Peter Ehlers wrote: On 2010-07-22 14:40, David Winsemius wrote: On Jul 22, 2010, at 4:24 PM, Peter Ehlers wrote: On 2010-07-22 11:44, Marcus Liu wrote: Hi everyone, I am plotting a boxplot with main title as main = bquote(paste(.(ts.ind[s]), : , bar(zeta

Re: [R] function return

commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] function return

On Jul 22, 2010, at 5:34 PM, David Winsemius wrote: On Jul 22, 2010, at 5:27 PM, Daniel Hocking wrote: I am sorry if this question is vague or uninformed. I am just learning R and struggling. I am using the book Hierarchical Modeling and Inference in Ecology and they provide examples of R

Re: [R] Filtering in R

: The dataframe is id salary 100500 101600 102700 103800 how can i generate a subsets if salary600? ?subset subset(DF, salary 600) -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

Re: [R] Midpoint between coordinates

-4.698 59.197 -4.5480 57.7830 -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide

Re: [R] Midpoint between coordinates

, second is latitude) or a SpatialPoints* object p2 As above n The requested number of points on the Great Circle --- There is also a geospatial Task View and SIG mailing list that is very active. -- David Thank you Mafalda On 23 July 2010 13:30, David Winsemius dwinsem...@comcast.net wrote

Re: [R] Figures in Latex

place to ask such a question. But I do not know who to ask, http://lmgtfy.com/?q=latex+users+group so just try my luck here. Thank you in advance. Hannah [[alternative HTML version deleted]] -- David Winsemius, MD West Hartford, CT

Re: [R] randomness using runif

-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

Re: [R] converting a time to nearest half-hour

were 11:59:ss, I want to convert to 12:00:00). How to achieve this? Couldn't you just coerce to numeric, divide by 60(sec)*30(half-hour minutes), round to integer, multiply by 60*30, coerce to POSIXct? _ David Winsemius, MD West Hartford, CT __ R

Re: [R] decimal seperator

, cex=1.5) text(19,4.2, labels=newdata[19], col=black, cex=1.5);text(20,3.55, labels=newdata[20], col=black, cex=1.5) David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] converting a time to nearest half-hour

On Jul 23, 2010, at 11:35 AM, David Winsemius wrote: On Jul 23, 2010, at 11:20 AM, murali.me...@avivainvestors.com murali.me...@avivainvestors.com wrote: Hi folks, I've got a POSIXct datum as follows: Sys.time() [1] 2010-07-23 11:29:59 BST I want to convert this to the nearest half

Re: [R] converting a time to nearest half-hour

hth Stephen Sefick On Fri, Jul 23, 2010 at 11:00 AM, David Winsemius dwinsem...@comcast.net wrote: On Jul 23, 2010, at 11:35 AM, David Winsemius wrote: On Jul 23, 2010, at 11:20 AM, murali.me...@avivainvestors.com murali.me...@avivainvestors.com wrote: Hi folks, I've got a POSIXct

Re: [R] Error produced by read.zoo: bad entries

But the structure of my OrigData is exactly the same as of monthly. OK - OrigData always has a few NAs in Value coming first - but that's consistent for all brands. Any idea what might be wrong? Thanks a lot! Just in case -attaching the actual file. No. Not attached. -- David Winsemius, MD

Re: [R] Error produced by read.zoo: bad entries

. But neither split = 1 nor split =2 work. And split =Brand does not work either. Why? D. On Fri, Jul 23, 2010 at 12:52 PM, David Winsemius dwinsem...@comcast.net wrote: ?read.zoo You didn't specify the index column correctly. On Jul 23, 2010, at 12:36 PM, Dimitri Liakhovitski wrote: Hello! I

Re: [R] Error produced by read.zoo: bad entries

what was expected, I am not in a position to know if this is success. And under warnings() it says: 1: In zoo(rval4[[i]], ix[[i]]) : some methods for “zoo” objects do not work if the index entries in ‘order.by’ are not unique On Fri, Jul 23, 2010 at 1:13 PM, David Winsemius dwinsem

Re: [R] Error produced by read.zoo: bad entries

, 98L, 95L, 97L, NA, 98L, 97L, 98L, NA, : series cannot be merged with non-unique index entries in a series In addition: There were 11 warnings (use warnings() to see them) On Fri, Jul 23, 2010 at 1:41 PM, David Winsemius dwinsem...@comcast.net wrote: On Jul 23, 2010, at 1:39 PM, Dimitri

Re: [R] Convert Row Names to data.frame column

with doubled \ but the second . is a regex for any character.) gsub(\\..*$,, c(Floodplain1.Tree, Floodplain1.Shrub)) [1] Floodplain1 Floodplain1 Use the same method with assignment to the column: ZoneCover.df$ID - gsub(\\..*$,, ZoneCover.df$ID) -- David Winsemius, MD West Hartford, CT

Re: [R] How to import simple java/mathematica expression to R

to minimize (with some constrains) here: http://analytic-products4you.com/target.txt Is it possible to import it as a function to minimize? What is the problem that you are trying to solve? We solve a problem for parameters estimation with ties. -- David Winsemius, MD West Hartford, CT

Re: [R] Trouble retrieving the second largest value from each row of a data.frame

On Jul 23, 2010, at 9:20 PM, mpw...@illinois.edu wrote: I have a data frame with a couple million lines and want to retrieve the largest and second largest values in each row, along with the label of the column these values are in. For example row 1 strongest=-11072 secondstrongest=-11707

Re: [R] Trouble retrieving the second largest value from each row of a data.frame

On Jul 24, 2010, at 4:54 PM, mpw...@illinois.edu wrote: THANKS, but I have one issue and one question. For some reason the secondstrongest value for row 3 and 6 are incorrect (they are the strongest) the remaining 10 are correct?? In my run of Wiley's code I instead get identical values

Re: [R] Trouble retrieving the second largest value from each row of a data.frame

On Jul 24, 2010, at 8:09 PM, David Winsemius wrote: On Jul 24, 2010, at 4:54 PM, mpw...@illinois.edu wrote: THANKS, but I have one issue and one question. For some reason the secondstrongest value for row 3 and 6 are incorrect (they are the strongest) the remaining 10 are correct

Re: [R] Trouble retrieving the second largest value from each row of a data.frame

On Jul 24, 2010, at 9:27 PM, David Winsemius wrote: On Jul 24, 2010, at 8:09 PM, David Winsemius wrote: On Jul 24, 2010, at 4:54 PM, mpw...@illinois.edu wrote: THANKS, but I have one issue and one question. For some reason the secondstrongest value for row 3 and 6 are incorrect

Re: [R] R equivalent of SAS proc freq

of this theory with RSeek pointed to the CrossTable function in gmodels.) Try also with... ?RsiteSearch RSiteSearch() Or with ... require(sos) ???proc freq sas -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https

Re: [R] 3d topographic map

On Jul 25, 2010, at 6:31 PM, sh...@ucar.edu wrote: Hi All- I would like to create a 3d topographic map using lat/lon and z(height). I have been scouring the R help pages and have not located the package I am looking for. Does anyone have a suggestion of package that will work for

Re: [R] Marginal effects from interaction regression model

-- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible

Re: [R] Marginal effects from interaction regression model

sure someone here can finish the job of showing you how setting up a call to predict will get you to that goal. -- David. Thanks again David Winsemius wrote: On Jul 25, 2010, at 10:24 PM, Guillem R. wrote: Dear all, I'd like to plot the marginal effect of a variable

Re: [R] Plot of a subset of a data.frame()

numbers. is there a way out? Best regards, /Steffen -- David Winsemius, MD Heritage Laboratories West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

<    1   2   3   4   5   6   7   8   9   10   >