similar to: How to extract se(coef) from cph?

Hi, How does R automatically load functions so that they are available from the workspace? Is it anything like Matlab - you just specify a directory path and it finds it? The reason I ask is because I found a really nice script that I would like to use on a regular basis, and it would be nice not to have to 'copy and paste' it into R on every startup:

Dear List Sorry if this question seems very basic. Is there a function to pro grammatically find number of processors in my system _ I want to pass this as a parameter to snow in some serial code to parallel code functions Regards Ajay Websites- http://decisionstats.com http://dudeofdata.com Linkedin- www.linkedin.com/in/ajayohri

Dear R gurus, I regularly come across a situation where I would like to apply a function to a subset of data in a dataframe, but I have not found an R function to facilitate exactly what I need. More specifically, I'd like my function to have a context of where the data it's analyzing came from. Here is an example: ### BEGIN ### func<-function(x){ m<-median(x$x) if(m > 2 &

guRus I have say a dataframe, d and I wish to do the following: 1) For each row, I want to take one particular value of the row and multiply it by 2. How do I do it. Say the data frame is as below: OPEN HIGH LOW CLOSE 1931.2 1931.2 1931.2 1931.2 0 0 0 999.05 0 0 0 1052.5 0 0 0 987.8 0 0 0 925.6 0 0 0 866 0 0 0 1400.2 0 0 0 754.5 0 0 0 702.6 0 0 0 653.25 0 0 0 348 0 0 0 801 866.55 866.55

Hi everyone, I am trying to make some publication-quality plots for use in Microsoft Word, but I am having trouble creating high-quality plots that are supported by Microsoft Word. If I use the R plot function to create the figure, the lines are jagged, and the picture is not of high quality (same with JPEG(), TIFF(), and PNG() functions). I have tried using the Cairo package, but it distorts

Hello- A basic question which has nonetheless floored me entirely. I have a dataset which looks like this: Type ID Date Value A 1 16/09/2020 8 A 1 23/09/2010 9 B 3 18/8/2010 7 B 1 13/5/2010 6 There are two Types, which correspond to different individuals in different conditions, and loads of ID labels (1:50)

Hello, I would like, if it is possible, to compare the effect of a variable across regression models. I have looked around but I haven't found anything. Maybe someone could help? Here is the problem: I am studying the effect of a variable (age) on an outcome (local recurrence: lr). I have built 3 models: - model 1: lr ~ age y = \beta_(a1).age - model 2: lr ~ age + presentation

Hi all! How is it possible to estimate standard errors for coef obtained from lme? Is there sth like se.coef() for lmer or what is the anaytical solution? Thank you!

Dear all, i have re-installed Debian on my computer and for some reason I cannot install some packages anymore using the install.package command from the R console. Some packages I can, and some I cannot (Design for instance) this is what I get: > install.packages('Design') Installing package(s) into ‘/usr/local/lib/R/site-library’ (as ‘lib’ is unspecified) --- Please select a CRAN

Dear Help List, I am trying to install a package (DiagnosisMed) on R (R version 2.10.1 (2009-12-14)) on Debian Lenny (architecture i386). I get the following error message: ~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~ * installing *source* package ‘DiagnosisMed’ ... ** R ** data ** preparing package for lazy loading Error in firstlib(which.lib.loc, package) : Tcl/Tk support is not available on this

coef() and vcov() have different dimensions if a model contains alised parameters as the following example illustrates. A search on "alised" gave noting as far as I could see. Is this a known bug? Bendix C ---------------------- Bendix Carstensen Senior Statistician Steno Diabetes Center Niels Steensens Vej 2 DK-2820 Gentofte Denmark tel: +45 44 43 87 38 mob: +45 30 75 87 38 fax: +45

Hi, Does anyone know how to extract fixed effects SE values from generalized linear mixed models estimated using the lmer function in the lme4 library? I searched attributes and structure with no luck. Thanks Frank A. La Sorte, Ph.D. Department of Fisheries and Wildlife Sciences University of Missouri Columbia, MO 65211 USA

Hello,   does anyone have an example on how to use restricted cubic splines function rcs within survfit.cph, if cph (Cox Proportional Hazard Regression) was done with restricted cubic splines (which I made to work)? Thank you. > [[alternative HTML version deleted]]

Is the following an inconsistency, programming glitch or a feature? One would expect that vcov(obj) was the variance-covariance of coef(obj), but apparently this is not the case for polr objects: > x <- rnorm( 100 ) > y <- rnorm( 100 ) > ff <- factor( sample( 1:4, 100, replace=T ) ) > pm <- polr( ff ~ x + y ) > coef( pm ) x y 0.21219010 0.03558506

Hello, I would like to get the likelihood ratio and score tests for specific variables in a multivariate coxph model. The default is Wald, so the tests for each separate variable is based on Wald's test. I have the other tests for the full model but I don't know how to get them for each variable. Any idea? David Biau. [[alternative HTML version deleted]]

Dear list, It looks like set.seed does not work for me (see code below). I am running Debian squeeze. Is this an issue with the distribution or with my understanding of the function? I'm pretty sure, but not certain, it used to work that way on windows > runif(5) [1] 0.6516557 0.5677378 0.1135090 0.5959253 0.3580500 > runif(5) [1] 0.42880942 0.05190332 0.26417767 0.39879073

Dear all R users, Can anyone tell me how I can get estimate of SE of coefficients from, lm() function? I tried following : x = 1:10 lm(x[-1]~x[-10]-1)$coefficients Here I got the est. of coefficient, however I also want to get some "automated" way to get estimate of SE. Regards, [[alternative HTML version deleted]]

Hi All! I am desperately needing some help figuring out how to calculate LD50 with a GLMM (probit link) or, more importantly, the standard error of the LD50. I conducted a cold temperature experiment and am trying to assess after how long 50% of the insects had died (I had 3 different instars (non significant fixed effect) and several different blocks (I did 4 replicates at a time)=

Dear all: I have used the cph function in the rms package. log10 was used to transform the variables, as follows: fit<-cph(pfsurv~log10(x1)+log10(x2),x=T,y=T,surv=T) after I run the nomogram function. I found "variable limits and transformations are such that an infinite axis value has resulted." How to add variable limits in the nomogram function? Thanks a lot *Yao Zhu*

Dear R help list, I am modeling some survival data with coxph and survreg (dist='weibull') using package survival. I have 2 problems: 1) I do not understand how to interpret the regression coefficients in the survreg output and it is not clear, for me, from ?survreg.objects how to. Here is an example of the codes that points out my problem: - data is stc1 - the factor is dichotomous