similar to: Generating unque patient IDs

Dear All, I want to create a dataset for a NONMEM simulation. I have a dataframe with individual PK parameters and want to create a dosing sceinario in a second dataframe. I want to merge them both so that every individiual's PK parameters are combined with the dosing scenario into one. I do not have a common variable to merge and cbind is giving error that lengths are different. I

Hi I am having some problems running a covariate analysis with my colleage using R with the NonMem program we are using for a graduate school project. R and NonMem run fine without adding in the covariates, but the program is giving us a problem when the covariate analysis is added. We think the problem is with the R code to run the covariate data analysis. We have the control stream, R code

Hello, I have made an R script that prepares a NONMEM dataset and I would like to start the NONMEM run right after the dataset is ready. I am using windows XP, R 2.9.1 and NONMEM 6. I have prepared a run.bat file that looks like this: ---------------------------------------- call K:\nmvi\NMdirectories.bat call K:\nmvi\nmfe6 "path\control.txt" "path\output.txt"

Dear Group, I am working with a large dataset where I need to select for each unique id the the unique lastpk row. Here is a sample subject: id wtdt wt lastpk 64050256 2010-09-18 275 2010-09-16 64050256 2010-09-19 277 2010-09-18 64050256 2010-09-20 272 2010-09-18 64050256 2010-09-21 277 2010-09-18 I want the result as:

Dear group, I need help on two problems: 1. I am trying to plot density plots for each individual in 8 occasions. I can do this by subject wiht the code below: par(mfrow=c(4,2)) plot(density(all8scenarios$SIMCONC[all8scenarios$ID==1&all8scenarios$WSEQ==0])) plot(density(all8scenarios$SIMCONC[all8scenarios$ID==1&all8scenarios$WSEQ==1]))

hi i am creating some variables from same data, but somewhere is different rouding. look: P = abs(fft(d.zlato)/480)^2 hladane= sort(P,decreasing=T)[1:10]/480 pozicia=c(0,0,0,0,0) for (j in 1:5){ for (i in 2:239){ if (P[i]/480==hladane[2*j-1]){pozicia[j]=i-1}}} period=479/pozicia > P[2]/334 [1] 0.0001279107 > hladane[1] [1] 0.0001279107 > P[2]/334==hladane[1] [1]

Dear R group, I am trying to plot uisng a for loop. Here is the script: for (i in unique(cfn$ID)){ plot(cfn$TIME[cfn$ID==i],cfn$DV[cfn$ID==i],pch=16)} I could access only the last plot of the series and cannot go back to see all the plots. I appreciate your help in resolving this issue and also please suggest how I can export the plots as a pdf or other format files for publication purposes.

Dear Group, I have a problem with time date formatting. I have Date and Time in different columns in a .csv file and want to have a column with Date and Time together. I could format the date into the right mode using strptime and as.Date functions. I am not able to do that with TIME column. I appreciate your help in solving this issue. Currently the TIME column is in character mode.

Hello T I want to order some calculation "result", there will be lots of "result" that need to calculate and order PS: the "result" is just a intermediate varible and ordering them is the very aim # problem: # For fixed NT and CT, and some pair (c,n). order the pair by corresponding result # c and n are both random variable CT<-6000 #assignment to CT

Dear R-people, I have to compute C - -(pnorm(B)*dnorm(B)*B + dnorm(B)^2)/pnorm(B)^2 This expression seems to be converging to -1 if B approaches to -Inf (although I am unable to prove it). R has no problems until B equals around -28 or less, where both numerator and denominator go to 0 and you get NaN. A simple workaround I did was C <- ifelse(B > -25, -(pnorm(B)*dnorm(B)*B

Dear Group, I have the following dataset: ID REPI DV CONC SS 1 1 156.84 116 0 1 2 146.56 116 0 1 3 115.13 116 0 1 4 207.81 116 0 1 5 129.53 116 0 1 6 151.48 116 0 1 7 158.95 116 0 1 8 192.37 116 0 1 9 32.97 116 0 1 10 151.66 116 0 I want to calculate the percentile of each CONC within ID=i and add as a column

Oops. Missed one line in the last patch.... Roger. -- ** R.E.Wolff@BitWizard.nl ** http://www.BitWizard.nl/ ** +31-15-2600998 ** *-- BitWizard writes Linux device drivers for any device you may have! --* * The Worlds Ecosystem is a stable system. Stable systems may experience * * excursions from the stable situation. We are currently in such an * * excursion: The stable situation does

Dear R-users, Do you know if we can use the function lme in R for log-normal distribution of parameters as used in Nonmem ? theta=theta0*exp(eta) In our model, the parameters follow the log-normal distribution so it's not reasonable to deal with normal distribution which gives us negative values in simulation Thanks for your help, Thu [[alternative HTML version deleted]]

Greetings, I used latex type code in my Rd files. The pdf version using R CMD Rd2dvi --output=PKtools.pdf --pdf --title="PKtools" PKtools man came out quite nice. However, my current HTML version does not tex the latex so there is latex code in the files and looks bad. example problem code: AIC$_c$ or \begin{itemize} \item NLME: \begin{itemize} \item population level:

[This email is either empty or too large to be displayed at this time]

Hi Jean-Marc, Alfred and Greg, Are you receiving the mails from IETF about draft-ietf-avt-rtp-speex The mails are not coming from AVT mailing list, but I think we are all 3 part of a minimal list (draft-ietf-avt-rtp-speex at tools.ietf.org) dedicated to latest discussion about the draft. I have answered some questions, but there are small changes and adaptation still required to the ietf

Dear Martin, Sad that the bug is beyond your ken... Fortunately, the error happens only rarely...The length of LYGH was 719 and there were only two such errors..I will just replace them with NA and make do. By the by, what if I send LYGH as an attachment to your actual mail ( not the r-help mail)? Will it help? Can you then pinpoint the cause? Or should I raise a bug

There's some sintyax available to use, with ldb* tools, multiple ''server'', eg use multiple '-H', or write server in a roud-robin style? Thanks. -- dott. Marco Gaiarin GNUPG Key ID: 240A3D66 Associazione ``La Nostra Famiglia'' http://www.lanostrafamiglia.it/ Polo FVG - Via della Bontà, 7 - 33078 - San Vito al Tagliamento (PN)

Hi Everyone: I am running NONMEM from R using following command. shell("nmfe72 TOT_EST_DOSE2.CTL DOSE2.OUT",wait=T) This command runs NONMEM successfully but does not show the progress. Usually, intermediate outputs are displayed on CMD console when this same command is run from CMD console. How can I view the intermediate outputs while running NONMEM from R?

Hello, I have a question about multichannel and the load-balancing algorithm. What kind of algorithm is it? (Roud-robin, dynamic round-robin, ...) What are the arguments the algorithm compute to load-balance the network traffic on the interfaces? I already know that if a 10Gb/s interface and a 1Gb/s interfaces are in the multichannel connection the 1Gb/s one will be ignored until the 10Gb/s goes