similar to: save as pdf

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

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

Dear All, I wonder whether anyone has an experience with FAME package written by Jeff Hallman. All my attempts to send him the following problem report did not succeed (the mail system says that my e-mail could not be delivered), so I turn for help to this list. I tried to use your FAME package written for R, but somehow I cannot get it working. I am using Windows XP and the newest R

Mich besch?ftig folgende Fragestellung. Ich kenne die Verteilung (lognormal) zus?tzlich weiss ich das 99%, das 90% und das 1% Quantil. Gibt es in R eine M?glichkeit die Lognormalverteilung zu finden, das heisst den korrespondierenden logmean und logsd? Vielen Dank f?r ihre Hilfe Gruss Yvonne

Dear List, I have problems with the function source() using a url of the kind: http://localhost:5984/path/fn.R I receive Fehler in file(file, "r", encoding = encoding) : kann Verbindung nicht ?ffnen Zus?tzlich: Warnmeldung: In file(file, "r", encoding = encoding) : ?ffnen fehlgeschlagen: HTTP Status war '502 cannotconnect' The url itself is ok since I can: -

Hi, after updating to foreign version 0.8-32, I experienced the following error when I tried to load a SPSS file: Fehler in inherits(x, "factor") : objekt "cp" nicht gefunden Zus?tzlich: Warning message: In read.spss("***l.sav", use.value.labels = TRUE, to.data.frame = TRUE) : ***.sav: File-indicated character representation code (1252) looks like a Windows

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

Dear R-help, We are using the Platform AMD 64 Bit Opteron. Our operating system is SUSE LINUX ENTERPRISESERVER SLES9.0 We are using R version 2.2.0 . This is our problem: None of the device drivers postscript() or pdf() is opening, e.g. the following message is coming: postscript("splines-Ex.ps") Fehler in postscript("splines-Ex.ps") : unable to start device PostScript

Hello, does anyone have experience with reading SPSS Version 15.0 files into R (version 2.4.1, WinXP)? I have long been sucessfully reading SPSS files with read.spss from the wonderful foreign package, but somehow after upgrading from SPSS14 to SPSS15 I seem to have problems. Trying a simple example, where test.sav is a SPSS 15.0 data file consisting of x1=c(1,2,3) and

Hi everyone, I want to make use of the save function but it did not work. I'm using vista and R 2.5.0, winzip is installed too. Here's the code (from example ?save): > x <- runif(20) > y <- list(a = 1, b = TRUE, c = "oops") > save(x, y, file = "xy.Rdata") Fehler in gzfile(file, "wb") : kann Verbindung nicht ?ffnen Zus?tzlich:

I cannot access internet through R. My IT-guys told me that I should set the proxy and the port and then everything will be fine. Where can I set them? I use R 2.6.2 under Windows 2000 and these settings are those of the Internet Explorer (no password). I read http://cran.at.r-project.org/bin/windows/base/rw-FAQ.html#The-Internet-download-functions-fail_002e but staring R by "C:\Program

R version 2.6.2 PowerBook G4 Hello R User, I try to perform an ANCOVA using the glm function. I have a dataset with continuous and categorical data (explanatory variables) and my response variable is also binary categorical. Fehler: NA/NaN/Inf in externem Funktionsaufruf (arg 4) Zus?tzlich: Warning messages: 1: In Ops.factor(y, mu) : - nicht sinnvoll f?r Faktoren (makes no sense for

I am looking for a way to get a warning message immediately after an evaluation within a function. To get error messages you can use geterrmessage(). But I found no function that allows me to check for warnings. Five years ago this questions has been posted but I haven't found any answer. Thanks for any help. Peter Wolf ------------ For illustration purpose a simple example follows: you

Full_Name: Christina Merz Version: R version 2.5.0 (2007-04-23) OS: mingw32 Submission from: (NULL) (213.70.209.132) R> version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status

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

Hi all, I can't come over a problem with formula. Suppose I have a coxmod model with the following formula: > somemod$formula Surv(lebzeit, tot == 1) ~ sex + (alter >= 65) + diff3k + zelltyp_k_c + q_nuc_3k + kar_k80_g80 + stadium and I want to drop the stadium explanatory variable from the model with > update(somemod, ". ~ . - stadium") I get the following messages:

Dear R-help list, I am working on an optimization with R by evaluating a likelihood function that contains lots of Gamma calculations (BGNBD: Hardie Fader Lee 2005 Management Science). Since I am forced to implement lower bounds for the four parameters included in the model, I chose the optim() function mith L-BFGS-B as method. But the likelihood often returns inf-values which L-BFGS-B

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