search for: exposures

Hello, I have a loop with 1000 repetitions which includes OR computation of an exposure factor and outcome. I compute OR like this: t<-table(exposure,outcome) oddsratio(t)$measure["estimate"] This gives me the estimates for exposure=0 and exposure=1 but exposure=0 is the reference group and i need only the estimate for exposure=1. I specified a matrix OR with 3 columns (for OR

Hi, Try either: tolerance <- read.csv("http://www.ats.ucla.edu/stat/r/examples/alda/data/tolerance1.txt") ?aggregate(exposure~male,data=tolerance,mean) ?# male exposure #1??? 0 1.246667 #2??? 1 1.120000 #or ?library(plyr) ?ddply(tolerance,.(male),summarize,exposure=mean(exposure)) #? male exposure #1??? 0 1.246667 #2??? 1 1.120000 #or

Dear R helpers exposure <- data.frame(id = c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20), ead = c(9483.686,50000,6843.4968,10509.37125,21297.8905,50000,706152.8354, 62670.5625, 687.801995,50641.4875,59227.125,43818.5778,52887.72534,601788.7937, 56813.14859,4012356.056,1419501.179,210853.4743,749961,6599.0862), pd =

Dear List, ? I want to fit a CPH survival model with a time-dependent cumulative exposure variable. More specifically, I have exposure measurements for all cohort members at discrete times (yearly measurements) during follow-up. ? I suppose that when estimating the model parameters the most correct would be to have the cumulative exposure for everybody at risk at each death time, as I cannot

R-helpers, I am hoping to find someone who uses both R and program Stata for GLMs. I am a beginner R user, finding my own way through; learning code etc. at the same time as learning the statistics I need to complete my project. What I have is the code from Stata and am trying to reproduce the same analysis in R - my program of choice. . glm count md ms rf sg, family(poisson)

Dear r-Community, Step1: I would like to calculate a NMDS (package vegan, function metaMDS) with species data. Step2: Then I want to plot environmental variables over it, using function envfit. The Problem: One of these environmental variables is cos(EXPOSURE). But for flat releves there is no exposure. The value is missing and I can't call it 0 as 0 stands for east and west. Therefore I

Hi, I'm hoping that someone will be able to help. I would like to compare how covariates associate with the risk of a binary outcome during two periods. Period 1 will be non-exposure to a treatment and period 2 will be exposure to a treatment. The same individuals will be examined in each group but I want to be able to compare the association of certain covariates between the two groups to

R-Experts: I am using the vglm function of the VGAM library to perform proportional odds ordinal logistic regression. The issue that I would like help with concerns the format in which the response variable must be provided for this function to work correctly. Consider the following example: ------ library(VGAM) library(MASS) attach(pneumo) pneumo # Inspect the format of the original dataset

Hi R help list I'm simulating logistic regression data with a specified odds ratio (beta) and have a problem/unexpected behaviour that occurs. The datasets includes a lognormal exposure and diseased and healthy subjects. Here is my loop: ors <- vector() for(i in 1:200){ # First, I create a vector with a lognormally distributed exposure: n <- 10000 # number of study subjects

...at gives an overview of one of my current issue. Although my working setting implies a different parametric model (which cannot be framed in the glm), I guess that what I'll get from the following example it would help for the next steps. Anyway here it is. Firstly I simulated from a series of exposures, a series of deaths (given a model, Gompertz, and a probability distribution, Poisson). Then a multiply both deaths and exposures by a constant. Finally I fitted with glm the simulated data sets and I compared the deviance residuals. They're different by a certain constant, but I could not get...

Hello I'm running Samba version 3.0.2a on Solaris 9 and can not get access to my defined shares. My config is below. I believe my configuration is good and I can see the server in the Windows Network browser but can not authenticate. Any one have any ideas? Thanks Steve [global] netbios name = f2z32-07 workgroup = LEVEL3 server string = %h # do not change

Dear all, I would like to calculate a cca (package vegan) with species and environmental data. One of these environmental variables is cos(EXPOSURE). The problem: for flat releves there is no exposure. The value is missing and I can't call it 0 as 0 stands for east and west. The cca does not run with missing values. What can I do to make vegan cca ignoring these missing values? Thanks a lot,

Hi, I just started with R and I found that there are many options to rearrange the data to do mixed models. I want to use the reshape function. I have 2 between subject variables and one within. I was able to change the data structure but still - the result of the aov functions are calculating everything as a within subject. the table looks like this: SerialNo breed treatment distance_1

...50 (analogous to LD50, but uses exposure concentration rather than dose) by fitting a Weibull model; but I can't seem to get it to work. From what I can gather, I should be using survreg() from the survival package. The survreg() function relies on time-to-event data; my data result from 96 h exposures (i.e., dead or alive after a fixed period; 96 h). I've tried the following (doesn't work): > conc <- c(10.3, 10.8, 11.6, 13.2, 15.8, 20.1) # Exposure concentrations > orign <- c(76, 79, 77, 76, 78, 77) # Original number of subjects > ndead <- c(16, 22, 40, 69, 78, 77) #...

Hello list, I'm trying to make sure that I'm performing a path analysis correctly using the sem package. the figure at http://flame.cs.dal.ca/~sstewart/regressDiag.png has a detailing of the model. The challenge I'm having is that reuse is an indicator (0/1) variable. Here's the code I'm using: corr =

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 ========================================================== == == Subject: Exposed clear text of domain machine == account password in debug logs (log == level >= 5) == CVE ID#: CAN_2006-1059 == == Versions: Samba Samba 3.0.21 - 3.0.21c (inclusive) == == Summary: The winbindd daemon writes the clear text ==

Greetings all, please consider the following data: #Build Data frame Slope<-c(1.291370, 12.208500, 2.110930, 0.578990, 5.019520, 0.807444, 0.554079 , 1.257080, 0.241504 , 0.184337 , 0.383044 , 0.342021) Exposure<-c(790.54, 1167.79 , 845.58 , 1082.47 , 1189.61 , 677.17 , 2058.56 , 469.09 , 112.02 , 803.31 , 254.14 ,1336.16) FwyDist<-c(11809.4222 ,10623.0458, 12279.6271,

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 ========================================================== == == Subject: Exposed clear text of domain machine == account password in debug logs (log == level >= 5) == CVE ID#: CAN_2006-1059 == == Versions: Samba Samba 3.0.21 - 3.0.21c (inclusive) == == Summary: The winbindd daemon writes the clear text ==

I just confirmed the following bug on my firefox. http://secunia.com/advisories/14820/ Mozilla/5.0 (X11; U; FreeBSD i386; en-US; rv:1.7.5) Gecko/20050219 Firefox/1.0 (I think my firefox is a month or two behind, from ports, but the advisary indicates both 1.0.1 and 1.0.2 are effected.) FreeBSD localhost 5.3-RELEASE FreeBSD 5.3-RELEASE #0: Fri Nov 5 04:19:18 UTC 2004

We specifically do not update the disk i_size if there are ordered extents outstanding for any area between the current disk_i_size and our ordered extent so that we do not expose stale data. The problem is the check we have only checks if the ordered extent starts at or after the current disk_i_size, which doesn''t take into account an ordered extent that starts before the current