search for: fzd

...oben <- rep(0, each=2) for (i in seq(1:2)) { n <- dim(daten)[1] ix <- sample(n,200) # producing samples samp_i <- daten[ix,] stichproben[i] <- samp_i # doesn???t works # Calculation of the model: posterior_i <- MCMClogit(y ~ fbl.ind + fekq3 + febitda4 + fuvs + fkru + fzd + fur3, data=samp_i, b0=prior, B0=precision, tune=0.5) # calculation ergebnisse[i] <- summary(posterior_i) # saving the results (works) } --snip-- I have a data set called "daten". I produce samples of the size 200. The samples are saved in samp_i. Question: How is the easiest way...

...0) samp_i <- daten[ix,] # draw samples y <- sum(samp_i$y) # number of defaults stichproben200[[i]] <- samp_i # saving the samples default200[i] <- y # saving the number of defaults # Modell berechnen: posterior_i <- MCMClogit(y ~ fbl.ind + fekq3 + febitda4 + fuvs + fkru + fzd + fur3, data=samp_i, b0=prior, B0=precision, tune=0.5) # calculation ergebnisse200[i] <- summary(posterior_i) # saving the results } # write out the solutions into an excel-file write.csv2(ergebnisse200, "ergebnisse.csv") --snip-- My solution has the following form: http://img296.i...

Dear all, I would like to use the Levenberg-Marquardt algorithm for non-linear least-squares regression using function nls.lm. Can anybody help me to find a a way to compute confidence intervals on the fitted parameters as it is possible for nls (using confint.nls, which does not work for nls.lm)? Thank you for your help Michael

...roducible code. ------------------------------ Message: 27 Date: Mon, 22 Nov 2010 09:13:19 -0600 From: Mike Lawrence <Mike.Lawrence at dal.ca> To: "r-help at lists.R-project.org" <r-help at r-project.org> Subject: [R] hierarchical mixture modelling Message-ID: <AANLkTi=fzD+Jz1-DnOY239U3xycCRfriyZ1bF-AJjjAh at mail.gmail.com> Content-Type: text/plain; charset=ISO-8859-1 Hi folks, I have circular data that I'd like to model as a mixture of a uniform and a Von Mises centered on zero. The data are from 2 groups of human participants, where each participant is me...

...zsk|mxSkTE9qcik!WNm)Q{WEo`M=2&3 z^OW<>(S<oQ3RKA~T`sQ)whq1~CH)w5LJzD><wHG9-3J*=nh$0}Wd6R?C^st`+EZuB z6`&W!4=bfUsHLTifKuAHFn(ug6Kj=}e>FziLej*di<H<xDfE_jU#`U7l5(%G<Aj)Z z;>5gwxX@@RI7TWN0|jxkWZ_ at SXo6e)WHCI&=UQ%^v`6dzN=70vO>AyWopy9S_xp_f zT~l1^hI~UI5cEx(piwfzdU+ at lBtmK3AEj4QQ0~Zo8&QQOY7JE=l=@|<mUI-qq at z^+ z4HT)i&9mIA at 6tvT+rzHa+yJQpjgk=^8Q8o!@8>e!zs-0*pYeWy-tW_>IiOB7CXIK6 z9?>Jb*~GfK%T?T4s5hsxZiTFacP0d>=xn{@PKlZngqma!+EG_$#TD~AZkOYtTaLgy zE?e7W|4h!}H<TKLQ|I8s^-LT@^w#Z&vP&cSG#Ej>nVY+4+C-nu62;J|!m~_Lx-sWP z47~T%{d...

dear list, I stumbled over the follwing strange behaviour/error when using `nls' which I'm tempted (despite the implied "dangers") to call a bug: I've written a driver for `nls' which allows specifying the model and the data vectors using arbitrary symbols. these are internally mapped to consistent names, which poses a slight complication when using `deriv' to

using X11(antialias="none") plot(1:10) produces an incomplete plot with 2.7.0. under MacOS 10.4.11: both y-axes are missing.the y-ticks (except the tick at `2' wich is missing) including labels are there. this happens with a 1280x854 labtop monitor. is this known behaviour/a cairo bug? joerg

I tripped over the following problem: due to (my?) unresolved problems with cairo graphics under MacOS (graphic window blocks unrecoverably) I modified my `.Rprofile' to issue `X11(type = "Xlib")` on each R startup (as a workaround to automatic opening of a `cairo' device when issuing a plotting command). I now noted that `R CMD INSTALL' obviously

I observed the following: match.arg("white", colors()) yields 'white', but match.arg("whit", colors()) yields: `Error in match.arg("whit", colors()) : 'arg' should be one of "white", "aliceblue", '... this message actually comes from `pmatch'. using a suitable subset of `colors()' works OK. the precise length

for unweighted fits using `nls' I compute confidence intervals for the fitted model function by using: #------------------- se.fit <- sqrt(apply(rr$m$gradient(), 1, function(x) sum(vcov(rr)*outer(x,x)))) luconf <- yfit + outer(se.fit, qnorm(c(probex, 1 - probex))) #------------------- where `rr' contains an `nls' object, `x' is the independent variable vector, `yfit'

dear list, my question concerns the use of `eval' in defining the model formula for `nls' (version 2.6.2.). consider the following simple example, where the same model and data are used to perform unweighted and weighted fits. I intentionally used very uneven weights to guarantee large differences in the results #================================CUT=========================== ln

I get the error Error : evaluation nested too deeply: infinite recursion / options(expressions=)? during a 'R CMD check ...' on one of my packages. The reason seems to be that this package is mutually dependent on another one (i.e. the DESCRIPTION files of package A lists package B under "Depends" and vice versa). this might be bad design (having bits in both packages needed

this is probably not really a R specific question, if so apologies for off-topic posting: I'm interested in the probability density function of the maximum values from repeated samples of size N from a normal distribution: smp <- rnorm(N, meanval, stdev) with some mean 'meanval' and standard deviation 'stdev'. I would like to know what is the frequency distribution of

2 days ago I asked this on r-help, but no luck... since this is actually a programming question, I post it here again: my question concerns the use of `eval' in defining the model formula for `nls' when performing weighted fits. (I use version 2.6.2., but according to NEWS there were no changes to `nls' in 2.7.0, so the problem is still present). in this scenario their

dear list members, I apologize in advance for posting a second time, but probably after one week chances are, the first try went down the sink.. my question concerns computation of confidence intervals in nonlinear fits with `nls' when weigthing the fit. the seemingly correct procedure does not work as expected, as is detailed in my original post below. any remarks appreciated. greetings

my understanding of `eval' behaviour is obviously incomplete. my question: usually `eval(expr)' and `eval(expr, envir=parent.frame())' should be identical. why does the last `eval' (yielding `r3') in this code _not_ evaluate in the local environment of function `f' but rather in the global environment (if `f' is called from there)?

dear list members, my question concerns computation of confidence intervals in nonlinear fits with `nls' when weigthing the fit. the seemingly correct procedure does not work as (I) expected. I'm posting this here since: (A) the problem might suggest a modification to the `m' component in the return argument of `nls' (making this post formally OK for this list) and (B) I got no

hi, I have two private packages, the first (`pkc') depending on the second one (`roiutils'). The source code and DESCRIPTION files describes the dependency as it should be ('Imports', `require'), at least I think so. now, running R CMD CHECK pkc yields the following output in which I have inserted my questions (lines starting with -->): * checking for working latex ...

I tried two times to get some help for this from the help list but to no avail. I hope, it's OK to post this here (once...). if not, please ignore (the rest of) this mail: for one of my private packages R CMD CHECK {package} givves me annoying warnings due to 'missing links' in the manpages for functions from another private package mentioned in the manpages. specifically, I have two

I initially thought, this should better be posted to r-devel but alas! no response. so I try it here. sory for the lengthy explanation but it seems unavoidable. to quickly see the problem simply copy the litte example below and execute f(n=5) which crashes. called with n != 5 (and of course n>3 since there are 3 parameters in the model...) everything is as it should be. in detail: