similar to: t-distribution values

Dear all, what function can be used to calculate weighted SD or/and SE? There's no such stuff in 'base' package. Best regards, A.A.Morkovin

I have a several questions about R graphic properties. I use a "barplot2" function for creating plots with error bars. My data include species' names in first column, and I need make plots for each species. I know how to select species for each plot: D<-read.table("FD_R.txt", h=T) Dens<-D[D[,1]=="Sit.eur",] but I want to make a cycle which will

Dear all, is there any functions which allow to calculate Student t-criterion using means, their SE and sample size? I've seek for, but bulit-in t-criterion works only with initial sample... Best regards, A.Morkovin

Dear all i want to estimated the parameter of the gamma density(a,b,d) f(x) = (1/gamma(b)*(a^b)) * ((x-d)^(b-1)) * exp{-(x-d)/a)} for x>d f(x) = Age specific fertility rate x = age when i run this in R by usling nls() gamma.asfr <- formula(asfr ~ (((age-d)^(b-1))/((gamma(b))*(a^b)))* exp(-((age-d)/a))) gamma.asfr1 <- nls(gamma.asfr, data= asfr.aus, start = list(b = 28, a = 1, d= 0.5),

Dear All Subject : Multiple Lines in Graph Could you please help me to draw two lines in a graph. plot(f1, t1, type ="b") and plot(f2,t2, type="b") where t1, t2, f1,and f2 are single dimensional matrix. I have these two graph. How can I draw these two lines in a single window? Thanks in advance Kind Regards Wesley C Mathew [[alternative HTML version deleted]]

Hi list, One of my co-workers found this problem with 'cor' in his code and I confirm it too (see below). He's using R 2.2.1 under Win 2K and I'm using R 2.3.0 under Win XP. =========================================== > R.Version() $platform [1] "i386-pc-mingw32" $arch [1] "i386" $os [1] "mingw32" $system [1] "i386, mingw32" $status

Pat Burns makes a good point. -Peter -------- Original Message -------- Subject: Re: [R] Using seq_len() vs 1:n Date: Fri, 12 Feb 2010 09:01:20 +0000 From: Patrick Burns <pburns at pburns.seanet.com> To: Peter Ehlers <ehlers at ucalgary.ca> References: <4B746AEF.10900 at ucalgary.ca> If you want your code to be compatible with S+, then 'seq_len' isn't going to work.

#I wish to find the "next non-NA" value within each row of a data-frame. #e.g. I have a data frame mydata. Rows 1, 2 & 3 have soem NA values. mydata <- data.frame(matrix(seq(20*6), 20, 6)) mydata[1,3:5] <-  NA mydata[2,2:3] <-  NA mydata[2,5] <-  NA mydata[3,6] <-  NA mydata[1:3,] #this loop accomplishes the task; I am tryign toi learn a "better" way for(i

Dear Rexperts: (R 2.0.1 on Windows XP Pro) Is the following problem unique to my setup? If it's a known problem, I didn't see it at http://bugs.r-project.org/ nor find discussion in the archives. plot(1:10) loc <- c(5, 6) text(loc, labels = "a") Produces expected results according to ?xy.coords. plot(1:10) loc <- list(x = 5, y = 6) text(loc, labels = "a") No

Hello, I have a simple question that I could not really figure out. I am plotting labels within a graph using the text function. I first plot the first label by specifying the x and y coordinates on the graph. Then to plot the second label next to it, I am using te strwidth function to get the width of the previous label in user coordinate then add the maximum width to the x value. However, this

This is probably simple, but I have a hard time finding the solution. Any help greatly appreciated.   I would like to use the results of fitdistr(z,densfun=dweibull,start=list(scale=1,shape=1)) for further processing.  How do I assign the values of scale and shape to b and a without manually entering the numbers?   TIA __________________________________________________________________

x<-c(1,4,9,20,3,7) i want to get a serie c(5,13,29,23,10). y <- c() for (i in 2:length(x)){ y[i-1] <- x[i-1]+x[i]} is there more simple way to get? [[alternative HTML version deleted]]

Full_Name: Cornell Gonschior Version: 2.10.0 OS: Linux Submission from: (NULL) (212.201.28.40) Hi, in the introduction to R, you can find the following sentence in the par() chapter: "Use tck=0.01 and mgp=c(1,-1.5,0) for internal tick marks." I thought that's nice, because I wanted to have tick marks and tick labels inside and the axis title outside. But: > plot(z, las=1,

mytest<-function(x,f){ sum(x*f)/sum(f) } cat(mytest,file="c:/mytest.R") Error in cat(list(...), file, sep, fill, labels, append) : argument 1 (type 'closure') cannot be handled by 'cat' how can i write the mytest function into a file c:/mytest.R with cat function? [[alternative HTML version deleted]]

Hello. Consider the response-variable of data.frame df is constant, so analytically perfect fit of a linear model is expected. Fitting a regression line using lm result in residuals, slope and std.errors not exactly zero, which is acceptable in some way, but errorneous. But if you use summary.lm it shows inacceptable error propagation in the calculation of the t value and the corresponding

Readers, I tried to the following commands: plot(y~x,ylab=expression(A[1]~B[2],xlab=expression(C~D)) mtext(expression(A[1]~B[2]),"additional text",side=3,line=1) I receive the text that I want, but the command terminal shows the following response: Warning message: NAs introduced by coercion in: mtext(expression(A[1]~B[2]),"additional text", side = 3, What is my mistake

I vote to 'fortunize' Doug Bates on Hierarchical data sets: which software to use? "The widespread use of spreadsheets or SPSS data sets or SAS data sets which encourage the "single table with a gargantuan number of columns, most of which are missing data in most cases" approach to organization of longitudinal data is regrettable."

First my excuses if I keep bugging everyone in this list, but I am a newbie, and tend to find some behaviour that looks unexpected to me; and I would really appreciate to be pointed to some location that allows me to understand more about this software. Here is my next question: > par(mfrow=c(1,1)) > qqnorm(rnorm(20)) > qqmath(rnorm(20)) > par(mfrow=c(3,4)) > for(i in

Hi, I want to estimate parameters of weibull distribution. For this, I am using fitdistr() function in MASS package.But when I give fitdistr(c,"weibull") I get a Error as follows:- Error in optim(x = c(4L, 41L, 20L, 6L, 12L, 6L, 7L, 13L, 2L, 8L, 22L, : non-finite value supplied by optim Any help or suggestions are most welcomed -- View this message in context:

Hi List I have the next code and the error. I have try with other codes and I have the same problem. > reut21578 <- system.file("texts", "crude", package = "tm") > (r <- Corpus(DirSource(reut21578), readerControl = list(reader = > readReut21578XMLasPlain))) A corpus with 20 text documents > (r <- Corpus(DirSource(reut21578), readerControl =