Re: [R] Using ggplot2 geom_path() in a grouped variable

2019-02-12 Thread David Winsemius
> On Feb 12, 2019, at 2:09 PM, Norberto Hernandez > wrote: > > Hi! I am trying to make a scatter/path graph in one variable that is > divided in two groups (variable Control), but when I use the > geom_path() option, the line continues from the group one to the group > two, and I wasn't able

Re: [R] Using ggplot2 geom_path() in a grouped variable

2019-02-12 Thread David Winsemius
To Norberto; Your code _probably_ would have succeeded if you had used color as a grouping argument, or you could have used group or linetype or probably others, but I find locating the listing of ggplot2 "behavioral" parameters rather frustrating. These details are not to be found in any of ?a

Re: [R] Saving and reloading function in a package

2019-02-15 Thread David Winsemius
Isn't this much more on topic with the package development list? -- David > On Feb 15, 2019, at 6:37 AM, Dimitrios Stasinopoulos > wrote: > > I would like to put a graphic background to a model diagnostic plot. > The background is created with plot()/lines() but it takes time. > My solution

Re: [R] Help formulating gamm with repeated measures and spatial autocorrelation using mgcv package

2019-02-26 Thread David Winsemius
Crossposting is explicitly advised against in the rhelp posting guide. Furthermore, this is the wrong list for this sort of question within the r-project hierarchy of mailing lists. I do not believe an answer is possible in the absence of the data (and note the request some would say deman

Re: [R] add points to lattice cloud plot (3D scatter)

2019-02-28 Thread David Winsemius
On 2/28/19 7:56 AM, Jeff Newmiller wrote: You are missing the point... lattice assembles the entire data set at once so it can adjust and synchronize all of the scales and then it generates an object that can be printed to a device. This approach is entirely incompatible with the base graphi

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

2019-03-10 Thread David Winsemius
Seems rather likely that 999 is not really a measured value but rather is a missing value indicator. -- David. On 3/10/19 1:54 PM, wong bowie wrote: I wish to calculate the weight of evidence of a variable x, which is positively skewed, with over 6000 of the observations are 999 but only 200

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

2019-03-10 Thread David Winsemius
ne to specify a value such as 999 to be missing but R needs to have it changed to NA is.na(Table$pdays) <- Table$pdays == 999 -- David > > David Winsemius <mailto:dwinsem...@comcast.net>> 於 2019年3月10日 週日 下午10:48寫道: > > Seems rather likely that 999 is not r

Re: [R] About wilcox test

2019-03-12 Thread David Winsemius
On 3/12/19 9:23 AM, Thomas Subia via R-help wrote: Javid wrote: "I have two set of data in excel: A column( 16.38, -31, -16.77, 127, -57, 23.44 and so on) B column ( -12, -59.23, -44, 34.23, 55.5, -12.12 and so on) If those are in an excel spreadsheet than you need to indicate exactly how

Re: [R] A general question about using Bayes' Theorem for calculating the probability of The End of Human Technological Civilisation

2019-03-19 Thread David Winsemius
Rhelp is not a forum for discussions of statistics. Instead it is for persons who have specific questions about the use of R. Please read the list info page where you started the subscription process. And do read the Posting Guide. Both these are linked at the bottom of this response. There

Re: [R] A general question about using Bayes' Theorem for calculating the probability of The End of Human Technological Civilisation

2019-03-19 Thread David Winsemius
e approves them), and if these 'statistics discussion' questions are inappropriate to the mission (as described), then...why isn't the 'moderator' (him/her/they) blocking on submission? On 3/19/2019 1:59 PM, David Winsemius wrote: Rhelp is not a forum for discussions o

Re: [R] A general question about using Bayes' Theorem for calculating the probability of The End of Human Technological Civilisation

2019-03-19 Thread David Winsemius
On 3/19/19 12:49 PM, Jeff Newmiller wrote: Highly off topic. Try StackOverflow. As it stands it's off-topic for SO. (You would just be making more work for those of us who know the rules but need 4 close votes for migration.)  Better would be immediately posting at CrossValidated.com (i.e

Re: [R] Rank ANCOVA

2019-03-21 Thread David Winsemius
If you are looking for robust multivariate or multivariable methods, then review the Robust Methods Task View: https://cran.r-project.org/view=Robust -- David. On 3/13/19 4:06 AM, Jackson, Daniel wrote: Hi Frank and Dennis I am in a similar situation but I would prefer to use a proportion

Re: [R] counting unique values (summary stats)

2019-03-21 Thread David Winsemius
On 3/21/19 3:31 PM, reichm...@sbcglobal.net wrote: r-help I have the following little scrip to create a df of summary stats. I'm having problems obtaining the # of unique values unique=sapply(myData, function (x) length(unique(x), replace = TRUE)) I just looked up

Re: [R] Regression Modeling Strategies and the R rms Package Short Course 2019

2019-03-24 Thread David Winsemius
Hi Graeme; I took the course about ten years ago. I did so after getting a Masters in Epidemiology from the University of Washington and doing very well in all my stats courses and submitting my thesis work on solving regression problems with stratified sampling using bootstrap methods. So

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

2019-03-28 Thread David Winsemius
On 3/27/19 3:43 PM, Bernard Comcast wrote: To follow on Jeff, is there a function to do 2-D (double) numerical integration in R? Packages pracma and cubature offer a variety of solutions to that task. -- David. Bernard Sent from my iPhone so please excuse the spelling!" On Mar 27, 2

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

2019-03-28 Thread David Winsemius
also has a "sliced" approach to disply of 2d densities that I found informative and attractive. -- David. Bernard McGarvey Director, Fort Myers Beach Lions Foundation, Inc. Retired (Lilly Engineering Fellow). On March 28, 2019 at 1:40 PM David Winsemius wrote: On

Re: [R] Help installing netReg

2019-04-14 Thread David Winsemius
On 4/14/19 3:20 PM, Spencer Brackett wrote: Good evening, I am having problems with downloading the package used to generate regression models on R. The following is the error message I received. I tried installing BiocManager instead as suggested, but this too did not work. Any ideas? The

Re: [R] packages across different versions of R

2014-07-22 Thread David Winsemius
rly named directory and then run: update.packages(checkBuilt=TRUE,ask=FALSE) # I added that second argument b/c I have a large collection (This doesn't work for packages only available on r-forge or github.) -- David Winsemius Alameda, CA, USA _

Re: [R] Importing random subsets of a data file

2014-07-23 Thread David Winsemius
ionaldiversity.org > > __ > 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-con

Re: [R] Retrieve Axis coordinates from map

2014-07-24 Thread David Winsemius
automatically to be able to > generate the legend within the function. > > Thanks > > Julien > > > > [[alternative HTML version deleted]] r-help is a plain text mailing list. > > PLEASE do read the posting guide http://www.R-project.org/posting-guide.h

Re: [R] Using R to Compute Covariance

2014-07-26 Thread David Winsemius
ultiplying it by 2 and then dividing by 3. However, when I get that I still do not get 41/16. What am I missing? -- David Winsemius, MD Alameda, CA, USA __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEA

Re: [R] Need(?) way to create a "picture" (or "plot") of a data.frame

2014-07-31 Thread David Winsemius
ns is sos::findFn. I'm pretty sure that Baptiste Augie has written a function or a package to do this, and I've seen several worked examples on StackOverflow. This is something I fairly quickly found: http://rpubs.com/baptiste/ftableGrob -- David Winsemius Alameda, CA, USA ___

Re: [R] keep information on the number of warnings

2014-07-31 Thread David Winsemius
> the following message: > > "There were 36 warnings (use warnings() to see them)" > > How can I store the number 36 in a variable in function Y? > > In other words, how can I extract the information on the number of warnings > generated? length(w

Re: [R] mutually exclusive events

2014-08-02 Thread David Winsemius
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 Alameda, CA, USA __ R-help@r-project.org ma

Re: [R] keep information on the number of warnings

2014-08-04 Thread David Winsemius
cio C2, 6º Piso, Sala 2.6.04/07 > Campo Grande > 1749-016 Lisboa, Portugal > Tel: +351 21 750 00 00 (ext: 22607) > Fax: +351 21 750 00 28 > > > > > > > > > > > [[alternative HTML version deleted]] > > _

Re: [R] big data?

2014-08-05 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 Alameda, CA, U

Re: [R] R script in batch mode — Echoing messages to windows shell

2014-08-05 Thread David Winsemius
tive HTML version deleted]] This was just seen on SO: "If my guess is correct, shell("echo Some test > CON") should work. – Harry Johnston 1 min" If you had read the Posting Guide, you might have seen that cross-posting to multiple sites is deprecated. -- David W

Re: [R] lattice scales format dates

2014-08-06 Thread David Winsemius
45) ) ) > > > > dat <- data.frame(var1=1:10, > date.val=as.POSIXct(seq(as.Date("1910/11/1"), as.Date("1911/8/1"), > "months") )) > > barchart(var1~date.val, data=dat, horizontal=F, > scales=list(x=list(rot=45, format="%Y-%b")))

Re: [R] problem with labeling plots, possibly in font defaults

2014-08-07 Thread David Winsemius
es-Bold" "Times-Italic" "Times-BoldItalic" > > > This is clearly not a big problem (at least for now) since I can put labels > on plots by running par(), but if it is indicative of a larger underlying > problem (or if there is a simple fix) I

Re: [R] how to process multiple data files using R loop

2014-08-08 Thread David Winsemius
the output is : > >> head(P_3_utr_source_data,2) > V1 > 1 1 > 2 1 >> head(P_5_utr_source_data,2) > V1 > 1 1 > 2 1 >> > . > > . > . > > > > Could anybody help me with this? > > Thank you very much for your time:) > [[alternativ

Re: [R] circlize package: some error message

2014-08-08 Thread David Winsemius
> Wonder anyone could please help me out with this? > > Thank you very much for your time.:) > > [[alternative HTML version deleted]] Your failure to post in plain text makes this message very difficult to read. Unless there i

Re: [R] Installing manual package problem

2014-08-08 Thread David Winsemius
> (...R/win-library/3.0) I see the glassbox folder there. > > I'm new to using packages not from the CRAN list so I'm trying to learn > fast. I tried some searching and this seems to be what I'm suppossed to > do, but perhaps I need to use dev mode

Re: [R] R Package for Text Manipulation

2014-08-09 Thread David Winsemius
ative HTML version deleted]] And do : > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. -- David Winsemius Alameda, CA, USA __ R-help@r-project.org ma

Re: [R] "Best" way to merge 300+ .5MB dataframes?

2014-08-10 Thread David Winsemius
) a b 1 1 2 > rbind( data.frame(a=1, b=2), data.frame(a=1, b=2) ) a b 1 1 2 2 1 2 > Now with that in mind, what is the idiomatic way? Do people usually do > something else because it is /faster/ (by some definition)? > > Kind regards, > -- David Winsemius Alameda, C

Re: [R] efficient way to replace a range of numeric with a integer in a matrix

2014-08-11 Thread David Winsemius
> [4,]1112 >>> >>> I wrote a function to do such thing: >>> >>> fun <- function(x) { >>> if (is.na(x)) { >>> NA >>> } else if (x > 0.0 && x <= 1.0) { >>> 1 >>>

Re: [R] "Best" way to merge 300+ .5MB dataframes?

2014-08-11 Thread David Winsemius
r implementation would work well on that system. Well, if > there were a version of R for it. It is a branded UNIX system which > was designed to be totally __and only__ POSIX compliant, with few > (maybe no) extensions at all. IOW, it stinks. No, it can't be > replaced. It is the z

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

2014-08-12 Thread David Winsemius
t trust. I've used it with the `cobs` function from the package of the same name to implement the monotonic constraint. I think there is a worked example in the quantreg package, but since I bought Koenker's book, I may be remembering from there. -- David Winsemius Alameda, CA, USA

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

2014-08-13 Thread David Winsemius
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 Winsemius > wrote: >>

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

2014-08-14 Thread David Winsemius
bolomics > Food Quality and Nutrition > Fondazione Edmund Mach > > > > On 08/14/2014 02:02 AM, David Winsemius wrote: >> >> On Aug 12, 2014, at 8:40 AM, Bert Gunter wrote: >> >>> PI's of what? -- future individual values or mean values? >&

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

2014-08-14 Thread David Winsemius
On Aug 14, 2014, at 9:06 AM, David Winsemius wrote: > > On Aug 14, 2014, at 7:17 AM, Jan Stanstrup wrote: > >> Thank you very much for this snippet! >> >> I used it on my data and indeed it does give intervals which appear quite >> realistic (script and

Re: [R] Generating a polygon around points

2014-08-15 Thread David Winsemius
;> Bob > > __ > 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] unable to use functions require DLL from package

2014-08-15 Thread David Winsemius
I'm using R3.1.0 on windows 8 and I've never used R before, so I might have > made some terrible newby error (I have programmed quite a bit before, but in > C > and Matlab). > > Thank you! > > Best, David Winsemius Alameda, CA, USA _

Re: [R] Generating a polygon around points

2014-08-15 Thread David Winsemius
On Aug 15, 2014, at 3:11 PM, Rolf Turner wrote: > On 16/08/14 05:10, Jeff Newmiller wrote: >> I use RSiteSearch regularly with no problems. Perhaps I have just had >> a lucky streak? I wonder what the odds are... :-) > > > > >> On August 15, 2014 8:55:40

Re: [R] unable to use functions require DLL from package

2014-08-15 Thread David Winsemius
On Aug 15, 2014, at 2:56 PM, Charlotte de Vries wrote: > David Winsemius comcast.net> writes: > >> >> >> On Aug 15, 2014, at 8:06 AM, Charlotte de Vries wrote: >> >>> Hey there! >>> >>> I'm having problems with the same

Re: [R] unable to use functions require DLL from package

2014-08-16 Thread David Winsemius
that prevent it from being a current package and working with the current mgcv. (I also checked to see if it had been submitted to CRAN in that era: http://cran.r-project.org/bin/windows/contrib/1.7/ But it's not there.) I looked at the code and it does not require compilation. You can pr

Re: [R] Making my own graphics device

2014-08-17 Thread David Winsemius
to the S language (at least on PC's) would be TIBCO. They do list the acronym `hpgl` on one page in the S+ documentation, but it is in the "deprecated functions" section (p 39): http://www.uni-koeln.de/themen/statistik/software/s/v81/functionguide.pdf Good luck. -- Dav

Re: [R] Making my own graphics device

2014-08-17 Thread David Winsemius
oducible code. >> > > -- > Dr Paul Murrell > Department of Statistics > The University of Auckland > Private Bag 92019 > Auckland > New Zealand > 64 9 3737599 x85392 > p...@stat.auckland.ac.nz > http://www.stat.auckland.ac.nz/~paul/ > > ___

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

2014-08-18 Thread David Winsemius
difference. > > > > On 08/14/2014 06:06 PM, David Winsemius wrote: >> >> On Aug 14, 2014, at 7:17 AM, Jan Stanstrup wrote: >> >>> Thank you very much for this snippet! >>> >>> I used it on my data and indeed it does give intervals which

Re: [R] DateTime wrong when exporting to csv in R

2014-08-20 Thread David Winsemius
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 Alameda, CA, USA __ R-help@r-project.org mailing list https

[R] Fwd: DateTime wrong when exporting to csv in R

2014-08-20 Thread David Winsemius
The image file I prepared and attached did not make it through. Trying again. -- David Begin forwarded message: > From: David Winsemius > Subject: Re: [R] DateTime wrong when exporting to csv in R > Date: August 20, 2014 1:48:14 PM PDT > To: Sneha Bishnoi > Cc: r-help &g

Re: [R] DateTime wrong when exporting to csv in R

2014-08-20 Thread David Winsemius
On Aug 20, 2014, at 1:48 PM, David Winsemius wrote: > > On Aug 20, 2014, at 4:53 AM, Sneha Bishnoi wrote: > >> Hi All! >> >> This seems to be trival but I am not able to find a solution for it. >> I have a dataframe with datetime columns in form of ("%

Re: [R] How to view the whole dataset that is imported through sasxport.get

2014-08-21 Thread David Winsemius
79 [2,]2468 10 > attr(x,"dim") [1] 2 5 > attr(x,"dim") <- NULL > x [1] 1 2 3 4 5 6 7 8 9 10 It also appears the there is a `label<-` function, so you could probably use that to NULL them out. -- David Winsemius Alameda, CA, USA

Re: [R] How to view the whole dataset that is imported through sasxport.get

2014-08-21 Thread David Winsemius
'labelled' atomic [1:1991] -1.1 0.8 1.5 3 4 ... > .. ..- attr(*, "label")= chr "Actual Relative Time" > If you don't want those to be "labelled"-class, then just: dfrm2 <- data.frame( lapply(dfrm1, unclass) ) > > > On Thu, Aug 2

Re: [R] Help on installing "R" packages in a Citrix

2014-08-22 Thread David Winsemius
-packages" for help, but was told this was > the correct Forum for info of this type. Have you installed R? Which OS? Which version? -- David Winsemius Alameda, CA, USA __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo

Re: [R] Help on installing "R" packages in a Citrix

2014-08-22 Thread David Winsemius
quot;,"rcom"),repos="http://rcom.univie.ac.at/download",lib=.Library) > library(rcom) > comRegisterRegistry() > > i get the following message: > > > Please can you help me? > > Many Thanks > > Chirag Patel > > > -Original

Re: [R] Help on installing "R" packages in a Citrix

2014-08-22 Thread David Winsemius
On Aug 22, 2014, at 9:13 AM, David Winsemius wrote: > > On Aug 22, 2014, at 7:37 AM, Chirag Patel wrote: > >> Hi David >> >> I have installed R Version 2.15.1 on the image but having problems install >> rcom and rscproxy. > > Those are commercial pack

Re: [R] converting dataframe into multiple time series

2014-08-23 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 Alameda, CA, USA

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

2014-08-24 Thread David Winsemius
ect.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. >>> >> >> >> >> -- &

Re: [R] Filled vector contours

2014-08-24 Thread David Winsemius
s. And this > seems to be the case also of R's 'contourLines'. > > But filling the areas requires a further post-processing (joining line > fragments, joining contours, set the correct direction to preserve holes). > And if I understand correctly, this is

Re: [R] Preventing loading of user packages

2014-08-25 Thread David Winsemius
;> which contains a line: >> >> R_LIBS_USER=${R_LIBS_USER-'~/Library/R/3.0/library'} >> >> but commenting this out does not prevent R from loading from that >> location. >> >> is there a way to prevent loading from that location? >> &

Re: [R] Plot survreg and intcox

2014-08-26 Thread David Winsemius
t; __ > 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 Winsemi

Re: [R] plot for "survreg" and "intcox" (rewritten)

2014-08-26 Thread David Winsemius
p=NA > > ______ > 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 commente

Re: [R] Mixed sorting/ordering of strings acknowledging roman numerals?

2014-08-26 Thread David Winsemius
t;chr V" > [9] "chr IV" "chr III" "chr II" "chr I" > >> sort(y) > [1] "chr I""chr II" "chr III" "chr IV" > [5] "chr IX" "chr V""chr VI" "chr VII&quo

Re: [R] NA's introduced by coercion

2014-08-26 Thread David Winsemius
er_) ) [1] 2 1 NA So it is not the NA's in that character vector but rather values that were coerce to NA because the conversion could not be accomplished. -- David Winsemius Alameda, CA, USA __ R-help@r-project.org mailing list https://stat.eth

Re: [R] lattice: packet.number() versus panel.number()

2014-08-26 Thread David Winsemius
, but I suspect the packet numbers are sequential increasing. -- David Winsemius Alameda, CA, USA __ 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/post

Re: [R] Where to find source of C_pbinom?

2014-08-27 Thread David Winsemius
oaded into R. Where can I find the source code of C_pbinom? www.r-project.org/doc/Rnews/Rnews_2006-4.pdf Pages 43-54 -- David Winsemius, MD Alameda, CA, USA __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE d

Re: [R] simple plotting wcor (Rssa) question

2014-08-27 Thread David Winsemius
not show up as F1,F2 in your example. It doesn't make much sense to request the at values of 30 and 40 since you only gave it a matrix that was 20 x 20. -- David Winsemius Alameda, CA, USA __ R-help@r-project.org mailing list https://stat.ethz.c

Re: [R] Issue with formula conversion

2014-08-27 Thread David Winsemius
, as shown below: > >> as.formula(ranEff) > ~1 + Age > > the "pdCompSymm" is lost in the conversion. Any solutions? as.formula(paste("~",ranEff)) ~pdCompSymm(~1 + Age) -- David Winsemius Alameda, CA, USA __ R-help@r-

Re: [R] Issue with formula conversion

2014-08-27 Thread David Winsemius
thout a tilde (~) at the beginning. > > On Wed, Aug 27, 2014 at 3:34 PM, David Winsemius > wrote: >> >> On Aug 27, 2014, at 11:19 AM, Gang Chen wrote: >> >>> A random effect formulation for R package nlme is read in as a string >>> of cha

Re: [R] Issue with formula conversion

2014-08-27 Thread David Winsemius
ndom=list(Block=eval(tt) )) > fm Linear mixed-effects model fit by REML Data: Oats Log-restricted-likelihood: -296.5209 Fixed: yield ~ nitro (Intercept) nitro 81.8722273.7 snipped rest of output. -- David > > works, I don't know the 'pdCompSymm

Re: [R] Issue with formula conversion

2014-08-27 Thread David Winsemius
ety-1)' >> >> The following >> >> (ranEff1 <- as.formula(ranEff)) >> ~Variety - 1 >> >> is not what I want. Even though >> >> fm <- lme(yield ~ nitro, data=Oats, random=list(Block=pdCompSymm(ranEff1))) >> >> works,

Re: [R] R packages for power analysis

2014-08-28 Thread David Winsemius
> and provide commented, minimal, self-contained, reproducible code. David Winsemius Alameda, CA, USA __ 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.

Re: [R] dont remenber my password

2014-08-29 Thread David Winsemius
>> From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- >> project.org] On Behalf Of Ernesto Villarino >> Sent: Friday, August 29, 2014 9:32 AM >> To: r-help@r-project.org >> Subject: [R] dont remenber my password >> >> >> snip David Winsemi

Re: [R] Question regarding the discrepancy between count model parameter estimates between "pscl" and "MASS"

2014-08-29 Thread David Winsemius
ispersion > parameter for Negative Binomial(1.1303) family taken to > be >>> 0.4549082) >>> >>> >Null deviance: 76.159 on 115 degrees of > freedom >>> Residual deviance: 35.101 > on 104 degrees of freedom >>>(296 observation

Re: [R] Unexpected behavior when giving a value to a new variable basedon the value of another variable

2014-08-30 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 Alameda, CA, USA

Re: [R] Unexpected behavior when giving a value to a new variable basedon the value of another variable

2014-08-30 Thread David Winsemius
On Aug 30, 2014, at 7:38 PM, David Winsemius wrote: On Aug 29, 2014, at 8:54 PM, David McPearson wrote: On Fri, 29 Aug 2014 06:33:01 -0700 Jeff Newmiller > wrote One clue is the help file for "$"... ?" $" In particular there see the discussion of character i

Re: [R] Depth vs Temp graph for different transects

2014-09-01 Thread David Winsemius
Cell: 082 080 9638 E-mail: tsonne...@gmail.com [[alternative HTML version deleted]] __ 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

Re: [R] depth of labels of axis

2014-09-01 Thread David Winsemius
It did not seem to me that padj=0 produced "top alignment". -- David Winsemius, MD Alameda, CA, USA __ 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

Re: [R] Detect expired RSQLiteConnection?

2014-09-02 Thread David Winsemius
Noodling through the S4 and then the S3 code I found: ?isIdCurrent as the test in `sqliteCloseConnection` > -- David Winsemius Alameda, CA, USA __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read th

Re: [R] Simulating from a Weibull distribution

2014-09-02 Thread David Winsemius
rvreg docs says: # Weibull parametrisation y<-rweibull(1000, shape=2, scale=5) survreg(Surv(y)~1, dist="weibull") # survreg scale parameter maps to 1/shape, linear predictor to log(scale) David Winsemius Alameda, CA, USA __ R-help@r-pr

Re: [R] depth of labels of axis

2014-09-04 Thread David Winsemius
xpression(E[g]), "E", expression(E[j]), >>> "E", expression(E[t]))) >>> >>> you may notice that the "E" within labels of axis(1) are not at the same >>> depth. So the vision of axis(1) labels is something like wave. >>> >>

Re: [R] Help with regression

2014-09-04 Thread David Winsemius
ery easy to post in plain-text from gmail. > > [[alternative HTML version deleted]] > > __ > 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

Re: [R] Operator proposal: %between%

2014-09-04 Thread David Winsemius
LEASE do read the posting guide http://www.R-project.org/posting-guide.html >> and provide commented, minimal, self-contained, reproducible code. > > > > -- > Gregory (Greg) L. Snow Ph.D. > 538...@gmail.com > > __ > R

Re: [R] Subset a column with specific characters

2014-09-04 Thread David Winsemius
la.edu/stat/stata/notes/hsb2.dta') > NROW( hsb2[ grepl("ca", hsb2$prog) & !grepl("ic", hsb2$race) , ] ) [1] 120 -- David Winsemius Alameda, CA, USA __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/l

Re: [R] depth of labels of axis

2014-09-05 Thread David Winsemius
5, xaxt = "n") axis(1, at = 1:5, labels = updown(lbl) ) abline(h=.7, xpd=TRUE, lty=3) > > Thanks a lot. > > Best regards, > Jinsong > > >> >> -Original Message- >> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On

Re: [R] depth of labels of axis

2014-09-05 Thread David Winsemius
On Sep 5, 2014, at 4:53 PM, David Winsemius wrote: > > On Sep 5, 2014, at 4:15 PM, Jinsong Zhao wrote: > >> On 2014/9/4 14:58, David L Carlson wrote: >>> The problem with this approach is that the horizontal positioning of the >>> labels is based on the

[R] Why does debugging print() change output of function?

2014-09-06 Thread David Winsemius
[[3]] [1] TRUE TRUENA FALSE FALSENA TRUE This actually seemed to be somewhat successful, but when ... Now if I take out the `print()` call for 'tests', I get an different answer: > mtest.data.frame <- + function (x, ..., drop=FALSE) + { tests

Re: [R] Why does debugging print() change output of function?

2014-09-06 Thread David Winsemius
eval(e, x, enclos) >if (!is.logical(r)) { >stop("'tests' must be logical") >} >r & !is.na(r) >})) >} >r > } > > used as: >> m2(testdata, group2=="UNC", Eleve

Re: [R] Solving equations

2014-09-07 Thread David Winsemius
tions based on data available in the > puzzle. > > How do I solve a set of such equations ? I need to find out the values of x > and w for each such equation. > > I know that here the equations are simple because the puzzle can be > simplfied. > > > Thanks, > Moha

Re: [R] Mixed sorting/ordering of strings acknowledging roman numerals?

2014-09-07 Thread David Winsemius
-9]*([eE][\\+\\-]{0,1}[0-9]+\\.{0,1}[0-9]*){0,1})", > 59,62d48 > > This difference is so small that the above could now be an option to > mixedorder() with minimal overhead added, e.g. mixedorder(y, > type=c("decimal", "roman")). One could even imagine addi

Re: [R] splitting data

2014-09-08 Thread David Winsemius
> > I initially tried to use "dcast" command but for no use. > This is typical use for base function aggregate: > with( dfrm, aggregate(dfrm[4:7], dfrm[1:2], FUN=sum)) year month A B C D 1 1961 1 78 102 138 258 >

Re: [R] splitting data

2014-09-08 Thread David Winsemius
On Sep 8, 2014, at 2:01 PM, David Winsemius wrote: > > On Sep 8, 2014, at 12:08 PM, eliza botto wrote: > >> >> > > This is typical use for base function aggregate: > >> with( dfrm, aggregate(dfrm[4:7], dfrm[1:2], FUN=sum)) > year month A

Re: [R] Why does debugging print() change output of function?

2014-09-09 Thread David Winsemius
On Sep 9, 2014, at 4:07 PM, peter dalgaard wrote: On 07 Sep 2014, at 00:31 , David Winsemius wrote: The goal: to create a function modeled after `subset` (notorious for its non- standard evaluation) that will take a series of logical tests as unquoted expressions to be evaluated in

Re: [R] How to plot soil moisture data as a contour plot

2014-09-10 Thread David Winsemius
. Thanks in advance! > > > Greetings, > -- > Thiago V. dos Santos > PhD student > Land and Atmospheric Science > University of Minnesota > http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm > Phone: (612) 323 > 9898 > [[alternative HTML versio

Re: [R] Margins to fill matrix

2014-09-11 Thread David Winsemius
: >>> >>> Hi, >>> >>> I have two vector of margins. Now I want to create "fill" matrix that >>> reflects the margins. >>> >>> seats <- c(17,24,28,30,34,36,40,44,46,50) >>> mandates <- c(107,23,24,19,112,19,2

Re: [R] how to run Fleming-Harrington weighted log-rank test in R

2014-09-12 Thread David Winsemius
ming-Harrington") [[alternative HTML version deleted]] Learn to search and learn to post in plain text. PLEASE do read the posting guide http://www.R-project.org/posting-guide.html -- David Winsemius, MD Alameda, CA, USA __ R-help@r-project.or

Re: [R] legend with math (greek letters) symbols

2014-09-13 Thread David Winsemius
nks, -Sergio. __ 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] ggplot y-axis labels are not continuous

2014-09-14 Thread David Winsemius
z.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 Alameda, CA, USA __ R-help@r-project.org mailing list https://stat.eth

Re: [R] simulation data in SEM

2014-09-14 Thread David Winsemius
xi<-mvrnorm(1,mu=c(0,0),phi) + #Generate the fixed covariates + co<-rnorm(1,0,1) snipped -- David Winsemius, MD Alameda, CA, USA __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] simulation data in SEM

2014-09-14 Thread David Winsemius
rix(data=NA,nrow=N,ncol=P); p<-numeric(P); v<-numeric(P) PLEASE read the Posting Guide. I will not respond to offlist messages from you in the future. -- David. > > > many thanks again > > On 14 September 2014 19:37, David Winsemius wrote: > > On Sep 14, 2014

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