search for: abitbol

Dear R Helpers, I need to read time from a .csv file which is formated as chartime (09:12:00) below. I need to add one minute (cf chartime2). Then I need to output the value just as 09:13 without the seconds for writing a csv file and input in another program. I get it with the following reproducible example but I can't help thinking that there must a less clumsy way to do that ! Thanks

...packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] SASxport_1.2.3 Hmisc_3.4-3 foreign_0.8-29 RWinEdt_1.8-0 loaded via a namespace (and not attached): [1] chron_2.3-24 cluster_1.11.11 grid_2.7.1 lattice_0.17-15 Jean-Louis Abitbol, MD Chief Medical Officer Trophos SA, Parc scientifique de Luminy, Case 931 Luminy Biotech Entreprises 13288 Marseille Cedex 9 France Email: jlabitbol at trophos.com ---- Backup Email: abitbol at sent.com Cellular: (33) (0)6 24 47 59 34 Direct Line: (33) (0)4 91 82 82 73-Switchboard: (33) (0)4 91 8...

Dear R Users, (and dear Marc) First of all many thanks for the answers to my previous questions. I would like to barplot the mean percent change of a variate with it's CI. Bars should start from the zero reference line to height (in barplot2). Is there a way to tweak barplot2, for example, to do that ? I have tried to see what the function was but unlike other functions was not able to

Dear All I am trying to graph a proportion and CI95% by a factor with ooplot (any other better solution ?) It works well until I try to add the confidence interval. this is the error message and and a description of the data: > dat1 PointEst TT1 1 3.6 TT2 2 5.0 TT3 3 5.8 TT4 4 11.5 TT5 5 7.5 TT5 6 8.7 TT7 7 17.4 > dat2

Dear all, I am trying to use layout with xyplot and save the plot in a file. The saved file contains only the last panel (used jpg, wmf and pdf options in the savePlot function) I am not familiar with the graphic file formats and don't know if they can handle multipage output. If not do I need to loop by ndose in the following script ? Here is my script:

...der if anyone has a function for doing non compartimental PK analysis by subject on a dataframe with subject, time and concentration variates. Thanks in advance for sharing and helping me not re-invent the wheel (which might be square given my R skills). Best wishes, Jean-Louis Jean-Louis Abitbol, MD Chief Medical Officer Trophos SA, Parc scientifique de Luminy, Case 931 Luminy Biotech Entreprises 13288 Marseille Cedex 9 France Email: jlabitbol at trophos.com ---- Backup Email: abitbol at sent.com Cellular: (33) (0)6 24 47 59 34 Direct Line: (33) (0)4 91 82 82 73-Switchboard: (33) (0)4 91 8...

...rature/scripts on simulation of clinical trials in general or more specific to my problem in R. Please answer to me directly to avoid OT inconvenience to the list and so that I am not black listed...I'll be happy to summarize if relevant to the list. Best regards, Jean-Louis Jean-Louis Abitbol, MD Chief Medical Officer Trophos SA, Parc scientifique de Luminy, Case 931 Luminy Biotech Entreprises 13288 Marseille Cedex 9 France Email: jlabitbol at trophos.com ---- Backup Email: abitbol at sent.com Cellular: (33) (0)6 24 47 59 34 Direct Line: (33) (0)4 91 82 82 73-Switchboard: (33) (0)4 91 8...

...sep=".")) which result in an error message such as "target of the assign does not belong to the language" (translation from French). I have looked in the FAQ but do not understand how to use assign in this context. Thanks for any help, Kind regards, Jean-Louis Dr. Jean-Louis Abitbol Directeur Recherche Clinique Trophos SA, Parc scientifique de Luminy, Case 931 Luminy Biotech Entreprises 13288 Marseille Cedex 9 Tel 04 91 82 82 73, Fax 04 91 82 82 89 Portable 06 24 47 59 34

Dear R Gurus Just started on R ! Using xYplot from Hmisc (R 1.9, W2K) I get a grey/blue background that I would like to change to white (ie no background) or may be to another color. Tried to do that with par(bg) but only changed the color of the trellis heading. What's the right command to do that ? Kind regards, JL PS if anyone has nice default settings for win device please let me

Dear All, I am trying to convert to POSIXct after pasting a date and a time in character format with strptime. It is probably obvious but I don't understand why I get an error message after bsamp$spltime<-strptime(test,format="%d-%B-%y %H:%M") whereas I can get what I want if I do it in 2 steps and rbinding ? Thanks and best regards, Jean-Louis This is the R console output

> From: Jean-Louis Abitbol > > Dear All, > > I would like to merge a data frame such as: > > > basetab > subject dose cmax > 1 1031 50 21.8 > .... > > with the result of a by processing such as: > > > tmax<-by(pkga,subject,f.tmax) > > tmax > subject: 10...

Dear R Gurus, Sorry if this has been asked before but I did not find it in the archives. I would like to add a horizontal display of mean and SEM on Hmisc ecdf plots done by group (ie variate~treatment). Has anyone written some code to do that ? Thanks and kind regards, Jean-Louis

Dear All, I have the following data coming out from s <- with(final, summarize(norm, llist(gtt,fdiab), function(norm) { n <- sum(!is.na(norm)) s <- sum(norm, na.rm=T) binconf(s, n) }, type='matrix') ) ie gtt fdiab norm.norm norm.norm2 norm.norm3 18

...uestion (Gabor Grothendieck) 29. How to remove columns of table in R recursively (Jason Lee) 30. turning comma separated string from multiple choices into flags (June Kim) 31. Script for dose proportionality assessment using power model and confidence interval criteria (Jean-Louis Abitbol) 32. Re: ariable Importance Measure in Package RandomForest (Liaw, Andy) 33. Re: turning comma separated string from multiple choices into flags (Peter Dalgaard) 34. Re: density estimate (Eik Vettorazzi) 35. Re: turning comma separated string from multiple choices into fla...

Dear R Helpers I am trying to get function smedian.hilow to work using Hmisc summarize on variable conc in dataframe pkindivmtd by time and dose using: attach(pkindivmtd) sconc <- summarize(conc,llist(time,dose),smedian.hilow) I get the error message Erreur dans 1:nrow(X) : argument NA / NaN Contents of pkindivmtd is > contents(pkindivmtd) Data frame:pkindivmtd 1296 observations