similar to: Problem with generated parameter estimates

Dear List: I am running into a memory issue that I haven't noticed before. I am running a simulation with all of the code used below. I have increased my memory to 712mb and have a total of 1 gb on my machine. What appears to be happening is I run a simulation where I create 1,000 datasets with a sample size of 100. I then run each dataset through a gls and obtain some estimates. This works

Hi R-helpers, My real data is a panel (unbalanced and with gaps in years) of thousands of firms, by year and industry, and with financial information (variables X, Y, Z, for example), the number of firms by year and industry is not always equal, the number of years by industry is not always equal. #reproducible example firm1<-sort(rep(1:10,5),decreasing=F) year1<-rep(2000:2004,10)

Dear all I need some help with R. How can I save the estimates of a regression model in a file? here is what I did: 1) this is my regression model: fit1 <- lm(logmilk ~ logdays + days, data=data2) 2) however, I want to get the parameters estimates for each individual (by group): so i did the following: by(data2, list(data2$V2),function(.data2) lm(logmilk ~ logdays + days, data= .data2)) 3)

Hello, I am trying to figure out how to use optimize() with array variables as inputs. I have a for loop in the function definition: SS <- function(int,slo,x,y){ for(i in 1:length(x)) ((int+slo*x[i])-y[i])^2->squares[i] sum(squares)->>sum_squares output_txt = c ("The sum of squares is", sum_squares) print(output_txt, quote=FALSE)} Even assuming I make x and y

I can't create new folder in a czech Win95 explorer on a unix using samba 1.9.18p7. The problem is that Windows exporer uses for the creation new folder a name containing characters above 128. The explorer creates new folder in a two phases: the first It creates directory named Nov<A0> slo<A7>zka and next it gives to user posibility to change this name. Crazy !!? The first phase

Two friend reported me a problem, which I can't solve: (I run R-1.0.0, Debian Linux) They hava a function "corr.matrix" (see end of mail), and when they create a 173x173 matrix with this function V <- corr.matrix(0.3, n=173) V1 <- qr.solve(V) reports: Error in qr(a, tol = tol) : NA/NaN/Inf in foreign function call (arg 1) For n < 173, qr.solve returns the correct

hi, - big and ugly mail. If you don't like them, delete it now :-) - I have collected and classified some information of: http://home.earthlink.net/~rwhron/kernel/bigbox.html And I observed that ext3 performance is worse than previous kernels(2.4.19...). -ac and -aa are here only as reference. Complete information is in the upper URL. dbench: Performance is worse. dbench (Numbers are in

Hello! I am trying to put two math expressions in the title of a plot. As you can see below, I can place correctly one expression at a time, but not both. Ideally I would like to have them separated by a comma. Any suggestions? > k <- 1 > n.eff <- c(20, 30) > ### this works > plot(0,0, main = substitute(n == k, list(k = k))) > ### this works > plot(0,0, main =

Hello List: Please find uploaded to CRAN a new package, PwrGSD The package is intended for the design and analysis of group sequential trials There are two main functions, (1) GrpSeqBnds: computes group sequential stopping boundaries for interim analysis of a sequential trial based upon a normally distributed test statistic. This can be done via the Lan-Demets procedure with

Hello List: Please find uploaded to CRAN a new package, PwrGSD The package is intended for the design and analysis of group sequential trials There are two main functions, (1) GrpSeqBnds: computes group sequential stopping boundaries for interim analysis of a sequential trial based upon a normally distributed test statistic. This can be done via the Lan-Demets procedure with

Hi all: I have two curve models: model1<-nls(result ~ exp(b0 + b1*(time)), start = list(b0 = 0, b1 = 5),trace=TRUE,data=data1) model2<-nls(result ~ exp(b0 + b1*(time)), start = list(b0 = 0, b1 = 5),trace=TRUE,data=data2) I wanna compare the two models to find out whether the difference between them is significant or not. How can I do then? Many thanks! My best [[alternative

Hello list, I've seen the recent discussions documenting problems with lm. I have encountered the following problem. I use WinXP Pro with service pack 1, and R 1.9.0, on a XEON 2GHz, with 1GB of RAM. > eff.fro std.dev mean NSTRDSP 7.403749e-01 1.215686e-01 CPFGEP 9.056763e+00 1.815686e+00 WSWOLF 4.703588e+05 1.112832e+05 NPILGRIM 1.017640e+06 2.134335e+05

Dear all, This problem came up initially while debugging a function, but it seems to be a more general problem of R. I hope I'm wrong, but I can't find another explanation. Let me illustrate with the raster package. For an object "RasterLayer" (which inherits from Raster), there is a method xyValues defined with the signature

I have a large dataset grouped by a factor and I want to perform a regression on each data subset based on this factor. There are many ways to do this, posted here and elsewhere. I have tried several. However I found one method posted on the R wiki which works exactly as I want, and I like the elegance and simplicity of the solution, but I don't understand how it works. Its all in the formula

I have 3 vectors x,y,z: x<-c(-2.88,-1.92,-.96,0,.96,1.92,2.88,3.84, -1.92,-.96,0,.96,1.92,2.88,3.84, -.96,0,.96,1.92,2.88,3.84, 0,.96,1.92,2.88,3.84, .96,1.92,2.88,3.84, 1.92,2.88,3.84, 2.88,3.84, 3.84) y<-c(rep(-3.84,8), rep(-2.88,7), rep(-1.92,6), rep(-.96,5), rep(0,4), rep(.96,3), rep(1.92,2), 2.88) z<-c(.65,1.78,2.4,2.54,3.04,2.22,2.97,3.56, .745,1.287,1.98,1.91,2.02,2.24,2.95,

In the course of exploring response prediction, I stumbled upon a small discrepancy between the CIs produced by predict.lm() and all.effects() require(mlmRev) require(effects) hsb.lm <- lm(mAch ~ minrty * sector, Hsb82) hsb.new <- data.frame( minrty = rep(c('No', 'Yes'), 2), sector = rep(c('Public', 'Catholic'), each = 2)) hsb.eff <-

R-helpers, I am hoping to find someone who uses both R and program Stata for GLMs. I am a beginner R user, finding my own way through; learning code etc. at the same time as learning the statistics I need to complete my project. What I have is the code from Stata and am trying to reproduce the same analysis in R - my program of choice. . glm count md ms rf sg, family(poisson)

Hi All, I am trying to use optim() to minimize a function with a penalty function term. This is a simple model bioeconomic model of a fishery. The penalty function constrains the amount of effort (f) at 9. This works fine. The code is: ********** nfleets<-2 M<-1 M<-array(M,dim=c(nfleets)) N<-1000 cost<-c(30,30) cost<-array(cost,dim=c(nfleets)) Price<-2

Hi, I have installed samba 2.0.6 under HP-UX 10.20. I downloaded the precompiler binaries for samba 2.0.6 and followed as per the instructions and installed the SAMBA from the HP-UX depot it created. And now when I start from /sbin/init.d/samba the smbd daemon starts and this is the message I see in the /var/opt/samba/log.smb : invalid argument

After an upgrade to Asterisk 12, I'm "collecting" channels. When I enter and then exit a conference room, I see: -- <CBAnn/207-0000067f;1> Playing 'confbridge-leave.slin' (language 'en') -- Channel CBAnn/207-0000067f;2 joined 'softmix' base-bridge <5edb1920-3774-4ba3-8c4d-23e8fd04519c> -- Channel CBAnn/207-0000067f;2 left