similar to: problem with nls()

> From: Siegfried.Macho > > Dear R-helpers, > > 3 questions: > 1. Is there a package that contains a routine for computing > Kendall's W > (coefficient of concordance), with and without ties ? Is that the same as Kendall's tau, as in cor(..., method="kendall")? > 2. Is there a package that contains a routine for computing > Goodman' s

Please do give reproducible example. The code you gave, which you claimed `works correctly' doesn't: > SDT.Optim <- function(crit = NULL, Hess = F) + { + q <- length(par); x <- data + if(hasArg(crit)) + cat("\n Crit present\n") + else +

Dear ladies and gentlemen, I have recently installed the latest version of R (4.3.3) for windows. Now I have the following problems with mutated vowels like ?, ?, etc. Here is an example: If I type the command: Dir <- "C/Users/macho/Documents/_LVn/Experimentelle ?bungen" in the R console there is no problem. However, if I put the same command into a source file (e.g. Test.r) and call

Hola! I am doing some minimal changes on the package SenSrivastava to make it fit for R2.0.0 (This package passed Rcmd check on earlier versions of R, and are passing the CRAN daily checks!) I am using rw2000dev from CRAN. Now i get (running in an Xemacs shell): installing R.css in C:/R/rw2000dev/src/gnuwin32/SenSrivastava.Rcheck ---------- Making package SenSrivastava ------------ adding

Hi there, I am trying to compare nonlinear least squares regression with AIC and anova. The simplest model is one nonlinear curve, and in the more complex model I have a categorical variable (producing parameter estimates for four curves). Both models run fine, but when I try to produce an AIC value for the second model I get the error: > AIC(pow.nls1) [1] 114408.3 > AIC(pow.nls2) Error in

Subject: nls, step factor 0.000488281 reduced below 'minFactor' of 0.000976563 Hi there, I'm trying to conduct nls regression using roughly the below code: nls1 <- nls(y ~ a*(1-exp(-b*x^c)), start=list(a=a1,b=b1,c=c1)) I checked my start values by plotting the relationship etc. but I kept getting an error message saying maximum iterations exceeded. I have tried changing these

Potential bug: I mistyped weights in the call ('weigths') and it did not produce any error= message. The coefs were exactly the same like without weights, so I was su= spicious and when weights(nls1) gave NULL, I saw my typo. Usually the function will say "Unused arguments", which shows you the error= , but not nls. Regards Stephen [[alternative HTML version deleted]]

  Hi, I fitted tree growth data with Chapman-Richards growth function using nls. summary(fit.nls) Formula: Parameters: Estimate Std. Error t value Pr Signif. codes: 0 ''***'' 0.001 ''**'' 0.01 ''*'' 0.05 ''.'' 0.1 '' '' 1 Residual standard error: 1.879 on 713 degrees of freedom Algorithm

Dear R Users, Here is another probelm/question. I would like to run some panel regressions with R. Therefore I have combined several time periods of data for different individuals in my database. I have already run pooled OLS but I would need to calculate a Fixed Effects Estimator (within estimator). Unfortunately I couldn't find anything like that in the RSearch and I suppose that lme

Hi, I am running glm logit regressions with R and I would like to test a linear combination of coefficients (H0: beta1=beta2 against H1: beta1<>beta2). Is there a package for such a test or how can I perform it otherwise (perhaps with logLik() ???)? Additionally I was wondering if there was no routine to calculate pseudo R2s for logit regressions. Currently I am calculating the pseudo R2

Dear list members, I have a question concerning multiple comparisons after using glm. My response variable is days until emergence of an insect species. The explanatory variables are sex (two levels), parasitoids added (two levels) and populations (34 levels). I would like to know now which populations are different in days until insect emergence. For this I used multiple comparisons as

Hello, I'm currently dealing with a rather strange problem when using the function "stepAIC" ("MASS" package). The setting is the following: From model learning data sets ("learndata"), I want to be able to build prediction functions (in order to save them in a file for further use). This is done by the function "pred.function" (see below). Therein,

Hi, Sorry, I accidentally rejected this message instead of approving it. Geoff. ---------- Forwarded message ---------- Date: Fri, 11 Mar 2005 17:33:06 +0100 From: Siegfried Wagner <esiggi@gmail.com> Reply-To: siegfried@esiggi.net To: icecast-dev@xiph.org Subject: Re: [Icecast-dev] patch for icecast-2.2.0 to add client maxtime Hi, I did your recommended changes on my patch. But with the

Hello, I had downloaded the samba 2.2.7a tarball, the pubkey and the samba-2.2.7a.tar.asc files from your website. The import of the samba-pubkey.asc works probably, but by the verify I get the errormessage that no path to this signature is defined and the signature can't verified. The distribution of Linux is the SuSE 8.0 and I used this lines. gpg --import samba-pubkey.asc

Hi, sorry I forgot to mention: > Siegfried, > This sounds good, and I will be trying it. I just wondered could this > timer be applied per connection or is it global per mount point ( > sorry havent had chance to look at the diff yet ). I tried to add the code so that it also could be defined per mount point but I didn't test this out. But of course this makes much more

Sigfried, I am not a statistician, but I have learned that according to the Central Limit Theorem (CLT) sums of random variables, regardless of their form, will tend to be normally distributed. CLT does not require the variables in the sum to come from the same underlying distribution. Ciao, Hannu Kahra Progetti Speciali Monte Paschi Asset Management SGR S.p.A. Via San Vittore, 37 IT-20123

Hi, David H.H. announced recently that jquery is going to be the default in Rails 3.1, and that Prototype helpers / RJS are going into a gem. What does that mean for the future? Should we progressively forget about things like javascript helpers, RJS, and all these fun parts of rails? And start coding with jquery/json in mind? I''m not against that at all, but I have to admit I find RJS

Dear Siegfried, I believe your boss is wrong saying that: >He also tried to explain me that the monthly means >(based on the daily measurements) must follow a >log-normal distribution too then over the course of a year. every statistician know that increasing the sample size the sample distribution of the mean is proxy to a gaussian distribution (Central Limit Theorem) independently

I'd like to use linear and quadratic trend analysis in order to find out a change in slope. Basically, I need to solve a similar problem as discussed in http://www.gseis.ucla.edu/courses/ed230bc1/cnotes4/trend1.html My subjects have counted dots: one dot, two dots, etc. up to 9 dots. The reaction time increases with increasing dots. The theory is that 1 up to 3 or 4 points can be counted

Hi all, I thought that readers of R-Help might find the following article at ScienceNews of interest: Odds Are, It's Wrong Science fails to face the shortcomings of statistics By Tom Siegfried March 27th, 2010; Vol.177 #7 (p. 26) http://www.sciencenews.org/view/feature/id/57091/title/Odds_are,_its_wrong Regards, Marc Schwartz