similar to: Nonlinear Regression Parameter Shared Across Multiple Data Sets

I'm attempting to calculate a regression in R that I normally use Prism for, because the formula isn't pretty by any means. Prism presents the formula (which is in the Prism equation library as Heterologous competition with depletion, if anyone is curious) in these segments: KdCPM = KdnM*SpAct*Vol*1000 R=NS+1 S=(1+10^(X-LogKi))*KdCPM+Hot a=-1*R b=R*S+NS*Hot+BMax c = -1*Hot*(S*MS+BMax) Y

I have a question regarding robust nonlinear regression with nlrob. I would like to place lower bounds on the parameters, but when I call nlrob with limits it returns the following error: "Error in psi(resid/Scale, ...) : unused argument(s) (lower = list(Asym = 1, mid = 1, scal = 1))" After consulting the documentation I noticed that upper and lower are not listed as parameter in

Hi there, I don't have much experience with fitting at all and I'd like to get some advice how to use the "weights"-argument with nls correctly. I have created some data with a sigmoidal curve shape. Each y-Value was generated by the mean of three values. A standard deviation was calculated too. Now, I'd like to weight the data points respective to its standard

Dear all, I am not sure if this mail is for R-help or should be sent to R-devel instead, and therefore post to both. While using nlrob from package 'robustbase', I ran into the following problem: For psi-functions that can become zero (e.g. psi.bisquare), weights in the internal call to nls can become zero. Example: d <- data.frame(x=1:5,y=c(2,3,5,10,9)) d.nlrob <-

Hi all, I'm trying to fit a nonlinear regression by "nlrob": model3=nlrob(y~a1*x^a2,data=transient,psi=psi.bisquare, start=list(a1=0.02,a2=0.7),maxit=1000) However an error message keeps popping up saying that the function psi.bisquare doesn't exist. I also tried psi.huber, which is supposed to be the default for nlrob: model3=nlrob(y~a1*x^a2,data=transient,psi=psi.huber,

Dear all, ? I implemented a non-linear model using nls(...) and it works just fine. I now tried to run the same model using nlrob(...) which basically does the same but uses a more robust estimation procedure. My problem: I cannot seem to get nlrob(...) running. Irrespective of how I try to call the function, I always get the error message "Error in is.null(x) : 'x' is missing".

Hi friends I just started using wine on my Linux-machine. Mostly because I want to switch finally from Windows to Linux on my computer at work. There is just one single application I want to "take with me" and that is "GraphPad Prism" (http://www.graphpad.com/prism/Prism.htm). So I installed wine0.9.20 on my Mandriva2007. Then I used wine to install Prism in its most recent

What am I failing to understand here? The script below works fine if the dataset being used is DNase1 <- DNase[ DNase$Run == 1, ] per the example given in help(nlrob). Obviously, I am trying to understand how to use nls and nlrob to fit curves to data using R. #package=DAAG attach(codling) plot(pobs~dose) #next command returns 'step factor reduced below min factor

Dear useRs: Is there a way to include a raster image (e.g., .gif, .jpg, .bmp) representing a company logo, a school logo, etc. into an R graphic? For example, it would be nice to be able to include the logo of the school into the charts that a student produces for her dissertation. Similarly, while working on a competitive analysis, it would be nice to be able to replace the names of the

Hello, I am kind of confused with this situation. I am attempting to build a PDC using TLS/SSL with the following version of software. Samba 3.0.23 OpenLDAP 2.3.19 Fedora Core 5 When I startup the Samba server via the "service" command (service smb start) I get the following errors in my logs. Using SSL: Jul 13 09:52:34 prism smbd[23161]: smbldap_search_suffix: Problem during

I'm fitting nonlinear functions to some growth data but I'm getting radically different results in R to another program (Prism). Furthermore the values from the other program give a better fit and seem more realistic. I think there is a problem with the results from the r nls function. The differences only occur with weighted data so I think I'm making a mistake in the weighting.

Hi there . Some one know how fix this??? I tried a lot of diffierent thinks, but nothing. and my Freebsd Access point keep crashing some time when I use cards prism 2.5 , but dosen't crash with prism 2 . The problem is I must to use prims 2.5 because it's a high power card, and work really well until crash. thanks!! -- Marcos Biscaysaqu Systems Administrator ThePacific.Net Ltd.

I have a prism 2.5 firmaware 1.5.6 but I can't make work this in 11Mbps on Freebsd 5.1 but the same card work perfect in Freebsd 5.0 help pls! this is my ifconfig wi0:flags=8947<UP,BROADCAST,DEBUG,RUNNING,PROMISC,SIMPLEX,MULTICAST> mtu 1500 ether 00:60:b3:7a:a3:1d media: IEEE 802.11 Wireless Ethernet DS/11Mbps <hostap> (DS/2Mbps <hostap>) status:

Hi there, I recently installed a plugin for SquirelMail (php imap client, squirrelmail.org) which adds up message sizes per folder to display a nice total size of each subscribed folder. After consulting the plugin maintainers at SquirrelMail in this thread: http://thread.gmane.org/gmane.mail.squirrelmail.plugins/12468/focus=12468 There is quite a discrepancy in the rfc822.size reported by

Is there a way to estimate a nominal response model? To be more specific let's say I want to calibrate: \pi_{v}(\theta_j)=\frac{e^{\xi_{v}+\lambda_{v}\theta_j}}{\sum_{h=1}^m e^{\xi_{h}+\lambda_{h}\theta_j}} Where $\theta_j$ is a the dependent variable and I need to estimate $\xi_{h}$ and $\lambda_{h}$ for $h \in {1...,m}$. Thank you, Mauricio Romero Quantil S.A.S. Cel: 3112231150

I'm attempting to insert variables as expressions into labels for my plots, but I'm running into the issue that characters in my labels aren't allowed in expressions, notably '[' and ']'. An input string of "[^3H]5-CT" needs to be displayed as [3H]5-CT with the 3 as superscript. But brackets aren't particularly liked in expressions. I know I can do

Hi There . Iam runnig a Freebsd5.0 Access Point with a PCMCIA card PRISM 2.5 chipset 200mw , I have got 32 clients on it but is not working fine, keep crashing and I have lot wi0: timesouts on the screen. I was using with same computer other prism card but prism2 and that card work very well but with a lot lower power and more noise because I need use an Ampl. With The prism2.5 the Access

Hi. I would like to add a horizontal line to each column of data on a stripchart (vertical=T) to indicate the mean(rather like the package Prism does - for those that have used this) for each dataset. I presume the best way to do this would be with "lines" but after much trying / playing about I can't figure it out. Can anyone help? Thanks

Dear useRs, Do you know, by any chance, if it is possible to write a string in the header/footer of a pdf file created during a plot session? To make sure that I know what R script generated a certain chart, I would like to add the path to the working directory into the header/footer of the corresponding pdf file. Thank you. Tudor -- Tudor Dan Bodea Georgia Institute of Technology School of

Hi! I have PC running FreeBSD 4.8-STABLE and installed Z-COM XI-626 PCI (Prism 2.5). When I try to directly connect notebook (running WinXP) with acx100 chip based card to FreeBSD I can't connect. I try connect running wi in AP and Ad-Hoc modes. NOP. Cards can see each other but can't connect. :( I is possible to connect these cards? Thanks. Andris