similar to: Own classes in "histogram"

Hi everyone I have the following problem: My dataframe has 3 variables: ID, Year and and an outcome variable. The dataframe contains repeated measurements because the subjects filled out a questionnaire every year. The time span covers 2 years. Now I want to check if there is a significant change in the outcome over the 2 years with a paired wilcox.test. The problem: Not every subject

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

Using R-2.12.1 and latticeExtra-0.6-14, I would like to understand why a rootogram displaying samples from the Poisson distribution looks like I expected it, whereas a rootogram using the normal distribution does not: library(latticeExtra) rootogram(~rpois(1000, lambda = 50), dfun = function(x) dpois(x, lambda = 50)) rootogram(~rnorm(1000), dfun = function(x) dnorm(x,mean(x),sd(x))) I

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

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

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

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

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

I have created an xyplot of a time series with the following code... win.graph(width = 10, height = 7) panel1 = function(x, y) { panel.loess(x, y, lwd=2.5, span=0.5, col="gray") panel.xyplot(x, y, pch=19, col="blue", cex=1.25) } xyplot(oneplusdensity ~ year, data=figdata, aspect="fill", cex=1.5, xlab=NULL, ylab=expression("Crabs per

Dear all I want to run a simple multiple linear regression model, say: mod <- lm(y~x1+x2+x3+x4) I have programmed a custom algorithm to calculate the p-values, because there is clustering present. I ordered the variables according to their p-value and want to stepwise eliminate the variables with the highest p-values. My question: How can I use "update" when I have only the

I need to do a somewhat complex rearrangement of an existing dataframe that involves an "on-the-fly" conditional summary as well. I've tried to make the example as simple as possible. Code to produce the object being acted on is at the end of the message. I hope this isn't too much for a single posting. Any guidance will be much appreciated. The input example looks like

Hi there, I hope this question is not as stupid as the one before ... I tried to shorten my histogram (because the distribution is quite skewed and I simply don't want to see the long tail but still use the histogram plot). How can I do something like this? (The example does not work but I don't know why...) data <- rnorm(100) # as example, of course this is not skewed... h <-

Hi to all I use a loop to plot 9 different histograms of 9 different transformations of one dataset (x). I want to label the histograms with the mathematical expression of each transformation (e.g. x^3). For that I've prepared a vector with the labeling names for "expression". > trans.expr <- c("x^3", "x^2", "x", "frac(1,x)",

Hello all, I have been looking through the internet and through the wiki, trying to find something which explains in details the available interfaces with which one might be able to insert some code to make some slight additions to Xen functionality. I am doing some research on the possibility of adding some extra functionality to Xen, and my supervisor has mentioned that these are things worth

Hi I am trying to get a table of means of parameter 1 across BINS of parameter 2. I am working in proteomics and a sample of my data is as follows cluster-age clock-rate(evolutionary rate) scopclass 0.002 10 A 0.045 0.1 B 0.13 15 A 0.15 34 D .... .... ....

TestHistogram_A is a function that plots a histogram, then returns 3. For me it fails to display the plot, at least in the order given below. TestHistogram_B is a function that plots a histogram, then exits. The only difference between the 2 functions is that B does not return 3. It displays the plot. This problem does not occur when histogram is replaced by boxplot Their use

Dear R users and helpers, I'm trying to find an example of a histogram plot as an inset (upper right or left corner) of another histogram. Anyone has an example of that? Thanks for your help, V?ctor. -- ----------------------------------------------------------------------- V?ctor Homar Santaner Grup de Meteorologia Edif. Mateu Orfila Tel: +34 971 17 1376