search for: exposur

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 colum...

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 with(tolerance,tapply(exposure,list(male),FUN=mean)) #?????? 0??????? 1 #1.246667 1....

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),...

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...

...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) exposure(effort) eform I am lost at the point of finding the equivalent code for 'exposure'. Having looked at a few forums and 'googled'. I thought 'offset', used as offset=(log(Eff)) or the equivalent +offset(log(Eff)) would produce the desired effect. Incidentally...

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 kicked all releves with missing environmental variables. Both, metaMDS and envfit then work without problems. Now I want to bring the releves with missing environm...

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 see if there is a treatment interaction. I've looked at case-crossover designs and...

...tion to work correctly. Consider the following example: ------ library(VGAM) library(MASS) attach(pneumo) pneumo # Inspect the format of the original dataset # Restructure the pneumo dataset into a different format pneumo2 <- data.frame(matrix(ncol=3, nrow=24)) colnames(pneumo2) <- c('exposure.time', 'severity', 'freq') pneumo2[,1] <- rep(pneumo[,1],3) pneumo2[,2] <- as.ordered(c(rep('normal',8),rep('mild',8),rep('severe',8))) pneumo2[1:8,3] <- pneumo[,2] pneumo2[9:16,3] <- pneumo[,3] pneumo2[17:24,3] <- pneumo[,4] pneumo2 #...

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 mean <- 6 sd <- 1 expo <- rlnorm(n, mean, sd) # Then I assign each study subject a p...

...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 ge...

...z32-07 workgroup = LEVEL3 server string = %h # do not change anything in the [global] section beyond this point. # Security settings to allow operation with Windows domain credentials. # Misconfiguration will impact the availability of this system # and is a severity 3 exposure. security = domain password server = * allow trusted domains = yes encrypt passwords = yes client use spnego = yes # We don't want Samba to become a master browser on the network, and # never act as the Primary Domain Controller. # Misconfi...

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, Kim -- f?r nur 19,99 Euro/mtl.!* http://port...

...the result of the aov functions are calculating everything as a within subject. the table looks like this: SerialNo breed treatment distance_1 distance_2 1 c57 dfp 235 3253 etc. I changed it to look like this: SerialNo breed treatment exposure distance 1 c57 dfp 1 235 1 c57 dfp 2 3253 etc. Then I do: dt<-aov(distance~(exposure*treatment*breed)+Error(SerialNo/exposure) + (treatment*breed), dataframe) what am I doing wrong? thank you. [[altern...

I am attempting to estimate LC50 (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...

...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 = hetcor(dat[,c('intent','exposure','benefit','norms','childBarrier','parentBarrier','knowBenefit','recuse')],use="pairwise.complete.obs")$correlations modMat = matrix(c( 'exposure -> intent', 'gam11',NA, 'benefit -> intent', 'gam12&...

-----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, 11571.2744, 15400.9318, 7786.4093, 1173.0635 , 2630.1493 , 4408.2271 ,6333.0374 ,8070.0090, 751.3506) RTotEmp<-c(...

-----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