similar to: Switch Help

Dear Rexperts, I am trying to add a '+' identifying the mean in a boxplot using the following sizelist <- split(size, grp) centers <- boxplot(sizelist, style.bxp = "att", medpch = "o", ylab = "Prostate Volume (cm3)") points(centers, unlist(lapply(sizelist, mean)), pch = "+") But, I get error Error in xy.coords(x, y) :

Dear Wiza[R]ds, I have a data.frame header that looks like this: v2FfaPre15 v2FfaPre10 v2FfaPre5 v2Ffa2 v2Ffa3 v2Ffa4 I need it to look like this, 15 10 5 2 3 4 i.e., with v2FfaPre and v2Ffa stripped off Any suggestions, Thanks in advance! -- Oscar Oscar A. Linares, MD Translational Medicine Unit LaPlaisance Bay, Bolles Harbor Monroe, Michigan 48161

Dear Wiz[R]ds, I am deeply grateful for the help from Duncan Murdoch, Gray Calhoun, and others. We are almost there. For whatever reason, I can't change the symbol from a circle to a triangle in the upright posture plots. Any ideas? I have included the problem in full. # tritiated (3H)-Norepinephrine(NE) disappearance from plasma # concentrations supine and upright # supine datasu <-

Dear Wiza[R]ds, I have the following data in a data.frame. I need to do an anova with multiple comparison but I don't know how to factor the groups for analysis. There are 3 groups, 1,2 and 3 labelled in column 1. Help appreciated with thanks in advance. Group SI Sif SG Io I2 lol2 1 9.08e-05 9.08e+00 0.060842287 1.798556446 32.574500

Dear Rxperts, Using R 2.8.1 and trying install.views("Cluster") getting error Error: could not find function "install.views" Please help:-( -- Oscar Oscar A. Linares Molecular Medicine Unit Bolles Harbor Monroe, Michigan [[alternative HTML version deleted]]

Is there a UseR! manual in the works for the outstanding Rcmdr? -- Oscar Oscar A. Linares, MD Translational Medicine Unit Bolles Harbor Monroe, Michigan [[alternative HTML version deleted]]

Dear Rxperts, I have a varaibles data file that looks like this p(1) 10 p(1) 3 p(1) 4 p(2) 20 p(2) 30 p(2) 40 p(3) 4 p(3) 1 p(1) 2 I cannot process these data with R because it does not like the parentheses. How can I get these to look like: p1 10 p1 3 p1 4 p2 20 p2 30 p2 40 p3 4 p3 1 p3 2 The data is in a tab delimited text file and I want to get it into a data.frame(). Many thanks in

Dear Rxperts, Running the following code: ======================================================= twlo=10; twhi=20; wt=154; vd=0.5; cl=0.046; tau=6; t=3; F=1; wtkg <- wt/2.2 # convert lbs to kg vd.pt <- wtkg * vd # compute weight-based vd (L) cl.pt <- wtkg * cl # compute CL (L/hr) k <- cl.pt/vd.pt # compute k (hr^-1) cave <-

Dear Rxperts I have a data.frame as follows A B C D 1 4 7 10 2 5 8 11 3 6 9 12 I want to convert it to a data frame with a single row (i.e., stack the columns without the heading) 1 2 3 4 5 6 7 8 9 10 11 12 Any suggestions please. Thanks in advance! -- Oscar Oscar A. Linares Molecular Medicine Unit Bolles Harbor Monroe, Michigan [[alternative HTML

Dear Rxperts, I have the following data: Study Study.Name C Category TC Time QC QO SD FSD Theta 1 NONE 0 P(22) 0 0.00 7.5596 0 0 8.0361e-03 0 1 NONE 6 G(50) 0 0.00 1.0000 0 0 0.0000e+00 0 1 NONE 2 F(02) 0 0.00 100.0000 0 0 0.0000e+00 0 1 NONE 3 F(03) 0 0.00 13.2280 0 0 1.6732e-02 0

Dear Wiza[R]ds, # I have the following experimentally observed data: csdata <- data.frame( time=c(0,1,3,9,20), conc=c(638.697,395.69,199.00,141.58,112.16) ) # weighting resp means response wt.MM<- function(resp, time,A1,a1,A2,a2) { pred <- A1*exp(-a1*time)+A2*exp(-a2*time) (resp - pred) / sqrt(pred) } # Fit using nls cs.wt <- nls( ~ wt.MM(conc, time,A1,a1,A2,a2),

After some processing... ct.df <- data.frame(time,conc) ct.df gives time conc 1 0 164.495456 2 1 133.671185 3 2 108.622975 4 3 88.268468 5 4 71.728126 6 5 58.287225 7 6 47.364971 8 7 38.489403 9 8 31.276998 10 9 25.416103 11 10 20.653462 12 11 16.783276 13 12 13.638313 14

Dear Wiza[R]ds, I am very grateful to Duncan Murdoch for his assistance with this problem. His help was invaluable. However, the problem has become a little more complicated for me. Now, in each plot, I need to plot the observed and fitted values of a supine and upright posture experiment. Here is what I have and how far I got. # tritiated (3H)-Norepinephrine(NE) disappearance from plasma #

Hello, i have a xls file. I will read it in r, what library-command i use for this?? any ideas?? Thanks Michele [[alternative HTML version deleted]]

Hola, Estoy intentando realizar un ajuste a un modelo setar con R y me surgen problemas a la hora de representar el modelo setar sin la constante. El código es el siguiente: # Estimacion TAR(2;2,2) con delta 1: mod.setar <- setar(x, m = 2, mL = 2, mH = 2, thDelay = 1) mod.setar summary(mod.setar) mod.setarc<-setar(x, m=2, mL=2, mH=2, thDelay=1, include=c("none")) ##eliminamos

Dear All, I like to simulate a physiologically based pharmacokinetics model using R but am having a problem with the daspk routine. The same problem has been implemented in Berkeley madonna and Winbugs so that I know that it is working. However, with daspk it is not, and the numbers are everywhere! Please see the following and let me know if I am missing something... Thanks a lot in advance,

Hi, I am using the mgcv package (version 1.7-22.) running the model works fine, but when I want to have a plot with residuals I get an error. fit29<-gam(IV~s(G3)+s(V3)+factor(AAR)+s(D3)+s(RUTE,bs="re"),data=subsf,gamma=1.4,method="ML") plot(fit29,residuals=T) Error in X[, first:last] %*% object$coefficients[first:last] : non-conformable arguments does some one know

While running a multiple VM testscase with each VM running iperf traffic between others the following kernel NULL pointer exception was seen. Race appears when the tun driver instance of one VM calls skb_array_produce (from tun_net_xmit) and the the destined VM's skb_array_consume (from tun_ring_recv), which could run concurrently on another core. Due to which the sock_wfree gets called again

While running a multiple VM testscase with each VM running iperf traffic between others the following kernel NULL pointer exception was seen. Race appears when the tun driver instance of one VM calls skb_array_produce (from tun_net_xmit) and the the destined VM's skb_array_consume (from tun_ring_recv), which could run concurrently on another core. Due to which the sock_wfree gets called again

Hi, I am using the gam function in the mgcv package, I have random effects in my model (bs="re") this has worked fine, but after I updated the mgcv package to version 1.7-19 I recive an error message when I run the model. > fit1<-gam(IV~s(RUTE,bs="re")+s(T13)+s(H40)+factor(AAR)+s(V3)+s(G1)+s(H1)+s(V1)+factor(LEDD),data=data5,method="ML") > summary.gam(fit1)