similar to: Tables extraction in R ?

> Dear All, > > I have built a survival cox-model, which includes a covariate * time interaction. (non-proportionality detected) > I am now wondering how could I most easily get survival predictions from my model. > > My model was specified: > coxph(formula = Surv(event_time_mod, event_indicator_mod) ~ Sex + > ageC + HHcat_alt + Main_Branch + Acute_seizure +

Hi All, The data set that I have is a cluster data, and I want to run a HLM mixed model with multi-level response. Here is my data set: response: - Level (num: 1, 2, 3, 4, 5 - 5 levels) Covariates: - Type (Factor: A, B, C - 3 levels) - yr (num: 2006, 2007, ...) - Male (num: 0=not Male, 1=Male - 2 levels) - Ethnicity (Factor: A, B, H, ..., - 7 levels) - ELL (num: 0, 1, - 2

Hi there, Could somebody help me disect this legacy R script I inherited at work, I have two questions: 1. I've tried to upgrade our R version from 1.6.2 (yeah, I know), to R 2.0, but some of the lines in this script are not compatible with R 2.0, could someone help me figure out where the problem is? 2. the jpeg generated (attached) seems to be off on some of the data, is there a better way

Dear R-helpers, I don't see a difference between the following two plots of effect objects, which I understand should be different. What am I missing? require(doBy) require(effects) data(budworm) m1 <- glm(ndead/20 ~ sex + log(dose), data=budworm, weight=ntotal, family=binomial) m1.eff <- all.effects(m1) plot(m1.eff, rescale.axis = FALSE, selection = 2, main = 'rescale =

Hello, This is Elaine. I am trying a path analysis using lavaan Package. There are three explanatory variables: X, Z, and M. The response variable is Y. A, b, and c have direct effects on Y. On the other hand, X and Z also have direct effects on M. In other words, X and Z have indirect effects on Y. I found the code example of lavaan package describes only one indirect effect as below. Please

Hello, I am trying to see whether there has been a significant difference in whether people experienced damages from wildlife in two different years. I therefore have two columns: year 1: yes no no no yes yes no year 2: no yes no yes I wanted to do a chisq.test, but if I enter it this way: chisq.test(year1, year2) I get the error saying the columns are two different lengths. So then I tried

I have a set of data that includes various data columns. One if the survival time and another if a continuous variable of ages. I want to put the ages into intervals so that I can then perform the Kalpan Meier test. I am trying to use the following code to build a column with the age group numbers in agecatagory<-c( ) for (i in 1:137) { { if(age[i]<=46) {agecat[i]<-1} if(age[i]>46

Hi, I have a large amount of data that I would like to create a histogram of and plot and do things with in R. It is pretty much impossible to read the data into R, so I have written a program to bin the data and now have a list of counts in each bin. Is it possible to somehow import this into R and use hist(), so I can, for instance, plot the probability density? I have looked at the help page

Dear R users, The package proptest (http://www.davidkraus.net/proptest/) has been recently released to CRAN. The package provides functions for testing the proportional hazards assumption in the Cox model for right censored survival data. Two types of tests for identifying nonproportional covariates are implemented: * data-driven Neyman type smooth tests (using both nested and all-subsets

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Release to fix regression in nouveau in 2.4.60, but brings lots of android changes and other fixes. Ben Skeggs (2): nouveau: restore check that avoids multiple user bos per kernel bo configure.ac: bump version to 2.4.61 for release Chih-Wei Huang (3): android: simplify the including rule of subdirs android: get rid of

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 libdrm 2.4.66 release, mainly for new nouveau API. lots of other changes in here as well though. Ben Skeggs (14): nouveau: import and install a selection of nvif headers from the kernel nouveau: move more abi16-specific logic into abi16.c nouveau: move object functions up, to avoid future foward decls nouveau: make it

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Adam Cheney (1): Add new DRM_MODE_CONNECTOR and _ENCODER defines Chris Wilson (1): intel: Avoid overcounting fences when emitting self-referential relocs Damien Lespiau (3): intel/skl: Add SKL PCI ids intel/skl: Add gen9 to the buffer manager init intel/skl: add gen9 to the CS decoding init Dave Airlie (1):

Chih-Wei Huang (1): android: add virgl to be a valid driver Daniel Vetter (7): headers: Remove _DRM_GEM tests/kms-steal-crtc: Use correct includes headers: Update drm_i915.h headers: Update drm_sarea.h headers: Update drm_fourcc.h headers: Update drm_mode.h headers: Update drm.h Emil Velikov (1): proptest: print signed values correctly Eric

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Bernd Kuhls (1): tests/etnaviv: link against libdrm Christian Gmeiner (1): etnaviv: remove struct etna_specs Emil Velikov (2): amdgpu: add amdgpu_bo_va_op_raw to the symbol check configure.ac: bring back pthread-stubs check Erik Faye-Lund (1): tegra: update symbol-check Junwei Zhang (1): amdgpu: add REPLACE and

Alex Deucher (2): amdgpu: add new marketing names amdgpu: add new marketing names Dan Willemsen (1): Convert to Android.bp Daniel Stone (1): libs: Tie DSO minor versions to libdrm version Elliott Hughes (1): readdir_r is deprecated. Enrico Weigelt, metux IT consult (1): Fix FTBS on undefined clock_gettime() and asprintf() Inseob Kim (1): Export

hello I like to comput median and mean point and intreval estimation by polya posterior when super population is normal(35,sd=5) with 1000 repeat. 500 and 25 are population size and sample size. how do i do? help me please [[alternative HTML version deleted]]

Hi, I am looking for some code for Moran's I. Has anyone previously done this? I have been unable to find it in the search engines. James