similar to: Remove missing observations

Hi everyone I am stuck on specifying my own maximum likelihood function for a special poisson model. My poisson model is as follow: O ~ Pois(b*N + b*RR*E) With O = observed cases b = constant (known) N = number of unexposed persons (known) E = number exposed persons (known) RR = relative risk (value is assumed under a scenario, e.g. RR=2.0) I used rpois to simulate the values of O for several

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

Hello, I just wanted to save a graphic in the pdf-format. But id failed: Fehler: Invalid font type Zus?tzlich: Warning messages: 1: font family not found in PostScript font database 2: font family not found in PostScript font database I use R 2.6.0 with all packages updated recently. Any idea? Thank you very much! Best, Denis __ Denis Aydin mailto:Denis.Aydin at stud.unibas.ch Montag, 22.

Hi, I try to make a histogram from a variable that contains the number of shoots from about 1000 individuals from a specific plant species (the range is 1-110). Those numbers are highly skewed to the right. My question is: how can I make my own classes with the lattice "histogram"? I tried it with "breaks=c(0,5,10,15,20,25,110)" but my "25-110"-class is presented

I would like to build a package for the HT HG-U133+ PM arrays from affy, but I can't find any good documentation on how to go about it. Naively using makecdfenv's make.cdf.package() causes R to seg-fault. I'm unfamiliar with the CDF format as such, but I'm guessing that it's changed somewhat because the PM arrays no longer have P/A and mismatches. I'm looking to build

Dear R-helpers I have a question related to logistic GAM models. Consider the following example: # Load package library(mgcv) # Simulation of dataset n <- 1000 set.seed(0) age <- rnorm(n, 50, 10) blood.pressure <- rnorm(n, 120, 15) cholesterol <- rnorm(n, 200, 25) sex <- factor(sample(c('female','male'), n,TRUE)) L <-

Dear all, I am trying to migrate one of our systems to dovecot. We were using uw-imapd in black-box-mode before. Since the home-directory of a user is not unique per definition on our systems, we have kind of a different setup. As an example: I got a user, let's say user1. It's home directory is /local/web/websites/5850844095/www. This directory has the following permissions: # ls -la

Hi everyone I have the following situation: In a longitudinal study, subjects fill out a questionnaire every year (repeated measurements over time). Also, the subjects are nested within departments. There is an intervention going on over time. The outcome variable is continuous. Now I'd like to analyse two things: 1. Is there a significant change over time? I think this is done by a

Hi I have a question regarding an output of a binomial lmer-model. The model is as follows: lmer(y~diet * day * female + (day|female),family=binomial) The corresponding output is: Generalized linear mixed model fit by the Laplace approximation Formula: y ~ diet * day * female + (day | female) AIC BIC logLik deviance 1084 1136 -531.1 1062 Random effects: Groups Name Variance

Hi everybody, I have a strange behavior of the as.Date() function. For example: as.Date(36525, origin="1900-01-01'") I would expect to get 2000-01-01. But R gives me [1] "2070-01-01" I tried with several origins but the results are still the same. It seems that the origin does not apply and remains 1970-01-01. Does anyone can give me a hint how to

Dear list, I don't know if this is by design or possibly a bug, but there is a difference when removing e-mail when a user is over filesystem quota when indexing is on or off. All is working fine on our servers when indexing is on, users can remove then purge e-mail when over quota, but when we turn off indexing (by setting INDEX=MEMORY) it doesn't work. Users can remove, but

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

Hi, I have two datasets stored in tab-separated format in the following way file1: country year1 year2 Germany var1 var1 Hungary var1 var1 file2: country year1 year2 Germany var2 var2 Hungary var2 var2 I can easily read in these files, but how can I merge them as a panel dataset? Thanks, Viktor

I ran the following code: And I run into problems with the last line of code (when it says hn<-......). I keep getting an error code: Error in distsamp(~hab ~ 1, peldist, keyfun = "halfnorm", output = "density", : object 'A' not found I would appreciate any and all help. rm(list=ls(all=TRUE)) #clear the computer's

Hi List, I'm trying to install aroma.affymetrix, which depend on many packages, affyio is one of them. The task is done on 32-bit Linux without problem. The same task was done on Solaris 9 with R-2.5.1 without problem. Following is the error message when installation affyio. I did try to put a file stdint.h I got from the web in the path, but got more error messages. Jun root@biostatsun1

Full_Name: Hamit AYDIN Version: R v2.1.1 OS: winxp Submission from: (NULL) (194.27.40.230) When I try to run sciviews the consol opens but when I try to import text, spss... it tries to load rmcdr pakage gives me "the pakage can not be find inthe Package or bundle 'Rcmdr' was not found in C:\PROGRA~1\R\rw2011\library. Would you like to install it now?" I know the pakage in

Dear R Experts, I am trying to produce a legend for a series of plots which are generated in a loop. The legend is supposed to look like this: 2000: gamma=1.8 where gamma is replaced by the greek letter and both the year and the value of gamma are stored in variables. Everything works fine as long as I have only one data series: year = 2001 g = 1.9 plot(1) legend('top',

I would like to know how I can assign a name to a name. I have a dataset that has different years in it. I am writing scripts using R and I would like to give a month a generic name and then use the generic name to do different analysis. The reason for the generic name would be so that I only have to change one thing if I wanted to change the year. For example. Year1 = 1999 datayear <-

I am graphing longitudinal data from three time points. I'd like to draw a solid line from point 1 to point 2, and then a dashed line from point 2 to point 3. It works if I do it in two steps: > first.vector <- c(mean(year1$variable1), mean(year2$variable1)) > second.vector <- c(NA, mean(year2$variable1), mean(year3$variable1)) > plot(first.vector, type="b",

I need to apply a yearly inflation factor to some wages and supply some simple sums by work category. I have gone at it with a brute force "for" loop approach which seems okay as it is a small dataset. It looks a bit inelegant and given all the warnings in the Intro to R, etc, about using loops I wondered if anyone could suggest something a bit simpler or more efficent? Example: