search for: biomed

I am fitting a logistic model to binary data. The response variable is a factor (0 or 1) and all predictors are continuous variables. The main predictor is LT (I expect a logistic relation between LT and the probability of being mature) and the other are variables I expect to modify this relation. I want to test if all predictors contribute significantly for the fit or not I fit the full

Hi All. I am a new R user. Trying to do scatterplot. Not sure how to resolve this error message A<-subset (ErablesGatineau, station=="A") > B<-subset (ErablesGatineau, station=="B") > > plot(diam ~ biom) > abline(lm(diam ~ biom), col = "red") > > goodcases <- !(is.na(diam) | is.na(biom)) > lines(lowess(diam[goodcases] ~

...cannot find a referee-proof citation for this. Probably SAS does it (I know what Bill Venables will say about this argument...). Can anybody provide me a citation on r-square for/against non-linear fits, or tell me that I am wrong? Dieter --------------------------------------- Dr. Dieter Menne Biomed Software 72074 T?bingen Tel (49) (7071) 52176 FAX (49) (7071) 55 10 46 dieter.menne at menne-biomed.de www.menne-biomed.de -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info&...

Hello, I have this command: x.axis <- seq(from=0.5, to=4.5, length.out=13112) How can I which of the x.axis components is the closest to a given value, for example 3.2? Best, Dani -- Daniel Valverde Saub? Grup de Biologia Molecular de Llevats Facultat de Veterin?ria de la Universitat Aut?noma de Barcelona Edifici V, Campus UAB 08193 Cerdanyola del Vall?s- SPAIN Centro de Investigaci?n

Hello, I am creating a plot and I would like to know how to put this expression to the y axis ?mol/10^6 cells I've tried some combinations using the expression() function, but none of them worked. Any idea? Best, Dani -- Daniel Valverde Saub? Grup de Biologia Molecular de Llevats Facultat de Veterin?ria de la Universitat Aut?noma de

Hello, I have a vector, lets say x <- 1:50 I would like it to be cut at certain points, being for example 1:5, 6:11, 12:17, ... How can I do it? I have tried the cut() function, but I don not know how to place the cutting points properly. Best regards, Dani -- Daniel Valverde Saub? Grup de Biologia Molecular de Llevats Facultat de Veterin?ria de la Universitat Aut?noma de Barcelona

Hi, I was just wondering if there are any books on R that have applications using Hidden Markov models or anyone with tutorials/exercises that can be shared. Ideally, these would have a bio/ecological slant though not neccesarily. Any suggestions would be appreciated. Thanks! [[alternative HTML version deleted]]

...l[1])<-"newempl" # no error message, yet .... names(Empl) #[1] "empl" "spo" # where's my newempl??????? If this line names(Empl[1])<-"newempl" was stupid, why is there no error message? --------------------------------------- Dr. Dieter Menne Biomed Software 72074 T?bingen Tel (49) (7071) 52176 FAX (49) (7071) 55 10 46 dieter.menne at menne-biomed.de www.menne-biomed.de -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info&quo...

Forget my last message, it was a really stupid problem. Sorry about the mess. Best, Dani -------- Missatge original -------- Assumpte: Problems when checking a package Data: Tue, 05 Feb 2008 11:20:28 +0100 De: Dani Valverde <daniel.valverde at uab.cat> A: R Help <r-help at r-project.org> Hello, I am tryining to build a package called NMRTools and when I check it using R CMD

Hello, I have two vectors, one with 13112 points and the other one with 10909. I wonder if there is a way to interpolate the data so the shorter vectors has the same number of points as the longer one. Best, Dani -- Daniel Valverde Saub? Grup de Biologia Molecular de Llevats Facultat de Veterin?ria de la Universitat Aut?noma de Barcelona Edifici V, Campus UAB 08193 Cerdanyola del Vall?s-

Hello, I have carried out an lda analysis using the lda function of MASS package. I have plotted the LD1xLD2 to represent the data. Now I would like to get the centroids for each group of data and plot it on the LD1xLD2 graph. How can I get the centroid value from the lda object? Best, Dani -- Daniel Valverde Saub? Grup de Biologia Molecular de Llevats Facultat de Veterin?ria de la

Hello, I would like to combine different label sizes in the same x-axis, for example 1 2 3 4 5 6 7 8 9 How can I do it? Best, Dani -- Daniel Valverde Saub? Grup de Biologia Molecular de Llevats Facultat de Veterin?ria de la Universitat Aut?noma de Barcelona Edifici V, Campus UAB 08193 Cerdanyola del Vall?s- SPAIN Centro de Investigaci?n Biom?dica en Red en Bioingenier?a, Biomateriales y

Hello! I have read somewhere (somehow, I can't seem to find it again, it's been a couple of months) that when analyzing factorial block design, the position where you put the block factor is important, even more when there are missing values. I understand that when using anova.lm, the order is sequential, so that if I want to check for a treatment effect, I should put my blocking factor

Hello, I have the next function call: lme(fixed=Error ~ Temperature * Tumour ,random = ~1|ID, data=error_DB) which returns an lme object. I am interested on carrying out some kind of lsmeans on the data returned, but I cannot find any function to do this in R. I'have seen the effect() function, but it does not work with lme objects. Any idea? Best, Dani -- Daniel Valverde Saub? Grup

Hello, I would like to convert an Rd object to a latex file, so that I can put it in my thesis. How can I do it? I tryed latex(), but it only works for code... Best, Dani -- Daniel Valverde Saub? Grup de Biologia Molecular de Llevats Facultat de Veterin?ria de la Universitat Aut?noma de Barcelona Edifici V, Campus UAB 08193 Cerdanyola del Vall?s- SPAIN Centro de Investigaci?n Biom?dica en

...of by index? The following gives the idea remove<-c("subj","drug") data[-remove] I found a solution with a few evals and substitutes, similar to that used in reshapeLong, but there must be an easier way out. Dieter --------------------------------------- Dr. Dieter Menne Biomed Software 72074 T?bingen Tel (49) (7071) 52176 FAX (49) (7071) 55 10 46 dieter.menne at menne-biomed.de www.menne-biomed.de -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info&quo...

Hello, I would like to place some text that should appear in the same position of the graphic device, preferably outside the plotting area, regardless the x and y axes limits. How can I do this? Best, Dani -- Daniel Valverde Saub? Grup de Biologia Molecular de Llevats Facultat de Veterin?ria de la Universitat Aut?noma de Barcelona Edifici V, Campus UAB 08193 Cerdanyola del Vall?s- SPAIN

...+Asym/(1+exp((xmid-x)/scal)),"x") expression({ .expr3 <- exp((xmid - x)/scal) ## Coorrect > deriv(~offv+Asym/(1+exp(-1.0*(x-xmid)/scal)),"x") expression({ .expr5 <- exp(-1 * (x - xmid)/scal) ## Correct --------------------------------------- Dr. Dieter Menne Biomed Software 72074 Tübingen Tel (49) (7071) 52176 FAX (49) (7071) 55 10 46 dieter.menne@menne-biomed.de www.menne-biomed.de -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info&quot...

...0) # ok demo.lm2<-lme(fma,data=demo,random=~1|subj) # looks ok, but isn't print(summary(demo.lm2)) #Fixed effects: Error in as.vector(x, "list") : cannot coerce to vector #predict(demo.lm2,newframe,level=0) # does not work --------------------------------------- Dr. Dieter Menne Biomed Software 72074 T?bingen Tel (49) (7071) 52176 FAX (49) (7071) 55 10 46 dieter.menne at menne-biomed.de www.menne-biomed.de -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info&quo...

Hello, I am trying to perform a glm analysis on a 68x13113 matrix (named data.spect). The first column corresponds to the predictor (data.spect[,1]) and the rest to the response variables (data.spect[,2:13113]). When I try this code glmObject <- glm(data.spect[,2:13113]~data.spect[,1]) I get the following error: Error: (subscript) logical subscript too long Could anyone help me on