Re: [R] Function try and Results of a program

2010-09-04 Thread David Winsemius
-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] Function try and Results of a program

2010-09-04 Thread David Winsemius
= \n,s, \n) out-list(s=s,d=d); sink() return(out) } View this message in context: http://r.789695.n4.nabble.com/Function-try-and-Results-of-a-program-tp2526621p2526822.html Sent from the R help mailing list archive at Nabble.com. -- David Winsemius, MD West Hartford, CT

Re: [R] What solve() does?

2010-09-04 Thread David Winsemius
-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 and provide commented, minimal

Re: [R] non-zero exit status error when install GenomeGraphs

2010-09-04 Thread David Winsemius
, 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-contained

Re: [R] Please explain do.call in this context, or critique to stack this list faster

2010-09-04 Thread David Winsemius
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-project.org/posting

Re: [R] Linear Logistic Regression - Understanding the output (and possibly the test to use!)

2010-09-04 Thread David Winsemius
of the variables has the influence on the outcome in the anorexia data? Please don't shout!! happy to be pointed to a reference but would prefer one in common english not some stats mumbo jumbo! Calum -- David Winsemius, MD West Hartford, CT

Re: [R] How can I fixe convergence=1 in optim

2010-09-04 Thread David Winsemius
at the answer? You are the one who know what f is and you are the one who has the option of increasing maxit. If the question is how to increase maxit, then the answer is perhaps as easy as: ?optim -- David Winsemius, MD West Hartford, CT __ R-help@r

Re: [R] Linear Logistic Regression - Understanding the output (and possibly the test to use!)

2010-09-05 Thread David Winsemius
On Sep 5, 2010, at 6:06 AM, st...@wittongilbert.free-online.co.uk wrote: David Winsemius wrote: 1. is glm the right thing to use before I waste my time Yes, but if your outcome variable is binomial then the family argument should be binomial. (And if you thought it should

Re: [R] converting string vector to integer/numeric vector

2010-09-05 Thread David Winsemius
to decimal ASCII [1] 49 50 51 32 116 104 105 115 32 105 115 32 97 32 115 116 114 105 110 103 -- 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] converting string vector to integer/numeric vector

2010-09-05 Thread David Winsemius
possible type identifiers? We need an example that has enough complexity to allow testing. -- David. On Sun, Sep 5, 2010 at 6:44 PM, David Winsemius dwinsem...@comcast.net wrote: On Sep 5, 2010, at 8:48 AM, rajesh j wrote: Hi, Is it possible to convert a string vector to integer

Re: [R] converting string vector to integer/numeric vector

2010-09-05 Thread David Winsemius
. The number of items in the vector is unknown. here's an example, a list has vectors INT 2 3 4 NUM 2.37 4.56 On Sun, Sep 5, 2010 at 6:56 PM, David Winsemius dwinsem...@comcast.net wrote: On Sep 5, 2010, at 9:22 AM, rajesh j wrote: for e.g., I get the following as a string vector int

Re: [R] converting string vector to integer/numeric vector

2010-09-05 Thread David Winsemius
-list(INT=c(1,2,3),NUM=c(2.34,4.56,6.78)) cc - lapply(cc, as.numeric) cc $INT [1] 1 2 3 $NUM [1] 2.34 4.56 6.78 -- David. On Sun, Sep 5, 2010 at 7:57 PM, David Winsemius dwinsem...@comcast.net wrote: So there is one item per line and the task is to recognize the strings INT and NUM

Re: [R] dirichlet models

2010-09-05 Thread David Winsemius
://finzi.psych.upenn.edu/R/library/VGAM/html/dirichlet.html (I'm looking for something like this: http://repec.org/bocode/d/dirifit.html , that allows for both dependent variables summing to 1 predictive variables of any sort.) Don -- Donald Braman -- David Winsemius, MD West Hartford, CT

Re: [R] Greek symbols (again but more complicated)

2010-09-05 Thread David Winsemius
-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] appending to a list

2010-09-05 Thread David Winsemius
On Sep 5, 2010, at 1:21 PM, Aks Ism wrote: Hi, I've looked at previous discussions and did not get anything. I want to be able to append to a list in a loop. Is this possible? Of course: ?c ?[[ -- David Winsemius, MD West Hartford, CT __ R

Re: [R] Need Help .RData in Mac OS X, Please

2010-09-05 Thread David Winsemius
have entered: defaults write com.apple.Finder AppleShowAllFiles YES . and then pt-click-hold on Dock-Finder-icon, choose relaunch Or you could open a Terminal window which will by default open in / Users/ctu/ and type: rm .RData -- David Winsemius, MD West Hartford, CT

Re: [R] appending to a list

2010-09-05 Thread David Winsemius
. Perhaps you mean something like this: ll - list(); for (i in 1:20) { my.element - scan(); ll - c(ll, my.element); if (is.na(my.element)){return(ll);break}} -- David. On Mon, Sep 6, 2010 at 1:01 AM, David Winsemius dwinsem...@comcast.net wrote: On Sep 5, 2010, at 1:21 PM

Re: [R] extracting x,y coordinates from a contour plot

2010-09-05 Thread David Winsemius
330 ... $ :List of 3 ..$ level: num 160 ..$ x: num [1:31] 270 263 262 260 260 ... ..$ y: num [1:31] 310 320 330 340 350 ... -- David. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

Re: [R] Correct coefficients from treatment contrasts?

2010-09-06 Thread David Winsemius
Ph: 1.250.359.5905 [[alternative HTML version deleted]] 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] Aggregating the matrices

2010-09-06 Thread David Winsemius
. Sergey __ 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

Re: [R] mac: lib/gtk.pkg

2010-09-06 Thread David Winsemius
that package, sometimes reads rhelp but I'm not sure on what schedule. You might see if he makes his email address available on those pages. -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

Re: [R] Aggregate certain rows in a matrix

2010-09-06 Thread David Winsemius
://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 mailing list

Re: [R] WriteXLS problem

2010-09-06 Thread David Winsemius
. __ 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

Re: [R] poisson distribution

2010-09-06 Thread David Winsemius
On Sep 6, 2010, at 1:13 PM, tamas barjak wrote: Hello! I need some help. How I know it to draw the formula of the poisson distribution? expr-expression(P(xi == k) == frac(lambda^k, factorial(k))*e^- lambda) --- not good ?plotmath (Do not see factorial as a plotmath function Try:

Re: [R] Failure to aggregate

2010-09-06 Thread David Winsemius
On Sep 6, 2010, at 12:15 PM, Dimitri Shvorob wrote: I have a (very big - 1.5 rows) dataframe with a (POSIXt POSIXlt) column h (hour). Surprisingly, I cannot calculate a simple aggregate over the dataframe. n.h1 = sqldf(select distinct h, count(*) from x group by h) Error in

Re: [R] how to change the xlab name?

2010-09-06 Thread David Winsemius
On Sep 6, 2010, at 2:07 PM, tooblue wrote: I simply put, plot(density(), main=, + xlab = XXX), it says that I have an unexpected = in it. It may be a case of a confused parser. You have an extraneous + in there: = rnorm(100) plot(density(), main=, xlab = XXX)

Re: [R] two questions

2010-09-06 Thread David Winsemius
talking about? Thanks Jason Dr. Iasonas Lamprianou -- 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] WriteXLS problem

2010-09-06 Thread David Winsemius
-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] how to combine several subsets?

2010-09-06 Thread David Winsemius
since you assigned a result to those names. (Greater clarity would occur if you offered at least str(NEVER) -- 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] Prediction and confidence intervals from predict.drc

2010-09-06 Thread David Winsemius
On Sep 6, 2010, at 7:54 PM, Brant Inman wrote: R-helpers, I am using the package drc to fit a 4 parameter logistic model. When I use the predict function to get prediction on a new dataset, I am not getting the requested confidence or prediction intervals. Any idea what is going on?

Re: [R] Saving fits (glm, nls) without data

2010-09-07 Thread David Winsemius
. -- 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] boundary correction - univariate kernel density estimation

2010-09-07 Thread David Winsemius
density estimation function (e.g. lambdahat in the spatialkernel package) with boundary corrections to allow univariate density estimation? When this question has been (multiply) posed in the past, the suggested answer has been to use package logspline. -- David Winsemius, MD West Hartford

Re: [R] a^c(1:3)

2010-09-07 Thread David Winsemius
, 4) [,1] [,2] [,3] [1,]1 81 125 [2,]8 16 1296 BTW, I have a 64bit R version (2.11) for Linux. Any advice would be appreciated. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

Re: [R] remove accents in strings

2010-09-07 Thread David Winsemius
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] remove accents in strings

2010-09-07 Thread David Winsemius
, I'll take this moment to remind anyone interested that R still has trouble with embedded zeros in character strings. I may be abusing terminology, but I think that makes R 8-bit dirty. -Matt On Tue, 2010-09-07 at 14:01 -0400, David Winsemius wrote: On Sep 7, 2010, at 1:35 PM, Matt Shotwell wrote

Re: [R] Saving fits (glm, nls) without data

2010-09-07 Thread David Winsemius
On Sep 7, 2010, at 2:53 PM, Johann Hibschman wrote: David Winsemius dwinsem...@comcast.net writes: On Sep 7, 2010, at 11:02 AM, Johann Hibschman wrote: Even so, I would prefer to only save the coefficients Have you read through the Value section of glm's help page? ...and ?coef I have

Re: [R] Saving fits (glm, nls) without data

2010-09-07 Thread David Winsemius
On Sep 7, 2010, at 3:16 PM, David Winsemius wrote: On Sep 7, 2010, at 2:53 PM, Johann Hibschman wrote: David Winsemius dwinsem...@comcast.net writes: On Sep 7, 2010, at 11:02 AM, Johann Hibschman wrote: Even so, I would prefer to only save the coefficients Have you read through

Re: [R] Help with decimal points

2010-09-07 Thread David Winsemius
(MatchedValues[Value2,Value]) [1] 4420 -- 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] change the for loops with lapply

2010-09-07 Thread David Winsemius
parameters some time. * Can anyone help me this this? Thanks so much! -- Sincerely, Changbin 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

Re: [R] problem with max in a function

2010-09-07 Thread David Winsemius
-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] R time series analysis

2010-09-07 Thread David Winsemius
series when se.fit = TRUE. -- View this message in context: http://r.789695.n4.nabble.com/R-time-series-analysis-tp2527513p2530595.html -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r

Re: [R] multiple graphs

2010-09-07 Thread David Winsemius
. -- 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] R time series analysis

2010-09-07 Thread David Winsemius
On Sep 7, 2010, at 9:33 PM, David Winsemius wrote: On Sep 7, 2010, at 7:51 PM, lord12 wrote: For each arima model, can you output an associated confidence interval for the predicted value at each time point? ?arima0 arima0 will return ... a list with components pred

Re: [R] problem with max in a function

2010-09-07 Thread David Winsemius
])]) s - with(measure_bkf_not_zero[bkf_min:bkf_max,], approxfun(bankfull_depths_m, measurment_num), ties=mean) int_bkf - s(0) On Tue, Sep 7, 2010 at 8:28 PM, David Winsemius dwinsem...@comcast.net wrote: On Sep 7, 2010, at 9:06 PM, stephen sefick wrote: s - 1.00 max(s) sprintf(%.2f, max(s

Re: [R] problem with max in a function

2010-09-07 Thread David Winsemius
[- which.max(vector)] ). So as Jim Holtman's tag line says: what problem are you trying to solve? Again, I am sorry for being vague. I get wrapped up in a problem and forget that I need to communicate. kindest regards, Stephen On Tue, Sep 7, 2010 at 8:48 PM, David Winsemius dwinsem

Re: [R] Regression using mapply?

2010-09-08 Thread David Winsemius
of the for-loop using mapply (or some other function)? -- 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] How to change font size in plot() function

2010-09-08 Thread David Winsemius
-font-size-in-plot-function-tp2531127p2531161.html Sent from the R help mailing list archive at Nabble.com. 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] Checking if the distribution follow a power law

2010-09-08 Thread David Winsemius
.html and I got results like this: Profiling... 2.5 % 97.5 % 2.393297 2.412650 What do these suggest? please tell me about this if someone knows. -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https

Re: [R] try-error can not be test. Why?

2010-09-08 Thread David Winsemius
-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] Replace NAs in one column with data from another column

2010-09-08 Thread David Winsemius
])} } Works - but takes too long time. I would appreciate alternative solutions. Best regards, Jakob -- 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] subbing a string vector for another string vector

2010-09-08 Thread David Winsemius
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-project.org/posting

Re: [R] Newbie cross tabulation issue

2010-09-08 Thread David Winsemius
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] coxph and ordinal variables?

2010-09-08 Thread David Winsemius
]] __ 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] Newbie cross tabulation issue

2010-09-08 Thread David Winsemius
for your help. -- Jonathan. [[alternative HTML version deleted]] -- 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] problem with outer

2010-09-08 Thread David Winsemius
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-project.org/posting

Re: [R] on error execute:

2010-09-08 Thread David Winsemius
messages and continue executing? ?try Sorry if this has been addressed already, but, even aided by the new awesome google i haven't been able to find it. -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https

Re: [R] try-error can not be test. Why?

2010-09-08 Thread David Winsemius
On Sep 8, 2010, at 11:46 PM, Philippe Grosjean wrote: On 08/09/10 19:25, David Winsemius wrote: On Sep 8, 2010, at 1:18 PM, telm8 wrote: Hi, I am having some strange problem with detecting try-error. From what I have read so far the following statement: try( log(a) ) == try-error

Re: [R] Determine Bounds of Current Graph

2010-09-09 Thread David Winsemius
('state','Indiana') map.axes() ?par bounds - par(usr) bounds [1] -88.12964 -84.77184 37.74583 41.82082 ??Something that lets me know the y-axis is from ~38 to ~42 and store this information into a vector Is there some way to query what the bounds of the current graph are? Thanks! David

Re: [R] Calculating with tolerances (error propagation)

2010-09-09 Thread David Winsemius
/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] Calculating with tolerances (error propagation)

2010-09-09 Thread David Winsemius
/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] Newbie cross tabulation issue

2010-09-09 Thread David Winsemius
Hmisc::summary.formula -- 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] Newbie cross tabulation issue

2010-09-09 Thread David Winsemius
On Sep 9, 2010, at 11:20 AM, David Winsemius wrote: On Sep 8, 2010, at 7:32 PM, Jonathan Finlay wrote: Thanks David, gmodels::Crosstable partially work because can show only 1 x 1 tablen CrossTable(x,y,...) I need something how can process at less 1 variable in X an 10 in Y. A further

Re: [R] Newbie cross tabulation issue

2010-09-09 Thread David Winsemius
that attitude and there are now quite a few worked examples of common task in which SAS, SPSS, Stata and R are applied using the same data. http://www.ats.ucla.edu/stat/dae/ -- David. HTH, Marc Schwartz David Winsemius, MD West Hartford, CT __ R-help

Re: [R] sequeeze a data frame

2010-09-09 Thread David Winsemius
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] URL error when trying to use help function in R [Sec: UNOFFICIAL]

2010-09-09 Thread David Winsemius
, 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] survfit question

2010-09-09 Thread David Winsemius
/+ groups -- 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] Calculating with tolerances (error propagation)

2010-09-10 Thread David Winsemius
On Sep 9, 2010, at 10:57 AM, David Winsemius wrote: On Sep 9, 2010, at 6:50 AM, Jan private wrote: Hello Bernardo, - If I understood your problem this script solve your problem: q-0.15 + c(-.1,0,.1) h-10 + c(-.1,0,.1) 5*q*h [1] 2.475 7.500 12.625 - OK, this solves

Re: [R] Data.frames : difference between x$a and x[, a] ? - How set new values on x$a with a as variable ?

2010-09-10 Thread David Winsemius
/ __ 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] convert 1, 10, and 100 to 0001, 0010, 0100 etc.

2010-09-10 Thread David Winsemius
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/listinfo

Re: [R] Solver in R

2010-09-10 Thread David Winsemius
. Is there any package can do it? If you minimize R0 - sum(r * Xt ), you should get your answer. -- David. Thanks. Peter David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] Solver in R

2010-09-10 Thread David Winsemius
On Sep 10, 2010, at 5:44 PM, David Winsemius wrote: On Sep 10, 2010, at 5:35 PM, Chien-Pang Chin wrote: Hi all: I'm looking for a package that similar to solver in MS. All I need is find a r to satisfy R0=sum( , where t are from 1 to n and Xt are come from another formula. The most

Re: [R] Where to find R-help options web page

2010-09-10 Thread David Winsemius
of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Where-to-find-R-help-options-web-page-tp2535123p2535148.html Sent from the R help mailing list archive at Nabble.com. -- David Winsemius, MD West Hartford, CT

Re: [R] Formatting of time strings

2010-09-10 Thread David Winsemius
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] Greek letter included in a character vector

2010-09-10 Thread David Winsemius
(3,4,6), text.col = green4, lty = c(2, -1, 1), pch = c(-1, 3, 4), merge = TRUE, bg = 'gray90') -- 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] convert 1, 10, and 100 to 0001, 0010, 0100 etc.

2010-09-10 Thread David Winsemius
-and-100-to-0001-0010-0100-etc-tp2535023p2535296.html Sent from the R help mailing list archive at Nabble.com. -- 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] Is there a bisection method in R?

2010-09-10 Thread David Winsemius
in advance. Gregory -- 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] convert 1, 10, and 100 to 0001, 0010, 0100 etc.

2010-09-10 Thread David Winsemius
On Sep 10, 2010, at 10:05 PM, Peng, C wrote: I mean to display 001,010, ..., as there are. In other words, whether there is a function, say func(), such that func(001,010) displays 001, 010. Not hard to construct one, but does not behave properly in the sub- unity decimal range. Not

[R] acceptable filetypes was Re: coxph and ordinal variables?

2010-09-10 Thread David Winsemius
On Sep 10, 2010, at 10:32 PM, Paul Johnson wrote: Hi, everybody On Wed, Sep 8, 2010 at 5:43 PM, Min-Han Tan minhan.scie...@gmail.com wrote: David said my R code text attachment got rejected by the mailing list. Pooh. I don't think that's nice. I don't see anything in the posting

Re: [R] Argument lib is missing

2010-09-11 Thread David Winsemius
at csdb.cn is 2.7.1 from May 2009. There have been several revisions since that date.) Warning message: In getDependencies(pkgs, dependencies, available, lib) : package ‘gplots’ is not available any suggestion will be appreciable. Thank you Khush David Winsemius, MD West Hartford, CT

Re: [R] Latex fonts in R graphics

2010-09-11 Thread David Winsemius
that you do not have Ghostscript installed properly. While you are at it, you should probably updated R as well. -- David Any help? -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r

Re: [R] for loop

2010-09-11 Thread David Winsemius
On Sep 11, 2010, at 8:39 AM, Peng, C wrote: or: k=0 for (i in 1:k) if(k0) print(i) Because of the way the : operator works, I would have tested k =1 k=0.5 for (i in 1:k) if (k0){print(i)} [1] 1 But Gabor's suggestion to use seq_len(k) is cleaner, anyway. -- David Winsemius, MD West

Re: [R] Argument lib is missing

2010-09-11 Thread David Winsemius
On Sep 11, 2010, at 8:02 AM, David Winsemius wrote: On Sep 11, 2010, at 6:29 AM, khush wrote: Dear all, I have installed R using yum install R-2.9. I am able to use R for general functions but when I installed some library say gplots I am getting the below error

Re: [R] bootstrapping

2010-09-11 Thread David Winsemius
the elements that were not included in the list? I realize i can do this with a for loop by going through each element and checking if it's in the list but I am wondering if there is a faster way. ?setdiff ?%in% -- David Winsemius, MD West Hartford, CT

Re: [R] approxfun-problems (yleft and yright ignored)

2010-09-11 Thread David Winsemius
survival_2.35-8 plyr_1.1 [25] lattice_0.18-8 loaded via a namespace (and not attached): [1] tools_2.11.1 -- David Winsemius I get completely different (and correct) results, by the way the *same* you have in the bug report you've submitted (https://bugs.r-project.org/bugzilla3/show_bug.cgi?id=14377

Re: [R] xlab with text and expression

2010-09-11 Thread David Winsemius
=expression(abc~x = 1)) I want proper sign of weak inequality not just '=' will appreciate! robert -- 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] confidence bands for a quasipoisson glm

2010-09-11 Thread David Winsemius
. Best regards Maik 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] scalable delimiters in plotmath

2010-09-11 Thread David Winsemius
a workaround, grid.text(expression(bgroup(,atop(x,y),))) Error in bgroup(, atop(x, y),) : invalid group delimiter Regards, baptiste sessionInfo() R version 2.11.1 (2010-05-31) x86_64-apple-darwin9.8.0 David Winsemius, MD West Hartford, CT __ R-help@r

Re: [R] create a '3D line plot'

2010-09-12 Thread David Winsemius
achieve this using R, Karl -- Karl Brand -- 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] reshape matrix entities to columns

2010-09-12 Thread David Winsemius
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-project.org/posting

[R] How to print matrices in standard format was ... Re: How to define new matrix based on an elementary row oper

2010-09-12 Thread David Winsemius
,] -2 -4-3 -20 -18 for( i in 1:nrow(A) ) { cat(sprintf(%4.0f, A[i, ]), paste( ,if( i==3 ){+}else{ }, , sep=),sprintf(%4.0f,B[i, ]), paste( ,if( i==3 ){=}else{ }, , sep=), sprintf(%4.0f, (A +B)[i, ]), \n )} -- David Winsemius, MD West Hartford, CT

Re: [R] scalable delimiters in plotmath

2010-09-12 Thread David Winsemius
understand these codes might correspond to extended ascii characters whose boundaries and positions we want to borrow. Then again, maybe it's something else entirely. Any hints? Best wishes, baptiste On 12 September 2010 03:27, David Winsemius dwinsem...@comcast.net wrote: On Sep 11, 2010

Re: [R] scalable delimiters in plotmath

2010-09-12 Thread David Winsemius
Regards, baptiste sessionInfo() R version 2.11.1 (2010-05-31) x86_64-apple-darwin9.8.0 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] reshape matrix entities to columns

2010-09-12 Thread David Winsemius
On Sep 12, 2010, at 3:34 PM, Dennis Murphy wrote: Hi: Natasha said: I changed it so i hope it will look better now the matrix is like this: AgeNo. Age No. AgeNo. Center1 52 8 7 Center210 720 9 4

Re: [R] How to print matrices in standard format was ... Re: How to define new matrix based on an elementary row oper

2010-09-12 Thread David Winsemius
On Sep 12, 2010, at 12:24 PM, David Winsemius wrote: On Sep 12, 2010, at 11:27 AM, Cuckovic Paik wrote: I appreciate all you help. This is only for instructional purpose: A = matrix(c(0,1,1,-2,-3,1,2,-1,0,2,2,4,1,-3,-2,1,-4,-7,-1,-19), ncol=5, byrow=T) B = matrix(sample(c(0,1,1,-2

Re: [R] How to do a trig regression

2010-09-12 Thread David Winsemius
On Sep 12, 2010, at 10:23 PM, Aaditya Nanduri wrote: Hello All, I cant seem to do a trig regression in R. The equation is as follows : y = a+b*(sin((2*pi*x/360) - c))^2 a, b, c are coefs that I want. y, x are input vectors. The equation I put into R: lm(y ~ sin(2*pi*x/360)^2) This equation

Re: [R] Overlay of two graphs of different axes

2010-09-13 Thread David Winsemius
worked examples in the archives, as well as canned solutions in widely used packages. -- 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] Create a time-series from cross-sectional data that has each year as a separate column

2010-09-13 Thread David Winsemius
, 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-contained

Re: [R] Saveing plot to multiple locations

2010-09-13 Thread David Winsemius
at Nabble.com. __ 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

Re: [R] lattice: Set x-axis in italics only

2010-09-13 Thread David Winsemius
: ..., scales=list(x=list(font=3)), In this case (and I suppose many others) lattice is much more flexible than base. I initially tried what appeared to be valid strategies with boxplot and ended up tied in knots. Then I noticed your subject line and it was trivial. David Winsemius, MD West

Re: [R] Axis break with gap.plot()

2010-09-13 Thread David Winsemius
, 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

<    3   4   5   6   7   8   9   10   11   12   >