search for: davidbickel

...of the value returned by glm. (I'm particularly interested in logistic regression.) If not, what is a good way to use R to compute Hessians or other partial derivatives of log likelihoods? I would appreciate any guidance. David _______________________________________ David R. Bickel http://davidbickel.com Research Scientist Pioneer Hi-Bred International (DuPont) Bioinformatics 7200 NW 62nd Ave.; PO Box 184 Johnston, IA 50131-0184 515-334-4739 Tel 515-334-4473 Fax david.bickel at pioneer.com This communication is for use by the intended recipient and ...{{dropped}}

Does anyone know how to use the equivalent of debug() on an S4 method? I would like R to enter the browser not for the generic function, but for the method of the class that I specify. Thanks, David _____________________________ David Bickel http://davidbickel.com Research Scientist Pioneer Hi-Bred International Bioinformatics & Exploratory Research 7250 NW 62nd Ave., PO Box 552 Johnston, Iowa 50131-0552 515-334-4739 Tel 515-334-6634 Fax david.bickel at pioneer.com, bickel at prueba.info This communication is for use by the intended recipient and...

Given a numeric vector of observations, does R have any generic way to estimate the parameters of commonly used distributions (exponential, gamma, etc.) without numerically optimizing the likelihood function? Thanks, David _______________________________________ David R. Bickel http://davidbickel.com Research Scientist Pioneer Hi-Bred International Bioinformatics & Exploratory Research 7250 NW 62nd Ave., PO Box 552 Johnston, Iowa 50131-0552 515-334-4739 Tel 515-334-6634 Fax david.bickel at pioneer.com, bickel at prueba.info This communication is for use by the intended recipient and...

...ountered a problem and needs to close. We are sorry for the inconvenience." I downloaded 1.9.0 and installed it, but when I launch it, I get this message: "Fatal error: Invalid HOMEDRIVE." I would appreciate any assistance. -Dave ______________________________ David Bickel http://davidbickel.com This communication is for use by the intended recipient and ...{{dropped}}

...n of the error, I would like browser() to be called at the level of numeric.sapply(), so that I can examine x. I'm wondering if this can be done by modifying the default error handling. Using try() with browser() didn't work. Thanks, David _____________________________ David Bickel http://davidbickel.com Research Scientist Pioneer Hi-Bred International Bioinformatics & Exploratory Research 7250 NW 62nd Ave., PO Box 552 Johnston, Iowa 50131-0552 515-334-4739 Tel 515-334-6634 Fax david.bickel at pioneer.com, bickel at prueba.info This communication is for use by the intended recipient and...

...and for any positive integer n? I need to perform the integration within a function with q and n as arguments. The function integrate() is giving me this message: "evaluation of function gave a result of wrong length" David _______________________________________ David R. Bickel http://davidbickel.com Research Scientist Pioneer Hi-Bred International (DuPont) Bioinformatics and Exploratory Research 7200 NW 62nd Ave.; PO Box 184 Johnston, IA 50131-0184 515-334-4739 Tel 515-334-4473 Fax david.bickel at pioneer.com, bickel at prueba.info This communication is for use by the intended recipient a...

What are some effective ways to leverage the strengths of R and Mathematica for the analysis of a single data set? More specifically, are there any functions that can assist with any of the following? 1. Calling an R function from Mathematica. 2. Calling a Mathematica function from R. 3. Using XML or another reliable data format to pass vectors, matrices, and/or lists from one environment to

...ensive command while I execute other R commands. I want to do this in the same R session, rather than with multiple calls to R, so that all my modified objects will stay in the same workspace. I would appreciate any help with this. Thanks, David ______________________________ David Bickel http://davidbickel.com Research Scientist Pioneer Hi-Bred International Bioinformatics & Discovery Research 7250 NW 62nd Ave., PO Box 552 Johnston, Iowa 50131-0552 706-736-9151 Home 515-334-4739 Work 515-334-6634 Fax david.bickel at pioneer.com, bickel at prueba.info This communication is for use by the intend...

...this: > plot.postscript.file(file = 'graph.ps') If no such function exists, I'd be interested in a way to use existing R functions to do this under UNIX or Windows, preferably without a system call to GhostView (gv). Thanks, David _____________________________ David Bickel http://davidbickel.com Research Scientist Pioneer Hi-Bred International Bioinformatics & Exploratory Research 7250 NW 62nd Ave., PO Box 552 Johnston, Iowa 50131-0552 515-334-4739 Work 515-334-6634 Fax david.bickel at pioneer.com, bickel at prueba.info This communication is for use by the intended recipient and...

...cs (or sample quantiles or extreme values) for a given distribution such as a normal distribution? Both exact and asymptotic calculations are of interest. I am already aware of the 'quantile' function of 'stats'. David _______________________________________ David R. Bickel http://davidbickel.com Research Scientist Pioneer Hi-Bred International (DuPont) Bioinformatics and Exploratory Research 7200 NW 62nd Ave.; PO Box 184 Johnston, IA 50131-0184 515-334-4739 Tel 515-334-4473 Fax david.bickel at pioneer.com, bickel at prueba.info This communication is for use by the intended recipient a...

Is there a function that can shade the region between two arbitrary curves? For example, how can I fill in the area between these two plotted curves? > x <- 1:10 > y <- x + rnorm(10) > z <- x + 10 + rnorm(10) > plot(x, y, type = "l", ylim = c(0,20)) > lines(x, z) I would appreciate any help. David

....list, function(x){x[1]}) vec2 <- sapply(f.list, function(x){x[2]}) The values of vec1 should be as iid as possible and the values of vec2 should be as iid as possible. Any help with what I need to add to f would be greatly appreciated. David _____________________________ David Bickel http://davidbickel.com Research Scientist Pioneer Hi-Bred International Bioinformatics & Exploratory Research 7250 NW 62nd Ave., PO Box 552 Johnston, Iowa 50131-0552 515-334-4739 Tel 515-334-6634 Fax david.bickel at pioneer.com, bickel at prueba.info This communication is for use by the intended recipient and...

...ternative may be to have R and an R package behave as a stand-alone application that can be called from the UNIX command line. Is there any documentation on this? S-PLUS can do something like this. I will be using Red Hat Linux. Thank you, David ______________________________ David Bickel http://davidbickel.com Research Scientist Pioneer Hi-Bred International Bioinformatics & Discovery Research 7250 NW 62nd Ave., PO Box 552 Johnston, Iowa 50131-0552 706-736-9151 Home 515-334-4739 Work 515-334-6634 Fax david.bickel at pioneer.com, bickel at prueba.info This communication is for use by the intend...

Does anyone know how to interrupt R in RedHat? Neither control-c nor Esc is working. What I don't want to do is close the window or kill the entire R process. Thanks, David This communication is for use by the intended recipient and ...{{dropped}}

Neither biocLite nor the GUI menus can install packages on my system. Here is relevant output: > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2 minor 11.1 year 2010 month 05 day 31 svn rev 52157 language R version.string R version 2.11.1 (2010-05-31) > source("http://bioconductor.org/biocLite.R") BioC_mirror =

What are some reliable R functions that can compute the value of a convergent series? David -- David R. Bickel, PhD Associate Professor Ottawa Institute of Systems Biology Biochem., Micro. and I. Department Mathematics and Statistics Department University of Ottawa 451 Smyth Road Ottawa, Ontario K1H 8M5 http://www.statomics.com Office Tel: (613) 562-5800 ext. 8670 Office Fax: (613) 562-5185