similar to: standard errors of fitted values are different S-plus survival pa ckage and R

I wonder if there is a package that will estimate a Zero Inflated Poisson Model (ZIP), and also if there is a package that will estimate what is called the Tobit model: that is a combination of censored and observed values in the same sample. Georgina Bermann Biostatistics AstraZeneca R&D M?lndal -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing

We are presenty using SAMBA 2.2 w. Windows 2000 sp1 and will be upgrading to Windows 2000 sp3. Are there any known or suspected problems with the combination of Windows 2000 sp3 and SAMBA 2.2. We are using Solaris 7 on the Unix side. /ola Ola Engstr?m Technical Computing & Information Services AstraZeneca R&D M?lndal S-431 83 M?lndal Sweden

If I have a set of data comprising a list of numbers of eggs on mosquito guts that range from 1 to 157. How can I get R to draw a barchart of the distribution of the data (i.e. x axis= number of eggs on a gut, y axis=number of mosquitoes found with that number of eggs)? Data is listed below. Thanks, Georgina Number of eggs on each mosquito gut: [1] 1 1 1 1 4 2 1 12 12 1 2

Hi! I have the data in a tab delimited text file titled "ken_data_try_anova." I tried to import it into R entering > read.delim(ken_data_try_anova) but received the error message Error in read.table(file = file, header = header, sep = sep, quote = quote, : object 'ken_data_try_anova' not found I have another file called 10423nad.txt. I tried

Does anyone know how to retrieve the R squared value for a line on a graph? Thanks Georgina

Can anyone advise a good transformation for this data below to produce a normalised distribution? Many thanks, Georgina No. of eggs on mosquito gut: 1 12 1 12 6 17 54 1 12 2 22 27 1 27 1 1 6 24 10 54 12 5 27 68 1 4 6 27 1 1 1 1 68 1 7 1 10 5 4 1 7 9 3 19 22 10 4 PhD Student Division of Infection and Immunity B5-29, GBRC 120 University Place Glasgow G12 8TA Tel: 0141 330

I have the attached data set (csv) and I want to run an analysis of variance on the wingsize data (comparing infected vs non-infected) within and between experiments. Can anyone help me with the command I should use? Many thanks Georgina PhD Student Division of Infection and Immunity B5-29, GBRC 120 University Place Glasgow G12 8TA Tel: 0141 330 5650

I have a set of data that comprises genome numbers in single eggs from three different parasite clones - 3D7, HB3, and MIX. I can draw a boxplot of the genome numbers for each clonefed but how do I carry out a t test or ANOVA to compare if the means are signifcantly different? (Data is listed below) Many thanks, Georgina Humphreys clonefed genomes HB3 21.3 HB3 23.5 HB3 25.9 3D7 27.2 HB3 28.1

The time has come to review the current Contributor and Core contributor grants for ZFS. Since all of the ZFS core contributors grants are set to expire on 02-24-2009 we need to renew the members that are still contributing at core contributor levels. We should also add some new members to both Contributor and Core contributor levels. First the current list of Core contributors: Bill

I have recently started a new job at an NHS hospital in Scotland. Since I took up this post 6 months ago I have had an ongoing dispute with the IT secutiry dept. who refuse to install R on my computer. I previously worked in another branch of the NHS where R was widely used and yet there is nothing I can say which will persuade the IT dept here to even visit the website! With some help from our

Hi! I am trying to work out the code to get a Fibonacci sequence, using the while() loop and only one variable. And I can't figure it out. Fibonacci<-c(1,1) while (max(Fibonacci)<500){ Fibonacci<-c(Fibonacci, (max(Fibonacci) + ?(Fibanacci))) } How can I tell R to take the value one before the max value? (Without defining another variable) (Probably super easy... I am a

Any ideas what is the problem with this code? > N <- 2; c(Sys.Date(), sprintf('N = %d', N)) [1] "2012-01-10" NA Warning message: In as.POSIXlt.Date(x) : NAs introduced by coercion Best regards, Ryszard Ryszard Czerminski AstraZeneca Pharmaceuticals LP 35 Gatehouse Drive Waltham, MA 02451 USA 781-839-4304 ryszard.czerminski@astrazeneca.com

Currently as the first column in a data frame I have string values in the format xx_yy - I want to create a new column with just the substring xx (for each row in turn). Three possible ways to do this might be (1) split the string by '_' using strsplit and paste the first of the resulting variables into a new column, but I have been unable to do this for each row of my data frame in turn

Is there any rational explanation for the bizarre seq() behavior below? > seq(2,8.1, lenght.out=3) [1] 2 3 4 5 6 7 8 > help(seq) > seq(2,8,length.out=3) [1] 2 5 8 > seq(2,8.1,length.out=3) [1] 2.00 5.05 8.10 Except maybe that it is early in the morning :) Best regards, Ryszard Ryszard Czerminski AstraZeneca Pharmaceuticals LP 35 Gatehouse Drive Waltham, MA 02451 USA 781-839-4304

Any ideas what is wrong? > strsplit("a.b", ".") # generates empty strings with split="." [[1]] [1] "" "" "" > strsplit("a b", " ") # seems to work fine with split=" ", and other characters... [[1]] [1] "a" "b" > > R.Version() $platform [1]

I am relatively new to R and am trying to do some basic data manipulation. Basically I have a table (csv - table 1) of data for a set of samples (rows), and a second table (table 2) of information about a subset of samples of particular interest. I want to pull out the data from table 1 for the samples in table 2, either by: * Merging the two tables based on a common identifier (SampleID - may

Is there an easy way to control smoothness of the contour lines? In the plot I am working on due to the undersampling the contour lines I am getting are jugged, but it is clear "by eye" these should be basically straight lines. In maps package I found smooth.map function, but maybe there is a more generic way of accomplishing the same thing. Ideally there would be an option to control

RF trains fine with X, but fails on prediction > library(randomForest) > chirps <- c(20,16.0,19.8,18.4,17.1,15.5,14.7,17.1,15.4,16.2,15,17.2,16,17,14.1) > temp <- c(88.6,71.6,93.3,84.3,80.6,75.2,69.7,82,69.4,83.3,78.6,82.6,80.6,83.5,76 .3) > X <- cbind(1,chirps) > rf <- randomForest(X, temp) > yp <- predict(rf, X) Error in predict.randomForest(rf, X) : subscript

This is what happens when i try to get acts_as_ferret ...."nothing much".... Please help me and excuse me if its really dumb, i''m new to this! thanks C:\rails\app>gem install ferret Attempting local installation of ''ferret'' Local gem file not found: ferret*.gem Attempting remote installation of ''ferret'' Updating Gem source index for:

I am currently summarising a data set by collapsing data based on common identifiers in a column. I am using the 'aggregate' function to summarise numeric columns, i.e. "aggregate(dat[,3], list(dat$gene), mean)". I also wish to summarise text columns e.g. by concatenating values in a comma separated list, but the aggregate function can only return scalar values and so something