search for: biosciences

...n('r[adj]=0.40'),expression('F[3,23]=15.43'),expression('P < 0.001')),col=black, bty='n',cex=0.9) Luego al tratar de colocar las cusriva utilizando par(font=3) tampoco funciona. ¿Alguna idea de un script para hacer esto? Gracias Eric Flores De Gracia School of Biosciences University of Exeter, Cornwall Campus Penryn, Cornwall TR10 9EZ United Kingdom Mobilephone: +044 07578724705 http://biosciences.exeter.ac.uk/cec/staff/postgradresearch/index.php?web_id=Eric_Flores_de_Gracia&tab=profile

Hello group, For my PhD I'm working on a spatial sampling grid. I do have two data sets which I'd like to compare using cross-correlograms or cross-variograms. Is this an option in one of the R-packages? I've been searching the R-help archive and the available package-documentations, but I can't find how to do this. Thanks in advance, Ren?.

...ado planteo el script tengo algo como: modelo<-coxp(surv(deaths,day)~size+strata(block)) El factor bloque es random pero no se si esta es la nomenclatura correcta para testar mi hipotesis, o simplemente debo omitir el termino -strata- Gracias,saludos Eric Flores De Gracia School of Biosciences University of Exeter, Cornwall Campus Penryn, Cornwall TR10 9EZ United Kingdom Mobilephone: +044 07578724705 http://biosciences.exeter.ac.uk/cec/staff/postgradresearch/index.php?web_id=Eric_Flores_de_Gracia&tab=profile [[alternative HTML version deleted]]

...n error message "can't run the model with multiple levels", or something alike.. Is there a way around this or should I start looking fo another way of analysing the data? Thank you for your time. Pekka Kontiainen Ph.D. Candidate Bird Ecology Unit Department of Environmental and Biosciences P.O. Box 65 (00014) University of Helsinki FINLAND Phone +35844-5496858 ----- Edelleenl?hetetty viesti p??ttyy ----- Pekka Kontiainen Ph.D. Candidate Bird Ecology Unit Department of Environmental and Biosciences P.O. Box 65 (00014) University of Helsinki FINLAND Phone +35844-5496858

?typeof? is your friend here: > typeof(`[`) [1] "special" > typeof(mc[[1]]) [1] "symbol" > typeof(mc2[[1]]) [1] "special" so mc[[1]] is a symbol, and thus not a primitive. - Lukas > On 28 Mar 2017, at 14:46, Michael Lawrence <lawrence.michael at gene.com> wrote: > > There is a difference between the symbol and the function (primitive >

I've noticed that I get a warning message every time a package masks some functions from Hmisc. The warning message says : Warning message: In identical(get(., i), get(., lib.pos)) : ignoring non-pairlist attributes This happens with eg: library(plyr) library(xtable) I think I've seen this passing by before, but I'm not sure any more. Just thought I'd mention it. Cheers Joris

Im developing an application with R and Gtk+. It's just a simple GUI which helps new users to interactuate with R. Thing is, when you do a statistical analysis, I also want to provide a HTML report, but HTMLStart doesnt work propperly when executing from TinnR. It does create the file but not empty, I've tried some examples from different websites, and it's always the same.. it works

Can someone help me with performing leave-out-one cross validation using R (model built is a Cox model)? Thanks. --------------------------------------------- David Verbel, MPH Senior Biostatistician Aureon Biosciences 28 Wells Avenue Yonkers, NY 10701 Phone: (914) 377-4021 Fax: (914) 377-4001 --------------------------------------------- [[alternative HTML version deleted]]

Dear gurus, I was utterly surprised to learn that one of my examples illustrating the need of match.fun() doesn't give me the expected result. center <- function(x,FUN) FUN(x) center(1:10, mean) mean <- 4 center(1:10, mean) Used to give me the error message "could not find function FUN". Now it just works, even though I didn't expect it to. I believe this is at least

...ution, I need to do this with multiple data.frames that slot-wise are identical, but each slot length is different between data frames. Seems like there should be an easy solution, but I just have not stumbled across it in the documentation. Thanks, John C. Sweval Database Architect Illumigen Biosciences, Inc. Email: jsweval@illumigen.com Phone: 206-378-0400 Fax: 206-378-0408 [[alternative HTML version deleted]]

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 I have the following data: > y [1] 0.000 0.004 0.008 0.016 0.024 0.032 0.044 0.064 0.072 0.088 0.108 0.140 [13] 0.156 0.180 0.208 0.236 0.264 0.296 0.320 0.360 0.408 0.444 0.472 0.524 [25] 0.576 > x [1] 100 200 300 400 500 600 700 800 900 1000 1100 1200 1300 1400 1500 [16] 1600 1700 1800 1900 2000 2100 2200 2300 2400 2500 I'd

Dear all, Before I file this as a bug, I wanted to check if I didn't miss something. The help page of lag() says that the function returns a time series object. It actually does return something that looks like a ts object (the attribute tsp is set). But when using a vector, the class "ts" is not added to the result: > avec <- 1:10 > lag(avec) [1] 1 2 3 4 5 6 7 8

Hi, My R version is old, and I would like to update it. How can I update R and keep all the libraries? Thanks Best, Yunfei Li -------------------------------------------------------------------------------------- Research Assistant Department of Statistics & School of Molecular Biosciences Biotechnology Life Sciences Building 427 Washington State University Pullman, WA 99164-7520 Phone: 509-339-5096 http://www.wsu.edu/~ye_lab/people.html [[alternative HTML version deleted]]

Hi, I would like to have an index for a column in a matrix encoded in a cell of the same matrix. For example: x = matrix(c(11,12,13,1, 21,22,23,3, 31,32,33,2),byrow=T,ncol=4) In this case, column 4 is the index. I then access the column specified in the index by: > for (i in 1:3) print(x[i,x[i,4]]) [1] 11 [1] 23 [1] 32 > > for (i in 1:3) {x[i,x[i,4]] <- x[i,x[i,4]] + 5} > x

Apparently, as.POSIXlt takes one o'clock as the start of the day : > as.POSIXlt(0,origin="1970-01-01") [1] "1970-01-01 01:00:00 CET" > as.POSIXlt(0,origin="1970-01-01 00:00:00") [1] "1970-01-01 01:00:00 CET" > as.POSIXlt(0,origin="1970-01-01 23:59:59") [1] "1970-01-02 00:59:59 CET" Cheers -- Joris Meys Statistical

I would like to prepare the data for barplot. But I only have the data frame now. x1=rnorm(10,mean=2) x2=rnorm(20,mean=-1) x3=rnorm(15,mean=3) data=data.frame(x1,x2,x3) If there a way to put data within a specific range? The expected result is as follows: range x1 x2 x3 -10-0 2 5 1 (# points in this

Dear all, I run a gamm with following call : result <- try(gamm(values~ s( VM )+s( RH )+s( TT )+s( PP )+RF+weekend+s(day)+s(julday) ,correlation=corCAR1(form=~ day|month ),data=tmp) )" with mgcv version 1.6.2 No stress about the data, the error is not data-related. I get : Error in isS4(x) : object 'VM' not found What so? I did define the dataframe to be used, and the

Dear all, I want to apply a weighted median on a huge dataset, and I remember a function from the package R.basic that could do this using an internal sorting algorithm qsort. This speeded things up quite a bit. Alas, I can't find that package anywhere anymore. There is a weighted.median function in the package limma too, but I didn't use that before. Anybody who knows what happened to

...d but that did not work either. The S Programming book mentions that a library Rdll.lib needs to be built from source, but I am not sure if this is still the recommended route. I am using Microsoft Visual C++ Toolkit in Windows XP Pro. Thank you for any advice. Christos Hatzis, Ph.D. Nuvera Biosciences, Inc. 400 West Cummings Park Suite 5350 Woburn, MA 01801 Tel: 781-938-3830 www.nuverabio.com <http://www.nuverabio.com/> [[alternative HTML version deleted]]

Dear all, I'm updating a package regarding a new type of models, and I'm looking to extend the formula interface with two functions (L() and R() ) for construction of these models. I want to use as much of the formula interface as possible, and hoped to do something similarly to I(). I know the I() function does nothing more than add the class "AsIs". I've been browsing the