similar to: combination of different vector values

Hello all, I have a very simple set of data and I would like to analyze them with probit analysis. The data are: X Event Trial 100 8 8 75 8 8 50 6 8 25 4 8 10 2 8 0 0 8 I want to estimate the value of X that will give a 95% hit rate (Event/Trial) and the corresponding 95% CI. Anyone can offer some help? Thanks!! -

Hello everyone, Does anyone know if I can do Breslow Day Test for the homogeneity of odds ratio in R? Thanks! - Jacqueline

Hi all, I am trying to run the k-means cluster analysis using the function kmeans in the package cluster. The data are: x = c(-0.26, -0.23, -0.05, -0.20, 0.30, -0.84, -0.10, -0.12, 0.10, -0.31, -0.19, 0.18, -0.26, -0.23, -0.37, -0.23) I've got two different solutions when I ran this function over a few times: kmeans(x, centers=2) The first solution gives the following: $cluster [1]

Dear list, I have a simple question. If I need to run linear regression on a subset of the dataset, will the following two commands give the same output: lm(y ~ as.factor(x1) + as.factor(x2), data, subset = z==Z[1]) and lm(y ~ as.factor(x1) + as.factor(x2), data[data$z==Z[1],] ) I have got different results running the two commands. Am I missing something here, or it there a bug on the

Hi all, Suppose I have a vector Fruits <- c('Apple','Orange','Pear','Banana','Mango') I want to print the statement "The fruits in the basket are: Apple, Orange, Pear, Banana, Mango" If I do: print(paste('The fruits in the basket are:',Fruits)) The output is: [1] "The fruits in the basket are: Apple" "The fruits in

How does R compute the p-value in the Wilcoxon Rank Sum Test? If I have the test statistics, can I get the p-value using the function pwilcox? Thanks -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the

Dear R helpers, I have a 3D scatter plot that I have generated from scatterplot3d (which looks great- thanks!) and I can see that the points in my graph fall in a plane. Following the example 5 from 3D scatter plot (below) I have fitted a regression plane. Now what I would like to do is a rotation so that my new co-ordinate system is about the fitted plane (by finding the normal to the plane

Dear R users, I'm interested in measuring the stability of a heirarchical clustering, of the overall clustering and finding sub clusters (from cutting the heirarchical clustering at different levels) which demonstrate stability. I saw some postings on the R help from a while back about bootstrapping for clustering (using sample and generating a consesus tree with a web based tool CONSENSE)

Dear list, I want to create a POSIX time vector as follows: day <- as.character("110809") time.t <- 1:3600 t.min <- time.t %/% 60 t.sec <- time.t-t.min*60 DATE <- as.POSIXct(strptime(paste(day,t.min,t.sec),"%y%m%d %M%S")) Tail(DATE) The problem is that the last element (3600) returns a NA and I don't understand why. 600, 1200, 2400

Dear all, I have a question on the dose-response estimation with clustered/ correlated binary data. I would like to estimate the hit rate for a certain test at various concentration levels. The test is used on 5 subjects, and each subject is tested 20 times. If we assume that the 100 samples are independent, the hit rate estimate is unbiased, but the variance is under-estimated. The other

Hi all, Does anyone know if there is an implementation of Brookmeyer & Crowley's confidence interval for the median survival time in R? Reference : Brookmeyer & Crowley, "A confidence interval for the median survival time" (1982) Biometircs Thanks in advance for your help, Jacqueline [[alternative HTML version deleted]]

Dear R list, I have a package downloaded and looked up a function in there. Now I find that it uses C code (.C call) to do part of its job. However, when I wanted to look that part up by using getAnywhere() I was told that no object of that name could be found. Then I tried typing C.(" funct.name") which only caused R to crash. Is there any way to look up the function? I'm not

Dear list, I have problem that I cannot solve and would like to ask your opinion. I tried to ask a few days ago already but got no answer and all my attempts to solve it by myself since then failed. Sorry for repeated posting! Here the problem broken down a bit. My problem basically is, that I want to use the elements of a character vector as names for objects and by recalling only the

> -----Original Message----- > From: Jacqueline Lee [mailto:jlee@isdomaininc.com] > Sent: Friday, September 23, 2005 11:46 AM > To: asterisk-users@lists.digium.com > Subject: channel offhook state > > > We are using a digium card (TDM400) with asterisk for our access to the > PSTN. Initially when the server starts, all the zap channels on the card > are in the

I need help with replacing NaN with zero (the value '0') in my dataset. The reason is that I can't get it to graph because of the NaN in the dataset. I have tried: data[is.nan(data)] <- 0 that others have suggested in the help archives but this does nothing so I am not sure what I am doing wrong. Thanks in advance for your help. [[alternative HTML version deleted]]

Hi all I have a multiple datasets that look like this MM Freq 1 30 2 35 3 54 4 33 5 27 6 13 7 25 8 29 9 40 10 32 11 36 12 23 I am plotting this using barchart (there is probably something better but it gives me the results I want) and I would like the x-axis to have the names of the months instead of the month numbers. So I have

I have a data set that I have imported (not sure if that makes a difference) and I would like to calculate the sum of only specific columns. I have tried >colSums(dataset, by=list(dataset$col5), dims=1) and I get an error of unused arguments I have also tried >aggregate(dataset, by=list(dataset$col5), sum) and I get the error that sum is not meaningful for factors. I want to only calculate

Hi all I have a bit of a problem. I want to make a barplot of some data. My data is of a score that is separated by year and by a limit (above 3 and below 3 to calculate the score). Year Limit HSS 1999 ALT 0.675 1999 VFR 0.521 2000 ALT 0.264 2000 VFR 0.295 I would like to have a barplot with year on the x axis and HSS on the y axis and the two limits as two different colors to show the

Hello R users I have a dataset that has different types of records with different dates and times pertaining to each. I would like to have a bar graph of a count of the types(ie. The number of types) of recods by hour grouped by year. So the count of the types would be the y axis, the hour on the x axis and then grouped by year for easy comparison. I think that I have to use barchart however I

I am using barchart to make charts for some data with a lot more functions and labels and such in the command. barchart(Freq ~ factor(HH), data = dataset1, group= year) So I have my data grouped by year and I get a legend at the top of graph, which is great cause I need the legend for the different years but it is a weird spot. So how can I manipulate the legend, ie. Move it, shrink it, do