similar to: differences between survival models between STATA and R

  Hi r-users,   I would like to use gev and my data (annual rainfall ) is as follows:   > head(dat,20) A B C D E F G H I J 1 45.1 41.5 58.5 50.1 46.0 49.1 37.7 49.1 59.8 54.0 2 50.3 39.8 49.4 56.4 49.4 48.8 42.1 49.8 49.4 58.3 3 41.7 39.3 44.6 39.1 35.7 41.5 40.8 40.8 38.5 45.6 4 50.7 33.9 48.4 28.2 35.5 39.1 61.4 17.0 30.7 38.3 5 39.3 30.6 46.9 23.8 25.8

Hello, I have the following data.frame structure(list(Study = structure(c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L,

Hi, Suppose I have the following vector (data points): > x [1] 36.0 57.3 73.3 92.0 300.4 80.9 19.8 31.4 85.8 44.9 24.6 48.0 [13] 28.0 38.3 85.2 103.6 154.4 128.5 38.3 72.4 122.7 123.1 41.8 21.7 [25] 143.6 120.2 46.6 29.2 44.8 25.0 57.3 96.4 29.4 62.9 66.4 30.0 [37] 24.1 14.8 56.6 102.4 117.5 90.4 37.2 79.6 27.8 17.1 26.6 16.3 [49] 41.4 48.9 24.1

Hi All, I have a small dataset named das (43 cases) in which I am trying to create a binary outcome (1/0) based on the following code: if (das$age>65 && das$bmi>30) {das$danger<-1} else das$danger<-0 I am setting a flag called 'danger' to 1 of the subject is over 65 and has a BMI > 30. I find that my statement evaluates the first record in the data.frame and

I'm trying to read a comma delimited dataset that uses '.' for NA. I found that if the last field on a line was a missing '.' it was not read as NA, but just a '.', and the life variable was made a factor. The data looks like this, income,imr,region,oilexprt,imr80,gnp80,life Afghanistan,75,400.0,4,0,185.0,.,37.5 Algeria,400,86.3,2,1,20.5,1920,50.7

Dear R list members, I'm encountering a problem with expand.model.frame(): Suppose that I define the following simple function (meant just to illustrate the problem): > fun <- function(model){ + expand.model.frame(model, all.vars(formula(model))) + } > and I have the following model, created with an explicit data argument: > mod Call:

Hi, I want to remove the entire row from the dataset if any of the row contains blank or NA my dataset look like State Year Y X2 X3 X4 X5 X6 State1 1960 27.8 397.5 42.2 50.7 78.3 65.8 State2 1960 29.9 413.3 38.1 52 79.2 66.9 State1 1961 29.8 439.2 40.3 54 79.2 67.8 State2 1961 30.8 459.7 39.5 79.2 69.6 State1 1962 31.2 492.9 37.3 54.7 77.4 68.7 State2 1962 528.6 38.1 80.2 73.6 State1 1963

Dear all, I am trying to create a consensus skyline plot using the "ape" package(newbie). I have a nexus file that looks like the one from above containing many trees. (my_file.trees) #NEXUS begin trees; [Treefile generated by sim_coal.exe (Laurent Excoffier)] tree true_tree_1 = [&U] (((13.1:6, 43.1:5):12, 28.1:14):284, (((((21.1:0, (20.1:1, 4.1:0):0):1, ((37.1:0,

Dear useRs, i want to split a matrix having 1116rows and 12 columns. i want to split that matrix into 36 small matrices each having 12 columns and 31 rows. The big matrix should be splitted row wise. which means that the first small matrix should copy values which are in first 31 rows and 12 columns of the big matrix. similarly 2nd small matrix should contain values from 32nd to 63rd row of the

Hello, I teach statistics and use R Commander for teaching. I have 2 students out of 169 that can't get scatterplots or plot to work. I have had them update packages, restart R/R Commander/their computers and even reinstall R/R Commander. One is using Windows 7 on a new pc and the other is a pc user (not sure the OS). They are both using R2.11.1 and R Commander 1.6-0. The data look like

Hello, dear R-community. This is a question about TukeyHSD between factor combinations of a Three-Way ANOVA, which is - since it is a multi measure ANOVA - not a simple ANOVA but a Generalised Linear Mixed Model (GLMM), calculated with "lmer". > growth <- groupedData(length~meas|box_id,outer=~spec*comp*water,data=all.spec) > model <-

Hello all! I've been struggling with is for many hours today; I'm close to getting what I want, but not close enough... I have a dataframe consisting of two date-time columns followed by two numeric columns. what I need is the max value (in the first numeric column) based on the 2nd date-time column, which is essentially a factor. But, I want the result to provide both date-time values

Hi Jeff, Thank you for the suggestions -- I appreciate your help. Unfortunately, the result2 has two problems... (1) there are now 3 date columns (it looks like 2 cols are merged into 1 col) (2) the output rows should not have any of the basistime dates repeated (maybe I misstated the problem); I need the max fcst value by basistime, but also list the date value for that row; for example:

Hi, my question is really precise. I'm interested in drops of stomacal temperature (of seal pups) following milk ingestion. Temperature starts at an average of 37.8oC (sd=0.2) and drops rapidly to around 36.0oC and then stabilizes slowly to a temperature silimar to the previous one. It looks a bit like an inverse skewed bell curve. I want to find different points on that curve, and I have a

Colleagues, I am fitting an Emax model using nls. The code is: START <- list(EMAX=INITEMAX, EFFECT=INITEFFECT, C50=INITC50) CONTROL <- list(maxiter=1000, warnOnly=T) #FORMULA <- as.formula(YVAR ~ EMAX - EFFECT * XVAR^GAMMA / (XVAR^GAMMA + C50^GAMMA)) ## alternate version of formula FORMULA <- as.formula(YVAR ~ EMAX - EFFECT / (1 + (C50/XVAR)^GAMMA)) FIT <-

Hi, I have generated a skyline plot of a tree in newick format using "ape". How can I plot the error bars for this graph? I only have the the tree data. un<-"((((8.1:0, 20.1:0):0, 6.1:3):123, (((((35.1:0, (22.1:0, (43.1:1, 29.1:0):0):0):4, 25.1:6):0, ((42.1:0, 21.1:0):3, (39.1:0, 2.1:0):3):1):8, (3.1:0, 7.1:0):8):48, (((((15.1:0, 14.1:0):0, 11.1:0):0, 37.1:0):0,

Hi Team Can someone advice me on how i can lower the load average on my asterisk server? dahdi-linux-2.1.0.4 dahdi-tools-2.1.0.2 libpri-1.4.10.1 asterisk-1.4.25.1 2 X TE412P Digium cards on ISDN PRI Im using the system as an IVR without any transcoding or bridging ************************************** top - 10:27:57 up 199 days, 5:18, 2 users, load average: 67.75, 62.55, 55.75 Tasks: 149

Hello all, I'm trying to create a 2x2 matrix, 32 times after unlist() so that I can convert the list to matrix. I've looked through the R archive but couldn't find the answer. There is what I've done. > f <- system("ls *.txt", intern=TRUE) > x <- lapply(f, read.table) > x [[1]] V1 V2 1 -27.3 14.4 2 29.0 -38.1 [[2]] V1 V2 1 14.4

Hi Everyone, I have a dataset like this mean sd 0% 25% 50% 75% 100% n BODY TEMPERATURE 36.41099 0.4015699 35.1 36.22222 36.5 36.66667 37.1 89 DIASTOLIC BLOOD PRESSURE 73.60079 9.4656186 50.0 67.00000 73.0 80.00000 95.0 253 HEIGHT 171.94000 9.2011670 153.5 166.50000 173.0 176.25000 190.0

Dear group members, I want to compare response variables ("logAUC") of two groups (treatment "Test", "Reference") of a subset ("period == 1") in dataframe "resp" (below): sequence subject period treatment AUC logAUC 1 RT 1 1 Reference 44.1 3.786460 2 RT 1 2 Test 39.1 3.666122 3 RT 2