similar to: 2.13 version doesn't load packages

HI there, I have been trying to use a code posted on R help to be able to calculate area under the curve for complicated data points and there seems to be an issue with the code: no "b" object found. I am not a good R user and can''t find were the problem is. Any help? Thanks!! This is the code ( as a test run I gave it this info because I know the answer: x<-seq(1:50)

Dear bugs @r; The following reproduces erroneous results on my system using R-gui: platform i386-pc-mingw32 arch x86 os Win32 system x86, Win32 status major 1 minor 4.1 year 2002 month 01 day 30 language R ##################################################### inter.gp.effect <- 0.1 animals <- data.frame(pos = 0:2) SI.model <- function(Animal.group) {

Hi there: I am a total beginner in R, and I have a simple question: I have a table with thousands of lines that represent locations, and two columns: latitude and longitude. I need to randomly sample 1000 lines. How do I do it? I know the command "sample", but it samples elements independently, not lines. If there is a better place for me to ask that type of question, please let me

Hi, everybody. Sorry to bring up the subject again but I have visited the revolution computing web page, but its not clear when or what will be the new release be. Does anybody have information about that?. Does anybody know if the version that will be released will be validated?. I have been reading the previous mails about validation and I work in a Pharma company and validation is a big

Hi all, there is a vector v with several NAs. I want to create a new vector n of the same length as v and the same NAs as in v and tried this: n <- vector(length=length(v), mode="numeric") replace(n, which(is.na(v)), NA) but this does't work, all values in n are 0. What went wrong? Thanks, Sven -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-

Hola que tal alguien de ustedes me podría auxiliar, lo que pasa es que estoy tratando de hacer un programa en R utilizando la siguiente tabla: x<1 km x>1km Estabilidad a c d f c d f A 213 440.8 1.94 9.27 459.70 2.094 -9.6 B 156 106.6 1.15 3.30 108.20 1.098 2.0 C 104 61.0 0.91 0.00 61.00 0.911 0.0 D 68 33.2 0.73 -1.70 44.50 0.516 -13.0 E 51 22.8 0.68 -1.30 55.40 0.305

Dear all, I don't understand the following behaviour: Running compar.gee (in library ape ) with and without the option 'data', it give me different results Example: .... Start R .... > load("eiber.RData") > ls() [1] "gee.na" "mydata" "mytree" > library(ape) > # runnig with the option data= mydata > compar.gee(alt ~ R,

Hi Paul: I am from Brazil and was looking for help to solve a memory problem using R. I found your answer to Melinda on https://stat.ethz.ch/pipermail/r-devel/2003-July/027094.html. You say there that you are not familiar with what to do when using Windows. I'm using Windows and would like to know if Melinda found a way to solve her problem. Can you give me her e-mail address?? Thanks...

Hello, I am trying to run the APE program COMPAR.GEE with a model containing a categorical response variable and a mixture of continuous and categorical independent variables. The model runs when I have categorical (binary) response and two continuous independent variables (VAR1 and VAR2), but when I include a categorical (binary) independent variable (VAR3), I receive the following output with

Tengo un archivo de texto donde cada línea es de la siguiente forma: "aa-mm-dd hh:mm:ss Nombre Apellido" ¿Hay alguna forma de usar read.table o algo similar para obtener directamente un data frame que tenga dos columnas donde una tenga la fecha y quede de tipo PosiX y la otra character con el nombre completo? Gracias, Sebastián.

Dear Listers, I am trying to install packages via install.packages() from MacOS 10.3.8. Installing work fine when run from the menu, but the following command (useful for setting up each computer of the student computer room) leads nowhere for some reasons: >

Greetings, I have a very, very, simple research question. I want to predict one dichotomous variable using another dichotomous variable. Straightforward, right? The issue is that the dataset has two issues causing some complications for me. 1) The subjects are not independent -- they are sibling pairs. Every person in the dataset has a sibling in the dataset. This needs to be treated a

Dear all, I need to fit a gee model with an auto-regressive correlation structure and I faced some problems. I attach a simple example: ####################################################### library(gee) library(geepack) # I SIMULATE DATA FROM POISSON DISTRIBUTION, 10 OBS FOR EACH OF 50 GROUPS set.seed(1) y <- rpois(500,50) x <- rnorm(500) id <- rep(1:50,each=10) # EXAMPLES FOR

Dear List: I found an error when I called the 'gee' function. I cannot solve and explain it. There are no errors when I used the 'geeglm' function. Both functions fit the gee model. The project supervisor recommends me to use the 'gee' function. But I cannot explain to him why this error happens. Would you help me solve this problem? I appreciate your help. In

Hi, I am trying to install gee on our server but I get the error below. I do not have root on this machine so no control on how R was installed itself. It looks like it cannot find blas libs, the only ones i can find on the machine are: /usr/lib64/libblas.so.3 -> libblas.so.3.0.3 /usr/lib64/libblas.so.3.0 -> libblas.so.3.0.3 /usr/lib64/libblas.so.3.0.3 and : $ R CMD config BLAS_LIBS

Hello, I am running the program COMPAR.GEE within the package APE. My dependent variable is binomial, while my independent variable is a multi-state categorical variable. The output reports an estimate for each state of the independent variable except the first one. For example, for the variable X with 3 states, the output is: intercept (estimate) X2 (estimate) X3 (estimate)

Hi, my dependent variable is a proportion ("prob.bind"), and the independent variables are factors for group membership ("group") and a covariate ("capacity"). I am interested in the effects of group, capacity, and their interaction. Each subject is observed on all (4) levels of capacity (I use capacity as a covariate because the effect of this variable is normatively

windows Vista R 2.10.1 Is it possible to get p values from gee? Summary(geemodel) does not appear to produce p values.: > fit4<- gee(y~time, id=Subject, data=data.frame(data)) Beginning Cgee S-function, @(#) geeformula.q 4.13 98/01/27 running glm to get initial regression estimate (Intercept) time 1.1215614 0.8504413 > summary(fit4) GEE: GENERALIZED LINEAR MODELS FOR

Hi, I downloaded both gee and geepack, and I am trying to understand the differences between the two libraries. I used the same data and estimated the same model, with a correlation structure autoregressive of order 1. Surprisingly for me, I found very different results. Coefficients are slightly different in value but sometimes opposite in sign. Moreover, the estimate of rho (correlation

Does anyone know how could I get the significant codes from mixed model output fitted with a GEE package? The output I got is the following: GEE: GENERALIZED LINEAR MODELS FOR DEPENDENT DATA gee S-function, version 4.13 modified 98/01/27 (1998) Model: Link: Logit Variance to Mean Relation: Binomial Correlation Structure: Exchangeable Call: gee(formula = bru