similar to: nlme Error: Subscript out of bounds

Hello. A quick question about incorporating variation due to study in the metafor package. I'm working with a particular data set for meta-analysis where some studies have multiple measurements. Others do not. So, let's say the effect I'm looking at is response to two different kinds of drug treatment - let's call their effect sizes T1 and T2. Some studies have multiple

Hello, I have a deriv function that I am feeding to nlme. It works, and I can use it in nls, but when I try to use it in nlme I get Error: subscript out of bounds. I can fit the model using SSasympOrig, instead of the deriv function, but I am trying to reproduce an earlier analysis (done in Splus) and I get slightly different results with SSasympOrig. Here are my calls: This does fit

Hi all, I post it on both r-help and r-finance since I don't know where is most appropriate for this topic. Sorry if it bothers you. I did garch fitting on S&P500 monthly returns with garchFit from fSeries. I got same coefficients from all cond.dist except normal. I thought that is probabaly usual for the data. But when I play with it, I got another question. I plot skew normal with

Hi. Background - I am working with a dataset involving around 750K observations, where many of the variables (8/11) are unordered factors. The typical model used to model this relationship in the literature has been a simple linear additive model, but this is rejected out of hand by the data. I was asked to model this via kernel methods, but first wanted to play with the parametric

Dear friends, After running the lm() model, we can get summary resluts like the following: Coefficients: Estimate Std. Error t value Pr(>|t|) x1 0.11562 0.10994 1.052 0.2957 x2 -0.13879 0.09674 -1.435 0.1548 x3 0.01051 0.09862 0.107 0.9153 x4 0.14183 0.08471 1.674 0.0975 . x5 0.18995 0.10482 1.812 0.0732 . x6 0.24832 0.10059 2.469 0.0154 * x7

Hello, I am trying to perform a logistic regression using counts. For example: cedegren <- read.table("http://www.cloudstat.org/index.php?do=/attachment/download/id_95 /", header=T) attach(cedegren) ced.del <- cbind(sDel, sNoDel) ced.logr <- glm(ced.del ~ cat + follows + factor(class), family=binomial("logit")) This works. However, if I change the family to

Hi all, I was trying garchFIt() of fSeries to fit volatility of monthly log returns of S&P500. I tried residuals of normal, student t, skew normal, skew t. But all innovations except normal got exaxtly same coefficients, even if I changed their parameters of skew and shape. Is this correct for the data or something wrong? I am attaching the code, thank you. Muster #GARCH analysis of

Hello, I have been looking for functions for calculating the within-groups variance and between-groups variance, for the case where you have several numerical variables describing samples from a number of groups. I didn't find such functions in R, so wrote my own versions myself (see below). I can calculate the within- and between-groups variance for the Sepal.length variable (iris[1]) in

The following behaviour (in R 0.62.3) is disturbing: > (0+0i)^2 [1] NaN+NaNi Is it deliberate?? Laimonis Kavalieris -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To:

The following behaviour (in R 0.62.3) is disturbing: > (0+0i)^2 [1] NaN+NaNi Is it deliberate?? Laimonis Kavalieris -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To:

I have a RedHat 7.0 system. I compiled wine (form CVS, mondya afternoon) with the transgaming patches. For the most part it works pretty good, I can use napster, winamp, and many other programs that are installed on my win98 partition. The problem is when I play games. I donm't have that many games that I want to play (about 5), but I can;t get any of them to work properly. The first game

Hello, In order to do a future forecast based on my past Time Series data sets (salespricesproduct1, salespricesproduct2, etc..), I used arima() functions with different parameter combinations which give the smallest AIC. I also used auto.arima() which finds the parameters with the smallest AICs. But unfortuanetly I could not get satisfactory forecast() results, even sometimes catastrophic

Hi list, thank you for Shorewall :) I''m trying to get a simple config to work but i can''t seem to work out how to gain access via ssh to the protected remote machine. But that doesn''t surprise me really as i have just spend well over an hour to find how to limit the lograte AND fill in the logburst in shorewall.conf. I have specified a logfile (not messages) in

I am running version 3.0.7 of Shorewall on a Debian Sarge system, but when I start Shorewall I get this: /usr/share/shorewall/firewall: line 204: 4: command not found I looked there and found this: # Run ip and if an error occurs, stop the firewall and quit # run_ip() { if ! ip $@ ; then if [ -z "$STOPPING" ]; then error_message "ERROR: Command \"ip

Dear All, I would like to be able to plot histograms/densities on a semi-log or log-log scale. I found several suggestions online http://tolstoy.newcastle.edu.au/R/help/05/09/12044.html https://stat.ethz.ch/pipermail/r-help/2002-June/022295.html http://www.harding.edu/fmccown/R/#histograms Now, consider the code snippet taken from http://www.harding.edu/fmccown/R/#histograms # Get a random

The blowfish crypt(3) mechanism supports the use of a "cost value" for password encryption. The cost value is encoded into the encrypted password that is stored in master.passwd. On OpenBSD, this cost value can be set in login.conf. FreeBSD does not currently support the cost value. The cost value is the base-2 logarithm of the number of rounds of encryption to use so

Is there an easy way to display the SQL in script/console (similar to what is done in development.log)? The following works but it''s ugly: $ script/console Loading development environment. >> class ActiveRecord::ConnectionAdapters::AbstractAdapter >> def logger= (logr) >> @logger = logr >> end >> end => nil >>

Some of you may have seen the pythag() part in the R API definition in "Writing R Extensions" (source = doc/manual/R-exts.texinfo). or followed the report and Prof. Brian Ripley's answer about pythag()'s availability from R's binary. As we say in above manual >> `pythag(A, B)' computes `sqrt(A^2 + B^2)' without overflow or >> destructive

Hi, Is lograte can be used to rotate icecast.log, ... ? IceCast seams to keep files open during is running time. So, when you try to move icecast.log, it still writing is the moved file. Moreover, if you delete the file, icecast doesn't write any more. Best regards, PE --- >8 ---- List archives: http://www.xiph.org/archives/ icecast project homepage: http://www.icecast.org/ To

Hi all, I use XP on my laptop PC, with 30 GB HD (single partition), and a 3.5" hard drive (not a stick), accessed through a USB encasing, partitioned as five logical drives. I have successfully set up syslinux to boot several Ubuntu versions from the primary USB hard drive partition, and no problems there. The problem is, I would like to have an entry in the menu, like "Boot from