similar to: count data with a specific range

Hi, I am thinking about IWLS vs ML estimation. When I use glm() for a 2-parameter distribution (e.g., Weibull), I can otain the MLE of scale parameter given shape parameter through IWLS. Because this scale parameter usually converges to the MLE. In this point, I am wondering: i) can you say that the direct MLE, which is obtained by maximizing a likelihood function, is equalvant to the indirect

Data has the first row for variable name and the first column for sample name. I want to take "Log" for all data, but how to compute without the first column for sample name. That is, column 1: sample ID column 2-10: data I want to find an average on each column (2-10) > apply(raw_data,2,mean) Error in Math.data.frame(list(sample_id = c(1L, 2L, 3L, 4L, 5L, 6L, 7L, :

Is there anyone knowing a function or way for standard normal cumulative distribution? ?(z=-0.1)=? also ?(z=?)=0.025 Thank you, -- View this message in context: http://www.nabble.com/calculation-for-standard-normal-cumulative-distribution-tp20282804p20282804.html Sent from the R help mailing list archive at Nabble.com.

Hi: I am coding some interactive interface using an xyplot from lattice. When a user clicks on a panel of the xyplot, a simulation is executed and the resulting data is plotted on the corresponding panel. The problem is that I start with an empty data frame and only fill it as the user request data from simulations. And, I have found it impossible to create a conditional plot that contains panels

I have an error for a simple optimization problem. Is there anyone knowing about this error? lambda1=-9 lambda2=-6 L<-function(a){ s2i2f<-(exp(-lambda1*(250^a)-lambda2*(275^a-250^a)) -exp(-lambda1*(250^a)-lambda2*(300^a-250^a))) logl<-log(s2i2f) return(-logl)} optim(1,L) Error in optim(1, L) : function cannot be evaluated at initial parameters Thank you in advance -- View this

Dear list, I would like to format the result of the 'cut' function to perform a subsequent frequency distribution table (fdt) suitable for publications. Below an reproducible example: set.seed(1) x <- c(rnorm(1e3, mean=10, sd=1), 50, 100) start <- 0 end <- 110 h <-10 c1 <- cut(x, br=seq(start, end, h), right=TRUE) levels(c1) # I get: # [1] "(0,10]"

Dear list, This weeks RCW has a slightly different format than before. After all, since what we want is less bloat, this RCW is less verbose. Please do let me know what you think off-list. *Interesting stuff on the Rails Timeline* _Fixed and commited to trunk:_ An issue with not null columns in SQLite3 migrations(<a

Hello, Considering 5 points in X-Y plain. Data is a 5*2 matrix (5 rows for samples & 2 columns for X and Y) With a distance from the origin, if a distance < 7, remove the row from the Data. After calculating the distance for each point, I can't forward because of this "Removing" problem. Anyone can help me? Data=matrix(1:10,5,2) a=rep(0,5); b=rep(0,5); for (i in

Hello, I have different results from these two softwares for a simple binomial GLM problem. >From Genmod in SAS: LogLikelihood=-4.75, coeff(intercept)=-3.59, coeff(x)=0.95 >From glm in R: LogLikelihood=-0.94, coeff(intercept)=-3.99, coeff(x)=1.36 Is there anyone tell me what I did wrong? Here are the code and results, 1) SAS Genmod: % r: # of failure % k: size of a risk set data

I use "while" loop but it produces an errro. I have no idea about this. Error in "[<-"(`*tmp*`, i, value = numeric(0)) : nothing to replace with The problem description is The likelihood includes two parameters to be estimated: lambda (=beta0+beta1*x) and alpha. The algorithm for the estimation is as following: 1) with alpha=0, estimate lambda (estimate beta0

V; Sent via BlackBerry from T-Mobile -----Original Message----- From: r-help-request at r-project.org Date: Tue, 30 Sep 2008 12:00:06 To: <r-help at r-project.org> Subject: R-help Digest, Vol 67, Issue 31 Send R-help mailing list submissions to r-help at r-project.org To subscribe or unsubscribe via the World Wide Web, visit https://stat.ethz.ch/mailman/listinfo/r-help or, via email,

Hi, I've found that the default options of gcc and g77 on alpha produce an executable that suffers from floating exceptions. There seems to be a problem with the configure scripts that prevents the --with-f77 option working correctly. The configure script tests that f77 is running by attempting to compile and link a dummy program - the link step seems broken on the alpha (it doesn't like

Should negative subscripts be outlawed in x[[subscript]] ? Currently, if subscript is a scalar then it can only be negative if length(x)==1 (otherwise [[ throws an error). If length(subscript)>1 then it gets treated as an attempt to recursively extract an element of a nested list. > list(10,20)[[-1]] # get the last element out of 2 [1] 20 > list(10,20,30)[[-(1:2)]] # get the last

Ok, I''ve flogged this issue probably longer than some of you can stand by now. (remember, I''m the nut trying to use a PPro200 to support ~500 users on a 10Mb internet link :o) To appease those who think I''m nuts, I am ordering a new firewall shortly to allow for future growth. (probably a Dell PE750 with P4/2.8 and dual GE nics.) However, since I have yet to prove

Hi All, I have a file with two columns, the first column has the names of the patients and the second column has the age. I am looking into creating an output file that looks like 1-10 10-20 etc Eric Chris Bob mat Andrew Suzan Where each column has the name of the patients in a given age category that is displayed in the header. For example in the output, the first

[Image: http://www.filmupdatelist.com/Saw.3D-MAIN.jpg ] Saw.3D.2010.CAM.READ.NFO.XViD-LKRG Code: Saw 3D (2010) DVDRip XvID AC3| 789 MB As a deadly director rages over Jigsaw's melton legacy, a group of Jigsaw greutert gathers to seek the support of self-help guru and fellow survivor Bobby Dagen, a man whose own dark secrets unleash a new wave of terror. Director: Kevin Greutert

I made this rather cool plot which I am quite pleased with: http://brainimaging.waisman.wisc.edu/~perlman/data/BeeswarmLinesDemo.pdf However, I feel there must be a better way to do it than what I did. I'm attaching the code to create it, which downloads the data by http so it should run for you if you have the current version of beeswarm installed (which was just updated today,

Jeff, Perhaps I was sloppy with my notation: I want groups >=0 <10 >=10 <20 >=20<30 ...... >=90 <100 In any event, my question remains, why did the four different versions of cut give me the same results? I hope someone can explain to me the function of include.lowest and right in the call to cut. As demonstrated in my example below, the parameters do not seem to alter

Found the solution to my 2.0.6 dying on me on SGI - SYSV_IPC needed to be #undef'ed in source/include/includes.h in favor of MMAP, seems to have fixed it. If I run the daemon in server mode, no prob. If I switch to domain mode (and use smbpasswd to join into the domain, per the instructions), I get the error indicating that the username/password pair was incorrect, as if I typed in a wrong

Hello, My question sounds simple, but as I am desperatly searchin for a solution I am asking you all. :) I try to filter out all non-numeric columns of a data frame using a for loop to go through all columns. My if clause looks like this:         for(j in 1:length(data)) {                 ...                 if(!is.numeric(data[j])) skip <- "not numeric")                 ...