Re: [R] read.fortran format

2016-05-27 Thread William Dunlap via R-help
acter vector. > >For a multiline record it should be a list with a character vector > >for each line. > > > > ​ > > > > ​I think (not sure) you need: > > > > > mydata<-read.frotran("G:/Journals/Disk1/12_restat_95/estimate/GROUPD.DAT", > > > list(c("1X",&q

Re: [R] read.fortran format

2016-05-28 Thread William Dunlap via R-help
uot;,"5F8.4","F10.6" > >> > /"1X","2F6.0","3E15.9","F8.0","F5.2","F5.3" > >> > /"1X","F7.0","2E15.9","F9.4","F5.3"), > >> > ​​

[R] matrix indexing/subset error

2016-05-30 Thread Carl Sutton via R-help
long Many thanks for your assistance. Carl Sutton CPA [[alternative HTML version deleted]] ______ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

[R] SEM GFI

2016-05-31 Thread VINAY KULKARNI via R-help
but not able to get the solution. Please let me know what might be the reason. Thanks,Vinay [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] SEM GFI

2016-06-01 Thread VINAY KULKARNI via R-help
ent ><-> brand1_Different,the11,NAbrand1_like_me  <-> >brand2_Def_Drink,the12,NAbrand1_Like <-> brand1_Love,the13,NAVal_brand1 <-> >Val_brand1,1,NA  > opt <- options(fit.indices = c("GFI", "AGFI", "RMSEA", >"NFI", &qu

Re: [R] Struggling trying to plot points on boxplot

2016-06-02 Thread William Dunlap via R-help
o be less than the number of boxplots. I have done this > previously, years ago, but can not find my notes and can not reconstruct > what I did. > > Help is appreciated. > > Regards, > Tom > > [[alternative HTML version deleted]] > > _

[R] Fw: i need to install "Rclimdex"........Senior Scientists , BARI

2016-06-05 Thread Ranjit Sen via R-help
how to install in my computer. With thank and very best regards. Sincerely,Dr Ranjit Sen, Senior Scientific Officer, Soil Science Division,  Bangladesh Agricultural Research Institute(BARI). [[alternative HTML version deleted]] ______ R-help@r

[R] Fw: i need to install "Rclimdex"........Senior Scientists , BARI

2016-06-05 Thread Ranjit Sen via R-help
). [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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] detecting if a variable has changed

2016-06-05 Thread William Dunlap via R-help
> > http://www.hep.by/gnu/r-patched/r-exts/R-exts_122.html > http://homepage.stat.uiowa.edu/~luke/R/references/weakfinex.html > > As far as I can tell, weakrefs are only available via the C API. Is > there a way to do what I want in R without resorting to C code? Is >

[R] Revolutions blog: May 2016 roundup

2016-06-07 Thread David Smith via R-help
Since 2008, Microsoft (formerly Revolution Analytics) staff and guests have written about R every weekday at the Revolutions blog: http://blog.revolutionanalytics.com and every month I post a summary of articles from the previous month of particular interest to readers of r-help. And in case

Re: [R] Power Calculation: Binomial Proportions (2 sided exact test for equivalence)

2016-06-07 Thread BONACHE Adrien via R-help
"two.sided")) Maybe you could use G*power :Lower case 1-beta = 0.9557076Upper one 1-beta= 0.9134817 Hope this help, Adrien Bonache. De : Munjal Patel À : r-help@r-project.org Envoyé le : Mardi 7 juin 2016 18h08 Objet : [R] Power Calculation: Binomial Proportions (2

Re: [R] Power Calculation:2-sided exact equivalence test for Binomial Proportions

2016-06-07 Thread BONACHE Adrien via R-help
Hi again, In fact the pwr package do not provide exact test power. So you should definitely use g*power. Sincerely, Adrien. De : Munjal Patel À : r-help@r-project.org Envoyé le : Mardi 7 juin 2016 18h26 Objet : [R] Power Calculation:2-sided exact equivalence test for Binomial

Re: [R] apply and cousins

2016-06-08 Thread William Dunlap via R-help
gt; Or does apply etc only operate on each row at a time, independently of > other rows? > > > Best wishes > > John > > John Logsdon > Quantex Research Ltd > +44 161 445 4951/+44 7717758675 > > __ > R-help@r-project

Re: [R] Reading and converting time data via read.table

2016-06-09 Thread William Dunlap via R-help
th Jeff’s comments. In fact, I learned most of them the hard way > by trial and error and realized that it’s difficult to separate time and > date using POSIXct and POSIXlt. > > > > Thanks againEK > > [[alternative HTML version deleted]] > > _____

Re: [R] create an empty data frame and then fill in it (and then evaluate the mean of semi-hourly data for each day)

2016-06-10 Thread William Dunlap via R-help
on(x) > round(mean(df_snow$snow[df_snow$day == x], na.rm=T))) > > This does not work. I do not understand why the class of df_snow_day$snow > is of type list either: > >day snow > NA NULL > NA.1 NULL > NA.2 NULL > > Where is my mistake? &

[R] Help: How to Convert Binary Data into Text Using R

2016-06-12 Thread Fahman Khan via R-help
on deleted]] ______ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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.

[R] PSM With Panel Data

2016-06-12 Thread Ahmed Raza via R-help
that how I can do this. Any help in this regard will be much appreciated With Best Regards,  Ahmed Arif  [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman

[R] Binary Value into Text

2016-06-13 Thread Fahman Khan via R-help
Any help would be appreciated. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-

[R] Dashed/dotted ecdf plots using lots of points

2016-06-13 Thread A A via R-help
short sample of code that illustrates my issue: #change to 1:100 to see the dotted line x = 1:1 e = ecdf(x) plot(e,do.p=FALSE, lty = "dotted") Is there a way to solve this problem? [[alternative HTML version deleted]] ______ R-help@r

Re: [R] Installation of package "rio" broken

2016-06-14 Thread William Dunlap via R-help
/win-library/3.2/feather' > Warning in install.packages : > running command '"C:/PROGRA~1/R/R-32~1.2/bin/x64/R" CMD INSTALL -l > "C:\Users\admin\Documents\R\win-library\3.2" > C:\Users\admin\AppData\Local\Temp\RtmpwNaGYz/downloaded_packages/feather_0.0.1.tar.gz&#

[R] Error in twitterR package

2016-06-14 Thread Seema Shelke via R-help
;-"xxx" api_secret  <- "xxx" access_token  <- " xxx" access_token_secret  <- "xxx" setup_twitter_oauth(api_key,api_secret,access_token, access_token_secret) R version : Ri386 3.3.0R studio:  0.99.902  Thanks,Seema [[alternative HTML versi

Re: [R] Closing FTP sessions with RCurl

2016-06-15 Thread William Dunlap via R-help
t;http://www.omegahat.net/RCurl/index.html";)) { > >> > >>curl = getCurlHandle() > >> > >>getURL("http://www.omegahat.net/RCurl/index.html";, curl = curl) > >> > >>#getCurlInfo(curl) # I skipped this step > >> > >>rm(curl) #

Re: [R] Closing FTP sessions with RCurl

2016-06-15 Thread William Dunlap via R-help
ship between handles and sessions, so I am probably missing >> something. >> >> >> >> >> >> >> >> #EXAMPLE from getCurl(), p. 39 >> >> >> >> if(url.exists("http://www.omegahat.net/RCurl/index.html";)) { >> &

[R] cv.glm problem

2016-06-15 Thread Amit Sengupta via R-help
outside allowable range") :   missing value where TRUE/FALSE needed Regards,Amit [[alternative HTML version deleted]] ______ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE d

[R] cv.binary problem

2016-06-15 Thread Amit Sengupta via R-help
appreciated. > cv.binary(fullmod) Error in sample.int(x, size, replace, prob) : invalid 'size' argument Regards,Amit [[alternative HTML version deleted]] ______ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.et

Re: [R] extracting coefficients from ar() output

2016-06-16 Thread William Dunlap via R-help
umeric,15 Numeric,36 "Yule-Walker" "data" > > > I get the table as shown above if I use head(). > > How can I extract the ar coefficients from this table? I have already > tried coef() and rollingarma$ar but both do not work. > What can I do? > > Tha

Re: [R] Writing R package that call Fortran codes

2016-06-16 Thread William Dunlap via R-help
; additional to useDynLib(NEpidemic). After doing that I couldn't build the > > package and it gave me another error: > > > > Error in library.dynam(lib, package, package.lib) : > > shared object ‘random_epi.so’ not found > > Error: loading failed > > Execu

Re: [R] what is the best way to process the following data?

2016-06-17 Thread William Dunlap via R-help
6/2016|02:11:15|Update > > |06/16/2016|02:11:15|Close > > |06/16/2016|02:11:15|BOP processing for 0 items has finished > > |06/16/2016|02:11:15|Step 007 > |06/16/2016|02:11:15|Initialization > > |06/16/2016|02:11:15|Filters > > |06/16/2016|02:11:17|Split Items &g

Re: [R] Problems with Mann-Kendall trend test

2016-06-17 Thread Morteza Firouzi via R-help
help.     [[alternative HTML version deleted]] ______ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide comme

[R] row.names c(NA, -length) vs. c(NA, +length)

2016-06-17 Thread William Dunlap via R-help
Bill Dunlap TIBCO Software wdunlap tibco.com [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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.

[R] Merging Issue

2016-06-17 Thread Farnoosh Sheikhi via R-help
;, "C", "C", "D", "A", "F", "G", "A", "F", "A", "D")df <- data.frame(Subject, dates, deps) ## Data Bloc<-c("CA","NY", "CA", "NY", "WA", "

Re: [R] Extract character strings from a vector

2016-06-21 Thread William Dunlap via R-help
","RTG", "GF TYHH", "HJN 89963", "KFTR1","RT 8") > > > > >From this vector, how can I extract the following character strings > (i.e., which contain 0 or 1 numeric value) > > > > [1] "RTG" &quo

Re: [R] issues with R raster temporary files

2016-06-22 Thread William Dunlap via R-help
just starts > with the date), but I'm not sure why this is happening or how to fix it. > > Any advice? > > Thanks, > Alex > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz

[R] Saving .eps files with Times New Roman font

2016-06-23 Thread A A via R-help
ver, when I right click the figure and select 'Save as postscript...', I get a .eps file of the plot, but without the font formatting that I applied. Is there any way to apply the formatting so that it gets saved? [[alternative HTML version deleted]] _

Re: [R] Subscripting problem with is.na()

2016-06-24 Thread William Dunlap via R-help
es, but you will have to do it > explicitly yourself if you code to 0. That seems to be asking for > trouble to me. > > As always, contrary views welcome. This discussion still seems on > (r-help) topic to me, but if not, please say so. > > Cheers, > Bert > > Bert Gunter &

Re: [R] What's box() (exactly) doing?

2016-06-24 Thread William Dunlap via R-help
't figure out how to solve the problem. Below is a minimal > >>> working (and a slightly bigger) example. > >>> > >>> Cheers, > >>> Marius > >>> > >>> ## MWE > >>> plot(NA, type = "n", ann = FALSE, axes = FALSE, xlim

[R] How I can calculate the value of response variable

2016-06-26 Thread rezvan hatami via R-help
 How I can calculate the value of response variable in a linear model of a matrix of several variables?Can somebody please answer me? Cheers Rezvan Hatami [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To

[R] Fw: How I can calculate the value of response variable

2016-06-26 Thread rezvan hatami via R-help
rom: David Winsemius To: rezvan hatami Cc: "r-help@r-project.org" Sent: Monday, 27 June 2016, 1:20 Subject: Re: [R] How I can calculate the value of response variable > On Jun 26, 2016, at 3:23 AM, rezvan hatami via R-help > wrote: > >  How I can calculate the

Re: [R] How I can calculate the value of response variable

2016-06-26 Thread rezvan hatami via R-help
what function I need to use for this. Thank you for calculating the rain, but it wasn't a part of my question. Cheers From: Jim Lemon To: rezvan hatami ; r-help mailing list Sent: Monday, 27 June 2016, 10:28 Subject: Re: [R] How I can calculate the value of response variable

Re: [R] Fw: How I can calculate the value of response variable

2016-06-26 Thread rezvan hatami via R-help
the right thing. Cheers From: Bert Gunter To: rezvan hatami Cc: "r-help@r-project.org" Sent: Monday, 27 June 2016, 10:34 Subject: Re: [R] Fw: How I can calculate the value of response variable Is this homework? This list tries to enforce  a no homework policy... If not, i

Re: [R] Fw: How I can calculate the value of response variable

2016-06-26 Thread rezvan hatami via R-help
ter now? From: Bert Gunter To: rezvan hatami Cc: "r-help@r-project.org" Sent: Monday, 27 June 2016, 10:34 Subject: Re: [R] Fw: How I can calculate the value of response variable Is this homework? This list tries to enforce  a no homework policy... If not, it looks to me a

Re: [R] How I can calculate the value of response variable

2016-06-26 Thread rezvan hatami via R-help
ter now? From: Jim Lemon To: rezvan hatami ; r-help mailing list Sent: Monday, 27 June 2016, 10:28 Subject: Re: [R] How I can calculate the value of response variable Hi Rezvan,This looks like a simple problem of solving linear equations (and very much like a homework exerc

[R] Shorting Matrices

2016-06-28 Thread Vasilis Bardakos via R-help
f data",22)] Is this even correct? Also, for some reason when I set var = Matrix[,2] somehow the values changed from rationals to integers. Thank you in advance [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To

Re: [R] t-test for regression estimate

2016-06-28 Thread BONACHE Adrien via R-help
, copy paste the code in a new function and save your code. Regards, Adrien. De : Steven Yen À : "Fox, John" Cc : R-help Envoyé le : Mercredi 29 juin 2016 0h43 Objet : Re: [R] t-test for regression estimate Thanks John. Reason is I am doing linear transformatio

Re: [R] [FORGED] Splitting data.frame into a list of small data.frames given indices

2016-06-29 Thread William Dunlap via R-help
;user system elapsed > >> 43.22 74.72 118.09 > > > > > > I'm not sure that I follow what you're doing, and your example is not > > reproducible, since we have no idea what "peaks" is, but on a toy example > > with 5e6 rows in the d

[R] Bestglm subset analysis

2016-06-29 Thread D Wolf via R-help
oug [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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.

[R] Extracting matrix from netCDF file using ncdf4 package

2016-07-02 Thread Hemant Chowdhary via R-help
ind(dX,idX)} } Is there more elegant/faster way other than to using a For Loop like this? It seems very slow when I may have to get much larger matrix where nD can be more than 1000. Thank you HC [[alternative HTML version deleted]] ______ R-help@r-p

[R] BCa Bootstrapped regression coefficients from lmrob function not working

2016-07-03 Thread varin sacha via R-help
;,index=2) Any help would be highly appreciated, S ______ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and pr

Re: [R] Extracting matrix from netCDF file using ncdf4 package

2016-07-03 Thread Hemant Chowdhary via R-help
  ncvar_get(nc=myNC, varid=“myVar”) and then use R’s subsetting abilities. You can do fancier subsetting of arrays in memory than you can to arrays on disk. HTH, -Roy > On Jul 2, 2016, at 3:43 PM, Hemant Chowdhary via R-help > wrote: > > I am working with a 3-dimensional netCDF fi

Re: [R] Extracting matrix from netCDF file using ncdf4 package

2016-07-04 Thread Hemant Chowdhary via R-help
things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Sat, Jul 2, 2016 at 3:43 PM, Hemant Chowdhary via R-help wrote: >  I am working with a 3-dimensional netCDF file having dimensions of X=100, >Y=200, T=365. > My objective is to extract t

Re: [R] BCa Bootstrapped regression coefficients from lmrob function not working

2016-07-05 Thread varin sacha via R-help
REDONDANCE + competitivite + innovation) boot.ci(results, type= "bca",index=2) Best, S De : peter dalgaard Cc : R-help Mailing List Envoyé le : Dimanche 3 juillet 2016 18h19 Objet : Re: [R] BCa Bootstrapped regression coefficients from lmrob function no

Re: [R] BCa Bootstrapped regression coefficients from lmrob function not working

2016-07-06 Thread varin sacha via R-help
r n'est pas multiple de la taille du remplacement In English it would be something like : number of items to replace is not a multiple of replacement length Best S De : Bert Gunter Cc : peter dalgaard ; R-help Mailing List Envoyé le : Mercredi 6 juillet

Re: [R] BCa Bootstrapped regression coefficients from lmrob function not working

2016-07-06 Thread varin sacha via R-help
calculate ? Best S - Mail original - De : peter dalgaard À : varin sacha Cc : Bert Gunter ; R-help Mailing List Envoyé le : Mercredi 6 juillet 2016 11h05 Objet : Re: [R] BCa Bootstrapped regression coefficients from lmrob function not working Offhand, I would suspect that the cause is

Re: [R] BCa Bootstrapped regression coefficients from lmrob function not working

2016-07-06 Thread varin sacha via R-help
Dear Duncan, Many thanks for your reply. What you propose is a very good idea yes. However, I am interested in the coefficients... Best, S De : Duncan Murdoch Cc : R-help Mailing List Envoyé le : Mercredi 6 juillet 2016 14h13 Objet : Re: [R] BCa

[R] R code for multilevel latent class analysis

2016-07-07 Thread Yahoo Mail via R-help
missed some passages? Thank you very much for your help! Cristina Cametti ______ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guid

[R] Revolutions blog: June 2016 roundup

2016-07-07 Thread David Smith via R-help
Since 2008, Microsoft (formerly Revolution Analytics) staff and guests have written about R every weekday at the Revolutions blog: http://blog.revolutionanalytics.com and every month I post a summary of articles from the previous month of particular interest to readers of r-help. And in case

Re: [R] 'ref' must be an existing level

2016-07-09 Thread William Dunlap via R-help
whether a > >>> person will go for a diploma, general or honors. > >>> > >>> The code below: > >>> m11$prog2<- relevel(m11$prog, ref = "honors" > >>> > >>> already loaded the nnet library. However i got the below error: > >>>

Re: [R] How to make the "apply" faster

2016-07-10 Thread William Dunlap via R-help
ster? > >> > > > > If dim(x)[3] << prod(dim(x)[-3]), > > > > new <- Reduce("+",lapply(1:dim(x)[3],function(z) x[,,z,]>=70)) > > > > will be faster. > > > > However, if you can follow Peter Langfelder's suggestion to use rowSums, > >

Re: [R] How to make the "apply" faster

2016-07-11 Thread William Dunlap via R-help
") {b <- (data <= threshold)} else {b <- (data >= >>> threshold)}b[b==TRUE] = 1 y <-rle(b) ans >>> <-length(subset((y$lengths[y$values==1]), (y$lengths[y$values==1])>=2)) >>> return(ans)}* >>> >>> Do you have any idea how

Re: [R] How to make the "apply" faster

2016-07-11 Thread William Dunlap via R-help
gt;>> >>>> >>>> >>>> >>>> >>>> >>>> >>>> *## Calculate spell durationfun.spell.deb <- function(data, threshold = >>>> 1, >>>> direction = c("above", "below")){ #co

Re: [R] use value in variable to be name of another variable

2016-07-11 Thread William Dunlap via R-help
; "AT1G01030" >>> And then, >>> AT1G01030 <- tTargTFS[-1,2] >>> >>> I want to do this up to tTargTFS[1, 2666], so I want to do this in a >>> script >>> and not manually. >>> tTargTFS is a list of 2: chr [1:265, 1:2666], but I a

Re: [R] Error:'subscript out of bounds'

2016-07-12 Thread William Dunlap via R-help
ionodd) : > object 'coordinationodd' not found > > > Could you please tell me how I can fix this problem? > > Best, > Samsad > [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https

[R] About "SpatioTemporal" Package

2016-07-13 Thread Elham Daadmehr via R-help
a(y1, datst1) :   Unable to find column 'obs$date', using 'obs[,2] 3: In createSTdata(y1, datst1) :   Unable to find column 'obs$ID', using 'obs[,3] I checked the example data "mesa.data.raw", I found that my "ID" has the same type (in

Re: [R] About "SpatioTemporal" Package

2016-07-13 Thread Elham Daadmehr via R-help
aste("ID",60101:60128,sep="") and see what happens. Jim On Wed, Jul 13, 2016 at 8:29 PM, Elham Daadmehr via R-help wrote: >  Dear all, > >  I want to build the data with the structure of STdata (using the >  function "createSTdata()"), I did as the follow

Re: [R] Forking and adapting an R package

2016-07-14 Thread BONACHE Adrien via R-help
function. After it copy and paste the function in R. Use it! Best wishes with your code! Adrien. De : "t...@timogrossenbacher.ch" À : "Adams, Jean" ; r-help Envoyé le : Jeudi 14 juillet 2016 8h40 Objet : Re: [R] Forking and adapting an R package On Jul 13, 201

Re: [R] failure with merge

2016-07-14 Thread William Dunlap via R-help
0.9815738 0.03881885 > > > Browse[3]> x[, seq_len(l.b), drop = FALSE] > select method > 2 TRUE GCV.Cp > 1 FALSE GCV.Cp > > and this line executes: > > Browse[3]> order(x[, seq_len(l.b), drop = FALSE]) > [1] 1 2 3 4 > > although nrow(x

Re: [R] Selecting 1st and last dates from a set of dates

2016-07-14 Thread William Dunlap via R-help
9 7 8 > 03/02/1997 9 6 8 > > > I want to select the data on the first date it can be 1st , 2nd or 3rd or > any and last date it can be 31st, 30th and /or29th. I don’t need time. > It would be great if you could help. > Regards, > Gaurang Meht

Re: [R] Selecting 1st and last dates from a set of dates

2016-07-14 Thread William Dunlap via R-help
ices & Enterprise Support Administration > Washington State Department of Social and Health Services > > > > -Original Message- > > From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Mehta, > > Gaurang > > Sent: Thursday, July 14, 2016 9:01 AM

Re: [R] How to group by and get distinct rows of of grouped rows based on certain criteria

2016-07-14 Thread William Dunlap via R-help
t; > > *Data shown below* > |ATP Group|Business Event|Category| > |02 |A |AC | > |02 |A |AD | > |02 |A |EQ | > |ZM |A |AU | > |ZM |A |AV | > |ZM |A

Re: [R] Please help: cannot import files from Windows into R

2016-07-14 Thread William Dunlap via R-help
a language of their own - some kind of old Irish. From "Enchanted > Woods", **William Butler Yeats* > > > > *It is in stone that the most abundant testimony is found of craftsmanship > in medieval Ireland. From Heritage of Ireland, **Brian de Breffny* > >

Re: [R] Selecting 1st and last dates from a set of dates

2016-07-15 Thread William Dunlap via R-help
t > not quite. > Regards, > Gaurang Mehta > -Original Message- > From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of William > Dunlap via R-help > Sent: 14 July 2016 17:32 > To: Nordlund, Dan (DSHS/RDA) > Cc: R-help Mailing List > Subject: Re: [R]

[R] R2WinBUGS with Multivariate Logistic Regression

2016-07-16 Thread Christopher Kelvin via R-help
"or.pam", "or.pah"), model1,n.iter = 11000, n.burnin = 1000, n.chains = 2, codaPkg = T, DIC = TRUE) ### Posterior summaries line.coda <- read.bugs(lineout1) summary(line.coda) __ R-help@r-project.org mailing list -- To UNS

Re: [R] means by year, month and day

2016-07-18 Thread William Dunlap via R-help
p working on it. > > > > Thank you for your kind response and examples. I will study them. > > > > Sincerely, Tom > > > > > > > From: jim holtman [jholt...@gmail.com] > Sent: Sunday, July 17, 2016 7:14 PM > To: Jianling

Re: [R] Concatenate two lists replacing elements with the same name.

2016-07-19 Thread William Dunlap via R-help
rns: > > u$a=6 > u$b=8 > u$c=3 > u$d=5 > u$e=7 > > I.e. it combines the 3 lists, but when names have the same value it > keeps the most recent one. > > Does such a function exists? > > Thanks for the help, > > Cheers, > Luca > > __

Re: [R] pairs: adjusting margins and labeling axes

2016-07-20 Thread William Michels via R-help
xt <- paste("r= ", txt, sep = "") > > > if(missing(cex.cor)) cex.cor <- 0.8/strwidth(txt) > > > text(0.5, 0.6, txt, cex=cex.cor * r) > > > > > > # p-value calculation > > > p <- cor.test(x, y)$p.value > > >

Re: [R] how to expand the dataframe

2016-07-20 Thread William Dunlap via R-help
> > > > > > > [[alternative HTML version deleted]] > > > > Really isn't it time you learned how to send plain text. You've > > posted many questions on Rhelp. It's really not that difficult on > gmail. I > > also have a gmail accou

Re: [R] how to expand the dataframe

2016-07-21 Thread William Dunlap via R-help
1990-02-14 > >>>> ... > >>>> > >>> > >>> Make a dataframe with a 'time' column using seq.Date and merge that > >>> dataframe with your df dataframe. > >>> > >>> > >>>>

Re: [R] about interpolating data in r

2016-07-21 Thread William Dunlap via R-help
A B C > 2009-01-013 4.5 > 2009-01-024 5 > 2009-01-033.3 6 > 2009-01-044.1 7 > 2009-01-054.4 6.2 5.4 > ... > > 2009-11-205.1 5.5 6.1 > 2009-11-215.4 4 > ... > 2009-12-314.5 6 > >

[R] Why the order of parameters in a logistic regression affects results significantly?

2016-07-21 Thread Qinghua He via R-help
how to obtain OR for all factors? All factors are categorical variables (i.e., 0 or 1). Thanks! Peter [[alternative HTML version deleted]] ______ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-

Re: [R] about interpolating data in r

2016-07-22 Thread William Dunlap via R-help
se to > end-points. > > Otherwise, you might find yourself in a huge mistake. > > > > As I mentioned in my first sentence, If you have a relationship between > > all columns or you have data for column C for other years (for instance, > > assume that you have data for column C for 20

[R] Firth logistic regression and odds ratio

2016-07-24 Thread Qinghua He via R-help
[[alternative HTML version deleted]] ______ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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,

[R] change the colour line in gamm4 plotting

2016-07-26 Thread Maria Lathouri via R-help
which is the default to another colour; I tried to use the col=" " function, at least for the main line, in the plot command but it is not working. Any advice/suggestion is welcome. Many thanks.Maria [[alternative HTML version deleted]] ______

Re: [R] Date Time in R

2016-07-26 Thread William Dunlap via R-help
names > > > > week_names <- wday(a1$date, label=TRUE) > > > > > > > > > > > > *From:* Shivi Bhatia [mailto:shivipm...@gmail.com] > > *Sent:* July 26, 2016 12:16 PM > > *To:* Tom Wright > > *Cc:* David L Carlson ; r-help >

Re: [R] Date Time in R

2016-07-26 Thread William Dunlap via R-help
ct so far. >>>> > >>>> > >>>> > >>>> > # convert strings (or factors) to a date object >>>> > >>>> > a1$date <- as.Date(a1$date, “%m-%d-%y”) >>>> > >>>> > >>>> > >>>> >

Re: [R] Reduce woes

2016-07-28 Thread William Dunlap via R-help
1,2,3,4,5), 96, > > accumulate=TRUE) > > [[1]] > > [1] 96 > > > > [[2]] > > [1] 96 1 > > > > [[3]] > > [1] 96 1 2 > > > > [[4]] > > [1] 96 1 2 3 > > > > [[5]] > > [1] 96 1 2 3 4 > > > > [[6]]

[R] store result of loop in df

2016-07-29 Thread Alain D. via R-help
VAR<-(binom.test(df[i,1],df[i,2],0.065))$estimate[[1]] } but bin does only contain the last result. What is wrong with my code? Can anyone help? Thank you in advance. Alain ______ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.e

Re: [R] Reduce woes

2016-07-29 Thread William Dunlap via R-help
gt;>> >> call > >>> >> > > to reduce() with a custom accumulator. Here's an example in > >>> Elixir: > >>> >> > > > >>> >> > > map = %{"one" => [1, 1], "three" => [3], "two" => [2, 2]} > >>> >&

Re: [R] about netcdf files

2016-08-02 Thread Marc Girondot via R-help
. Sincerely Marc __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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] Regression expression to delete one or more spaces at end of string

2016-08-02 Thread William Dunlap via R-help
lso, without the $ and ‘space' instead of ‘blank') without success — to > my surprise, characters other than the trailing space were deleted but the > trailing space remained. > > Guidance on the correct syntax would be appreciated. > > Dennis > > Dennis F

[R] Martingale residuals warning

2016-08-03 Thread Isaac Singini via R-help
y = 1:2, xlab = "Cox-Snell Residuals", ylab = "Survival Probability",      main = "Survival Function of Cox-Snell Residuals")curve(exp(-x), from = 0, to = max(surv_full_constr$var6mon_constr), add = TRUE, col = "red", lwd = 2) I will appreciate any assis

Re: [R] Fwd: only plot borders of a region in a scatter plot

2016-08-04 Thread William Dunlap via R-help
s > created by ferret, but I prefer to implement it by R. Anybody know how to > do it? Thanks a lot. > > Cheers, > > Zun Yin > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.eth

[R] difftime in years

2016-08-04 Thread Thomas Subia via R-help
10.10563 days How does one not output days? ______ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide com

Re: [R] difftime in years

2016-08-04 Thread William Dunlap via R-help
one of those places. Bill Dunlap TIBCO Software wdunlap tibco.com On Thu, Aug 4, 2016 at 11:05 AM, Thomas Subia via R-help < r-help@r-project.org> wrote: > Colleagues, > > age_days <- difftime(Date,DOM,units="days") > date_vals$age_yrs <- age_days/365.242 > >

[R] SAS file

2016-08-05 Thread Yuan Jian via R-help
gure it out? thanksyuan __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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] changing x and y ranges in a PCA plot created by library(labdsv)

2016-08-07 Thread William Dunlap via R-help
plot(o.pca, xlim=c(-2, 3), ylim=c(-2, 4))", nothing will > change. I really appreciate if somebody can help on this. > > Cheers > > [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list -- To

[R] Extracting dates to create a new variable

2016-08-08 Thread Farnoosh Sheikhi via R-help
3",  "2011-01-03", "2011-01-05", "2011-01-05", "2011-01-05" , "2011-01-08", "2011-01-08", "2011-01-08", "2011-01-11", "2011-01-11", "2011-01-11")final<-data.frame(Subject, dates, deps, newdate)final I really appreciate any help.  Best,Farnoosh [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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.

[R] Revolutions blog: July 2016 roundup

2016-08-08 Thread David Smith via R-help
Since 2008, Microsoft (formerly Revolution Analytics) staff and guests have written about R every weekday at the Revolutions blog: http://blog.revolutionanalytics.com and every month I post a summary of articles from the previous month of particular interest to readers of r-help. And in case

Re: [R] funnel plot asymmetry

2016-08-10 Thread BONACHE Adrien via R-help
a {meta} R Documentation Test for funnel plot asymmetry Description Test for funnel plot asymmetry, based on rank correlation or linearregression method. Usage | | | | Afficher sur artax.karlin.mff.cuni.cz | Aperçu par Yahoo | | | |   | Kind regards, Adrien. De : christos tb À 

Re: [R] Continuation-parsing / trampoline / infinite recursion problem

2016-08-10 Thread William Dunlap via R-help
Thomas Mailund wrote: > >> > >> make_thunk <- function(f, ...) f(...) > > > > Doh! It is of course this one: > > > > make_thunk <- function(f, ...) function() f(…) > > > > It just binds a function call into a thunk so I can delay its evalua

[R] Extracting dates to create a new variable

2016-08-10 Thread Farnoosh Sheikhi via R-help
lt;- data.frame(Subject, dates, deps)df The final data set should look like this:newdate<-c(" 2011-01-03",   "2011-01-03",  "2011-01-03", "2011-01-05", "2011-01-05", "2011-01-05" , "2011-01-08", "2011-01-08", "201

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