Re: [R] character set problem

2019-02-19 Thread Bert Gunter
Oops, you already seem to have done this! My bad. -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Tue, Feb 19, 2019 at 11:36 A

Re: [R] help getting a research project started on regulations.gov

2019-02-19 Thread Bert Gunter
Please search yourself first! "scrape JSON from web" at the rseek.org site produced what appeared to be several relevant hits, especially this CRAN task view: https://cran.r-project.org/web/views/WebTechnologies.html Bert Gunter "The trouble with having an open mind is that peo

Re: [R] R Software

2019-02-21 Thread Bert Gunter
Please stop these silly posts. R is open source software, and its open source licensing requirements are explained on its website and referenced links. As stated there, it comes with NO guarantees. The R Foundation is *not* a company. Bert Gunter "The trouble with having an open mind is

Re: [R] Cannot reproduce tutorial results

2019-02-22 Thread Bert Gunter
This is a plain text list and your html post below is pretty mangled and difficult to read. If you re-post in plain text, you are more likely to get a response. Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (ak

Re: [R] Newton-RaphsonMethod

2019-02-24 Thread Bert Gunter
This list has a no homework policy. Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Sun, Feb 24, 2019 at 9:22 AM malika yassa via R-help &l

Re: [R] How to change the number of bins?

2019-03-10 Thread Bert Gunter
You are asking the wrong question. The right question is, "why are so many values missing?" Is it because they were censored, not reported for some reason, due to instrument failure,...? Until you answer that question, any analysis you do is garbage. I strongly recommend you consult a competent

Re: [R] Sorting vector based on pairs of comparisons

2019-03-14 Thread Bert Gunter
This cannot be done unless transitivity is guaranteed. Is it? S L a b b c c a Bert On Thu, Mar 14, 2019, 4:30 AM Pedro Conte de Barros wrote: > Dear All, > > This should be a quite established algorithm, but I have been searching > for a couple days already without finding any satisfact

Re: [R] Sorting vector based on pairs of comparisons

2019-03-15 Thread Bert Gunter
If I understand correctly, the answer is a topological sort. Here is an explanation https://davidurbina.blog/on-partial-order-total-order-and-the-topological-sort/ This was found by a simple web search on "Convert partial ordering to total ordering" Btw. Please use search engines before posting

Re: [R] plot.xmean.ordinaly vs plot() in package "rms"

2019-03-20 Thread Bert Gunter
Please study the documentation to which you were referred. This list is not appropriate for comprehensive tutorials, which is what you need, although all help is of course in some sense a tutorial. Bert Gunter "The trouble with having an open mind is that people keep coming along and sti

Re: [R] System("source activate condaenv")

2019-03-20 Thread Bert Gunter
R is *not* RStudio. Please go to the RStudio site, not here, for help with that software. Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On

Re: [R] Solving Function using Conjugate Gradient

2019-03-20 Thread Bert Gunter
This list has a no homework policy. Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Wed, Mar 20, 2019 at 8:54 AM smart hendsome v

Re: [R] Solving Function using Conjugate Gradient

2019-03-20 Thread Bert Gunter
ttps://go.onelink.me/107872968?pid=InProduct&c=Global_Internal_YGrowth_AndroidEmailSig__AndroidUsers&af_wl=ym&af_sub1=Internal&af_sub2=Global_YGrowth&af_sub3=EmailSignature> > > On Thu, 21 Mar 2019 at 1:34 a.m., Bert Gunter > wrote: > This list has a no homework policy. > > Bert Gunter > > &quo

Re: [R] creating a dataframe with full_join and looping over a list of lists

2019-03-21 Thread Bert Gunter
received. Please read the ?data.frame or a tutorial on data frames. All columns *must* have the same length. So what it comes down to is probably how you want to fill with NA's. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and stick

Re: [R] loop through columns in a data frame

2019-03-25 Thread Bert Gunter
"Does anyone know how to use loop (or other methods) to create new columns? In SAS, I can use array to get it done. But I don't know how to do it in R." Yup. Practically all users of R know how, as this is entirely elementary. You will too if you make the effort to go through a basic R tutorial, o

Re: [R] Monte Carlo simulation for ratio and its CI

2019-03-25 Thread Bert Gunter
alculated, but that is another off topic discussion for which stats.stackexchange.com is a more appropriate venue. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bl

Re: [R] Substitution in expressions

2019-03-26 Thread Bert Gunter
I believe you're going about this the wrong way. You seem to want mathematical expressions. Fot this, see ?plotmath. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in hi

Re: [R] Substitution in expressions

2019-03-26 Thread Bert Gunter
Perhaps something like this (apologies if beating a dead horse): plot(NA,NA, xlim = c(-1,5),ylim = c(-1,5), xlab = "", ylab = "") for(i in 1:3) text(i,i,labels =bquote(2^.(i))) Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming

Re: [R] [FORGED] Re: Quantile Density Contours

2019-03-27 Thread Bert Gunter
You are missing a crucial point. The reals are well ordered; higher dimensions are not. Therefore 2d quantile contours are not unique. Of course assuming I understand your query correctly. Bert On Wed, Mar 27, 2019, 13:55 Bernard McGarvey wrote: > If I understand correctly the ContourLines fu

Re: [R] converting a character string to an object name

2019-03-29 Thread Bert Gunter
I think you want ?assign Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Fri, Mar 29, 2019 at 10:39 AM Assa Yeroslaviz wrote: > I am

Re: [R] Fwd: Potential Issue with lm.influence

2019-04-02 Thread Bert Gunter
Nothing was attached. The r-help server strips most attachments. Include your code inline. Also note that > 0/0 [1] NaN so maybe something like that occurs in the course of your calculations. But that's just a guess, so feel free to disregard. Bert Gunter "The trouble with h

Re: [R] Fwd: Potential Issue with lm.influence

2019-04-02 Thread Bert Gunter
I told you already: **Include code inline ** See ?dput for how to include a text version of objects, such as data frames, inline. Otherwise, I believe .txt text files are not stripped if you insist on *attaching* data or code. Others may have better advice. Bert Gunter "The trouble

Re: [R] Fwd: Potential Issue with lm.influence

2019-04-02 Thread Bert Gunter
Also, I suggest you read ?influence which may explain the source of your NaN's . Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Tu

Re: [R] Potential Issue with lm.influence

2019-04-03 Thread Bert Gunter
Second! Bert Gunter On Wed, Apr 3, 2019 at 9:35 AM Richard M. Heiberger wrote: > fortune nomination. > > > The lesson to me here is that if you fit a sufficiently unreasonable > model to data, the computations may break down. > > On Wed, Apr 3, 2019 at 10:18 AM Fox, John

Re: [R] function predict

2019-04-03 Thread Bert Gunter
This list has *no homework* policy. I would assume that the purpose of your "project" is for you to learn how to deal with exactly the sorts of issues you describe. (But you might get lucky with a response anyway). Bert Gunter "The trouble with having an open mind is that peo

Re: [R] R SIG mailing lists

2019-04-04 Thread Bert Gunter
lp/support/ "Healtrhcare" and "Biology" are far too vague and all-encompassing, imo. This might be of some use to you, however: https://cran.r-project.org/web/views/ -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into

Re: [R] Alternative to loops

2019-04-06 Thread Bert Gunter
I skipped pre-populating MyDF$C as unnecessary: > MyDF <- data.frame(A=c(1,2,3,4,5),B=c("aa ab ac","bb bc bd","cc cf","dd","ee"), + stringsAsFactors = FALSE) ## I think this does what you want: > choices<- sapply(MyDF$B, strsplit, split = " +") > nm <- names(MyList) > MyDF$C <- nm[sapply(choices

Re: [R] Unable to read csv files with comma in values

2019-04-07 Thread Bert Gunter
we need a loop: for(i in wh) txt[wh[i]] <- sub("\\[.+\\]",fixup[i],txt[wh[i]]) ## replace original bracketed text with fixed up bracketed text > txt [1] "Sam, [HadoopAnalyst-DBA-Developer], R46443 " [2] "Jan, DBA, R101" [3] "Mary, [HadoopAnalyst-DBA-Developer], t14"

Re: [R] Unable to read csv files with comma in values

2019-04-07 Thread Bert Gunter
quot;\\[|\\]") txt[wh] <- sapply(spl, function(y) paste0(y[1], gsub(" *, *","-", y[2]), y[-(1:2)])) > txt [1] "Sam, HadoopAnalyst-DBA-Developer, R46443 " [2] "Jan, DBA, R101" [3] "Mary, Stats-Designer-R, t14" Bert Gunter "The

Re: [R] Unable to read csv files with comma in values

2019-04-07 Thread Bert Gunter
","\\1",txt[wh]) ## before "[" txt2 <- gsub(" *, *","-",sub(".+(\\[.+\\]).+","\\1",txt[wh])) ## bracketed part txt3 <- sub(".*\\], *(.+?) *$","\\1",txt[wh]) ## after "]" txt[wh]<- paste(txt1, txt2, txt3, sep

Re: [R] How can I solve this prediction problem?

2019-04-09 Thread Bert Gunter
This looks like homework, and there is a no homework policy on this list. Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Tue, Apr 9, 2019

Re: [R] Define pch and color based on two different columns

2019-04-09 Thread Bert Gunter
1. I am quite sure that whatever it is that you want to do can be done. Probably straightforwardly. The various R graphics systems are mature and extensive. 2. But I, for one, do not understand from your post what it is that you want to do. Nor does anyone else apparently. Cheers, Bert Gunter

Re: [R] Define pch and color based on two different columns

2019-04-09 Thread Bert Gunter
}) Cheers, Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Tue, Apr 9, 2019 at 9:53 PM Peter Langfelder wrote: > Sorry for being

Re: [R] Bifactor model and infit statistics?

2019-04-11 Thread Bert Gunter
ll as searching at rseek.org . Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Thu, Apr 11, 2019 at 2:23 PM kende jan via R-help wrote:

Re: [R] Are fitted.values available in pglm?

2019-04-12 Thread Bert Gunter
?fitted ?predict ## This is what one usually does, but I have not checked pglm. You also need to get friendly with ?str ... and probably also spend time with an R tutorial or two to become familiar with R modeling conventions. Bert Gunter "The trouble with having an open mind is that p

Re: [R] Are fitted.values available in pglm?

2019-04-12 Thread Bert Gunter
ould be using residuals.?pglm also > includes residuals in its list of elements. However, str(mymodel) does not > mention residuals. Does that mean it’s just not there? > > - Simon > > On Apr 12, 2019, at 10:44 AM, Bert Gunter wrote: > > ?fitted > ?predict > ##

Re: [R] Syntax differences between aov and lmer for 2-way repeated measures design using a mixed model

2019-04-12 Thread Bert Gunter
/or lme4 packages are usually better posted on the r-sig-mixed-models list. ... and if this is homework, this list has a no homework poilicy. Cheers, Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley B

Re: [R] create

2019-04-13 Thread Bert Gunter
If the NA's are really 0's, replace them with 0 before doing the calculation. (see ?is.na). If they are not 0's, think again about doing this as the results would probably mislead. Bert Gunter "The trouble with having an open mind is that people keep coming along and sti

Re: [R] create

2019-04-13 Thread Bert Gunter
t that's for you to determine. Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Sat, Apr 13, 2019 at 9:01 PM Val wrote: > Sorry for th

Re: [R] Why is it not possible to cut a tree returned by Agnes or Diana by height?

2019-04-14 Thread Bert Gunter
Inline. Bert Gunter On Sun, Apr 14, 2019 at 4:12 PM Leszek Nowina wrote: > > asdf = data.frame(x=c(1,2,3), y=c(4,5,6), z=c(7,8,9)) > > cutree(agnes(asdf), h=100) > Error in cutree(agnes(asdf), h = 100) : > the 'height' component of '

Re: [R] Limiting the scope of RNGkind/set.seed

2019-04-16 Thread Bert Gunter
} > f() [1] "Super-Duper" "Inversion" > RNGkind() [1] "Mersenne-Twister" "Inversion" Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed

Re: [R] 1st el of a list of vectors

2014-07-22 Thread Bert Gunter
u tell vapply the return type and lapply always returns a list. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Tue, Jul 22, 2014 at 1:14 AM, aru

Re: [R] Partition of sums of squares (ANOVA)

2014-07-22 Thread Bert Gunter
Wrong list! Try stats.stackexchange.com for statistics questions. This list is about R programming related issues. Also, note that HTML does not work and should not be used here. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Inform

Re: [R] Using apply() with functions I wrote

2014-07-24 Thread Bert Gunter
ummm R is case sensitive! "fun" != "FUN" (Have you gone through any R tutorials yet? If not, please do so before posting further). Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. A

Re: [R] How to modify the body of a function?

2014-07-27 Thread Bert Gunter
would probably be a better place to look for thoughts anyway. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Sun, Jul 27, 2014 at 11:37 AM,

Re: [R] analyzing qualitative data sets

2014-07-29 Thread Bert Gunter
program in R and have not already done so, please read "An Introduction to R" or R web tutorial of your choice before posting here further. Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge i

Re: [R] Dependency Injection & Inversion of Control for Data

2014-07-29 Thread Bert Gunter
say so. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Tue, Jul 29, 2014 at 9:57 AM, William Dunlap wrote: > R is a functional language

Re: [R] DATA SUMMARIZING and REPORTING

2014-07-30 Thread Bert Gunter
Is this homework? There is a no homework policy here. And stop posting in HTML --- plain text only-- and learn to use ?dput to post example data. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowled

Re: [R] mutually exclusive events

2014-08-02 Thread Bert Gunter
Homework? There is a no homework policy here. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Sat, Aug 2, 2014 at 11:47 AM, Don McKen

Re: [R] Count number of change in a specified time interval

2014-08-04 Thread Bert Gunter
Or ?rle Bert Sent from my iPhone -- please excuse typos. > On Aug 4, 2014, at 8:28 AM, jim holtman wrote: > > Try this, but I only get 2 changes for CB27A instead of you indicated 3: > >> require(data.table) >> x <- read.table(text = "CASE_ID YEAR_MTH ATT_1 > + CB26A201302 1 > +

Re: [R] Prediction intervals (i.e. not CI of the fit) for monotonic loess curve using bootstrapping

2014-08-12 Thread Bert Gunter
nsult references or your local statistician for help if needed. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Tue, Aug 12, 2014 at 8:20 AM, David W

Re: [R] generating a sequence of seconds

2014-08-12 Thread Bert Gunter
Marc: You just need to be more patient -- this is already happening: http://en.wikipedia.org/wiki/Tidal_acceleration Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not w

Re: [R] Operating on the value from row i and row i+1

2014-08-14 Thread Bert Gunter
Your query is a bit unclear, but I suspect ?plot and a **careful read** about types "s" and "S" therein would address your problem. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. A

Re: [R] Operating on the value from row i and row i+1

2014-08-14 Thread Bert Gunter
?findInterval -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Thu, Aug 14, 2014 at 9:39 AM, Jaiprasart, Pharavee (HSC) wrote: > Hi Bert,

Re: [R] Negative values on output

2014-08-19 Thread Bert Gunter
Please follow the posting guide and post in plain text, not HTML, so that helpeRs do have to try to decipher this mess. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not w

Re: [R] anova

2014-08-21 Thread Bert Gunter
, not here. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Thu, Aug 21, 2014 at 3:13 AM, Lynn Govaert wrote: > I pressed enter to soo

Re: [R] print vectors with consecutive numbers

2014-08-22 Thread Bert Gunter
Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Fri, Aug 22, 2014 at 7:16 AM, James Wei wrote: > Hi Jorge, > > Thanks so much, it is working p

Re: [R] How can I let the dimension change via the circulation?

2014-08-23 Thread Bert Gunter
Huh? -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Sat, Aug 23, 2014 at 5:22 AM, 西风古道 wrote: > Here is the fraction of Matlab code for th

Re: [R] recording age into different variables

2014-08-25 Thread Bert Gunter
yourself. This one ships with R, but there are many other good ones on the web: http://cran.r-project.org/doc/manuals/R-intro.pdf Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainl

Re: [R] HANDLER FUNCTION PROBLEM

2014-08-29 Thread Bert Gunter
functions. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Fri, Aug 29, 2014 at 2:26 AM, Ashis Deb wrote: > Hello , > > > > I h

Re: [R] Testing general hypotheses on regression coefficients

2014-09-05 Thread Bert Gunter
y. But for statistics help, you should try a statistics list like stats.stackexchange.com instead. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll

Re: [R] splitting data

2014-09-08 Thread Bert Gunter
?tapply e.g. with(yourdata, tapply(A,list(year,month),sum,simplify=FALSE)) This assumes "sum them up" means summing each column separately. You were unclear as to exactly what you meant by this. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is

Re: [R] Efficient frontier

2014-09-15 Thread Bert Gunter
Please read the posting guide (link at bottom of message) to learn how to post coherently to get a useful response. I, at least, found your post to be unintelligible gibberish. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Informati

Re: [R] Changepoint analysis--is it possible to attribute changpoints to explanatory variables?

2014-09-16 Thread Bert Gunter
time series) and there probably is inertia in the system, too. So it may be complicated. That's why you need to spend time with someone who knows how to handle this. Econometricians tend to do this sort of thing I believe. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467

Re: [R] See the numeric codes of a factor

2014-09-19 Thread Bert Gunter
Talk to a local statistician or study a book on regression. You do not understand how regression works. In R, see ?contrasts . -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not w

[R] Fwd: Dynamic regex/sub changes to function

2014-09-22 Thread Bert Gunter
c(3,3,1)]] sin > class(body(f)) [1] "{" You should listen to your elders (Bill and Duncan) and **don't do this. ** Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly

Re: [R] division of col by the sum of the col

2014-09-23 Thread Bert Gunter
Neither. Unless I misjudge, you really really need to do your homework and read some basic R documentation -- e.g. An Intro to R, which ships with R (and also the Posting Guide) -- before posting here further. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data i

Re: [R] Writing .csv file

2014-09-24 Thread Bert Gunter
Read the help docs? ( pay attention to the row.names argument) -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Wed, Sep 24, 2014 at 2:00 PM,

Re: [R] adding rows

2014-09-25 Thread Bert Gunter
Inline. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Thu, Sep 25, 2014 at 3:28 PM, Sven E. Templer wrote: > see inline for another vecto

Re: [R] Ifelse statement on a factor level data frame

2014-09-28 Thread Bert Gunter
Inline. Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Sun, Sep 28, 2014 at 6:38 AM, Kate Ignatius wrote: > Strange that, > > I did

Re: [R] Plotting Categorical/Numerical Data?

2014-09-28 Thread Bert Gunter
om. You gave no details, but it may well be complex. 4. Or consult someone locally with knowledge of probability and statistics. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wis

Re: [R] Rolling window linear regression

2014-10-04 Thread Bert Gunter
Use ?loess instead. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Sat, Oct 4, 2014 at 12:09 AM, Grace Shi <1104271...@qq.com> wrote:

Re: [R] truncated normal

2014-10-05 Thread Bert Gunter
... yes. ... And do note that in sampling, truncated != censored. (They are often confused) Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll O

Re: [R] lattice add a fit

2014-10-07 Thread Bert Gunter
t;red" mod <- lm(y~x + w + z, data = data[subscripts,]) ... etc. } Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Tue,

Re: [R] lattice add a fit

2014-10-07 Thread Bert Gunter
Sorry, you do **NOT** need to pass the covariates if the data argument to lm is used. Or you can explicitly pass the covariates and subscript them in your panel function. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information i

Re: [R] lattice add a fit

2014-10-07 Thread Bert Gunter
Fit your model in the panel function using lm and plot the fits using ?panel.points, ?panel.lines, etc. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford St

Re: [R] understanding the no-label concept

2014-10-11 Thread Bert Gunter
No, you are wrong. Read the docs! -- start with "An Introduction to R" which ships with R. Please do not post further until after you have done your homework. x <- c(a=1,b=2,c=3) See also ?names. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374

Re: [R] plot hclust object

2014-10-28 Thread Bert Gunter
Copacetic cophenetics are a way To better see much genetics. ;-) -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Tue, Oct 28, 2014 at 1:2

Re: [R] Ways to get all function signatures of a library?

2014-10-29 Thread Bert Gunter
Perhaps the ?formals function in R is what you are looking for. Or maybe its (internal C) code. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll

Re: [R] MSE increased by increasing the sample size for Nadaraya-Watson kernel regression

2014-11-01 Thread Bert Gunter
essentially constant as an estimator of the population variance. Note: for nonparametric smoothers, mse is related to bandwidth choice also. This might change by default with different sample sizes. 3. In future, please post in plain text, not html, as the posting guide requests. Cheers, Bert Bert

Re: [R] Finding MLE

2014-11-02 Thread Bert Gunter
You do not appear to provide initial values for a and b , i.e. the "start" argument for mle. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford S

Re: [R] split a field in dependence of the semicolon

2014-11-03 Thread Bert Gunter
your data frame containing missings. There are various ways to do this, but "An Introduction to R" (ships with R) -- you have read it right? -- should provide the info you need. That should get you started. Others may provide a more complete solution, but that's much less fun. Cheer

Re: [R] Help with glm and glht for analysing count data

2014-11-04 Thread Bert Gunter
well to read and follow). Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Tue, Nov 4, 2014 at 6:37 AM, Mary Crossland wrote: > Dear all, &

Re: [R] how to determine power in my analysis?

2014-11-08 Thread Bert Gunter
Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Sat, Nov 8, 2014 at 3:49 AM, Kristi Glover wrote: > Hi R Users, > I was trying to determine

Re: [R] how to determine power in my analysis?

2014-11-08 Thread Bert Gunter
This is discussion is now off topic here. Either post elsewhere, e.g stats.stackexchange.com, or consult your local statistician for help, as I previously suggested. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not know

Re: [R] Working with data-frame

2014-11-09 Thread Bert Gunter
above functions already have this "built" in, so you don't need to do this explicitly, although my impression is that it may be a tad more efficient to do it the long way. But don't quote me on this! Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 46

Re: [R] Is it possible to define another kind of NA

2014-11-09 Thread Bert Gunter
Ouch! The values are **NOT** missing -- they are (left) censored, and need to be handled by appropriate censored data methods. I suggest you (all!) either read up on this or consult someone locally who has knowledge of such methods. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650

Re: [R] Error: (list) object cannot be coerced to type double; on running the following code in R ver 3.1.2

2014-11-11 Thread Bert Gunter
ts -- this is about automatic coercion -- before dispensing false advice. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Tue, Nov 11, 2014

Re: [R] (no subject)

2014-11-11 Thread Bert Gunter
Start reading: http://cran.r-project.org/doc/manuals/R-intro.pdf (or hire a consultant?) Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll O

Re: [R] Two dimensional likelihood surface plot

2014-11-11 Thread Bert Gunter
iently bad set of (least squares)starting values, as ?loess explicitly says. It can also be slow. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stol

Re: [R] how to compure R-squared in glm

2014-11-16 Thread Bert Gunter
ussion here. Consult a local statistician, do your glm homework, or post on stats.stackexchange.com for statistical follow-ups. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly no

Re: [R] Conditionally replace multiple rows in a dataframe

2014-11-20 Thread Bert Gunter
to whether this is actually preferable and why or why not) Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Wed, Nov 19, 2014 at 11:23 PM, Jef

Re: [R] Conditionally replace multiple rows in a dataframe

2014-11-20 Thread Bert Gunter
Sorry, just noticed the typo: it's is.na() of course. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Thu, Nov 20, 2014 at 4:04 AM, Ber

Re: [R] Interpreting ANOVA / quantile regression results

2014-11-21 Thread Bert Gunter
This is off topic here.Post to a statistics forum like stats.stackexchange.com instead -- or talk to your professor or TA (if you're a student). Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And kn

Re: [R] Generate random numbers under constrain

2014-11-22 Thread Bert Gunter
Well, if their sum must be < 1 they ain't random... But anyway... given n randnums <- function(n) { Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Sto

Re: [R] Generate random numbers under constrain

2014-11-22 Thread Bert Gunter
(Hit send key by accident before I was finished ...) Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Sat, Nov 22, 2014 at 7:14 AM, Bert Gunter wrote:

Re: [R] DOT PLOT help!!

2014-11-24 Thread Bert Gunter
://www.rcommander.com/ , is one that is well documented and well developed. Pestering us on this list to tell you how to do everything from the R command line is unfair to us and a foolish strategy for you. IMO only, of course. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374

Re: [R] Help me

2014-11-24 Thread Bert Gunter
No. See the posting guide link below for how to ask an intelligible question. Better yet, do your homework (is that what this is?) yourself: http://cran.r-project.org/doc/manuals/R-intro.pdf -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not inform

Re: [R] covariate or predictor

2014-11-25 Thread Bert Gunter
understand the use of contrasts (= dummy variables). As you said, either stackexchange or perhaps a local consultant is probably where she should be seeking advice, Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not

Re: [R] please please unsubscribe!!!!!!!!!!!!!!!!

2014-11-29 Thread Bert Gunter
Please please follow the mailing list link below and unsubscribe yourself there as directed, not here. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll

Re: [R] Creating submatrices from a dataframe, depending on factors in sample names

2014-12-01 Thread Bert Gunter
"3.1" "3.2" "3.3" > x[upper.tri(x)] [1] "1.2" "1.3" "2.3" > x[upper.tri(x,diag=TRUE)] [1] "1.1" "1.2" "2.2" "1.3" "2.3" "3.3" This gives you a vector all possible pairs (

<    7   8   9   10   11   12   13   14   15   16   >