similar to: covariance matrix under null

Hello R users, I am a fairly new R user, however, one of the problems I am having is the use of applying if-, and-, then-, else- statements in R against datasets & dataframe. Rtips mentions it, but without simple examples I am not quite sure what my syntax should look like. I am particularly interested in syntax that can be used with the transform function. for example if foo[,1] ==

Hello all, I am working with a 1000x1000 matrix, and I would like to return a 1000x1000 matrix that tells me which value in the matrix is greater than a theshold value (1 or 0 indicator). i have tried mat2<-as.matrix(as.numeric(mat1>0.25)) but that returns a 1:100000 matrix. I have also tried for loops, but they are grossly inefficient. THanks for all your help in advance. Lanre

I am trying to make a box plot and overlay it with a scatter plot from another data.frame. I was able to successfully create the boxplot, but when i tried using points(x~y...) the dots did not show up. example code aa<-(300,300,300,300,600,600,600,600,900,900,900,900) bb<-(13,12,14,11,56,44,34,75,22.,34,22,98,59,55,56) cc<-(13,12,14,11,56,44,34,75,22.,34,22,98,59,55,56)

hello all, I have a dataset where the subjects are duplicated. How do I subset such that I can get only I row/subject. aa<-c(1,1,2,2,3,3,4,4,5,5,6,6) bb<-c(56,56,33,33,53,53,20,20,63,63,9,9) cc<-data.frame(aa,bb) I would like to subset df(cc) such that I can get aa bb 1 56 2 33 3 53 4 20 5 63 6 9 I know this should be fairly easy but I can't figure how to do it in a dataframe and

Dear list, I am running some linear regressions through lapply, >lapply(c('EMAX','EC50','KOUT','GAMMA'),function(x)confint(lm(get(x)~RR0,dataset2))) I got results like [[1]] 2.5 % 97.5 % (Intercept) 0.6595789212 0.8821691261 RR0 -0.0001801771 0.0001489083 [[2]] 2.5 % 97.5 % (Intercept) -63.83694930

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 this data: Time AMP 0 0.2000000 10 0.1958350 20 0.2914560 40 0.6763628 60 0.8494534 90 0.9874526 120 1.0477692 where AMP is the concentration of this metabolite with time. If you plot the data, you can see that it could be fitted using a logistic regression. For this purpose, I used this code: AMP.nls <- nls(AMP~SSlogis(Time,Asym, xmid, scal), data

Dear all, I want to do a F test, which involves calculation of the degrees of freedom for the residuals. Now say, I have a nlme object "mod.nlme". I have two questions 1.How do I extract the degrees of freedom? 2.How is this degrees of freedom calculated in an nlme model? Thanks. Jun Shen Some sample code and data =================================================================

Folks I have modified an existing function to calculate 'ec/ld/lc' 50 values and their associated Fieller's confidence limits. It is based on EC50.calc (writtien by John Bailer) - but also borrows from the dose.p (MASS) function. My goal was to make the original EC50.calc function flexible with respect to 1) probability at which to calculate the expected dose, and 2) the link

Colleagues, I am fitting an Emax model using nls. The code is: START <- list(EMAX=INITEMAX, EFFECT=INITEFFECT, C50=INITC50) CONTROL <- list(maxiter=1000, warnOnly=T) #FORMULA <- as.formula(YVAR ~ EMAX - EFFECT * XVAR^GAMMA / (XVAR^GAMMA + C50^GAMMA)) ## alternate version of formula FORMULA <- as.formula(YVAR ~ EMAX - EFFECT / (1 + (C50/XVAR)^GAMMA)) FIT <-

Hello all I have done some fitting of pnorm functions to dose-response data, so I could calculate EC50 values (dose where the response is 0.5). I used the nlm function for this, so I did not get any information about the confidence intervals of the fitted parameters. What would be a good way to do such a probit fit, or is there a package which I could use? Best regards Johannes Ranke

Hello All,  I am working on migrating the vm from redhat to ubuntu. the emulator value is as below.  Ubuntu: <emulator>/usr/bin/qemu-system-x86_64</emulator> redhat: <emulator>/usr/libexec/qemu-kvm</emulator> When i try to migrate from redhat, i get the following error. It would be great help if you can suggest how to workaround the issue.  virsh migrate --live --verbose

Hello, We use drc to fit dose-response curves, recently we discovered that there are quite different standard error values returned for the same dataset depending on the drc-version / R-version that was used (not clear which factor is important) On R 2.9.0 using drc_1.6-3 we get an IC50 of 1.27447 and a standard error on the IC50 of 0.43540 Whereas on R 2.7.0 using drc_1.4-2 the IC50 is

Hi, I am trying to construct a logistic regression model from my data (104 patients and 25 events). I build a full model consisting of five predictors with the use of penalization by rms package (lrm, pentrace etc) because of events per variable issue. Then, I tried to approximate the full model by step-down technique predicting L from all of the componet variables using ordinary least squares

Dear R Users: I have a list function as: Flat: y = 0 Linear: y = -(1.65/8)d Logistic: y = 0.015 - 1.73/{1+exp[1.2(4-d)]} Umbrella: y= -(1.65/3)d + (1.65/36)d^2 Emax: y = -1.81d/(0.79+d) Sigmoid Emax: y = -1.70d^5/(4^5+d^5) And want draw the figure as attached (those material are extracted from a paper). Could anyone give me a sample code to do this? Thanks

Hi Team, After performing the gluster clean up and again doing gluster configuration the sub directories become invisible. Even after they are present the directories are invisible. If I create a new directory by the same name it will say the directory already exists. Thanks & Regards, Sangeeta Ramapure -------------- next part -------------- An HTML attachment was scrubbed... URL:

Classification: UNCLASSIFIED Caveats: NONE A similar question has been posted in the past but never answered. My question is this: for probit analysis, how do you program a 95% confidence interval for the LD50 (or LC50, ec50, etc.), including a heterogeneity factor as written about in "Probit Analysis" by Finney(1971)? The heterogeneity factor comes into play through the chi-squared

can you please describe a bit more about the steps taken to clean up and re-configure gluster? Regards, Vijay On Mon, Jun 12, 2017 at 12:08 PM, Sangeeta Ramapure < sangeeta.ramapure at globallogic.com> wrote: > Kindly somebody help me with this issue. > > > > Thanks & Regards, > > Sangeeta Ramapure > > > > *From:* Sangeeta Ramapure

This is to announce the release of Icecast 2.4 beta5 (aka 2.3.99.5). We're looking forward to testing and feedback. As this is a beta release please discuss it on the icecast-dev mailing list. (Please do NOT cross-post to both lists, do NOT reply to the icecast mailing list, PLEASE reply to icecast-dev at xiph.org only). * Updated web interface to be fully XHTML compliant and mobile

This is to announce the release of Icecast 2.4 beta5 (aka 2.3.99.5). We're looking forward to testing and feedback. As this is a beta release please discuss it on the icecast-dev mailing list. (Please do NOT cross-post to both lists, do NOT reply to the icecast mailing list, PLEASE reply to icecast-dev at xiph.org only). * Updated web interface to be fully XHTML compliant and mobile