similar to: Plot observed vs. fitted values (weighted nls)

Dear Wiza[R]ds, # I have the following experimentally observed data: csdata <- data.frame( time=c(0,1,3,9,20), conc=c(638.697,395.69,199.00,141.58,112.16) ) # weighting resp means response wt.MM<- function(resp, time,A1,a1,A2,a2) { pred <- A1*exp(-a1*time)+A2*exp(-a2*time) (resp - pred) / sqrt(pred) } # Fit using nls cs.wt <- nls( ~ wt.MM(conc, time,A1,a1,A2,a2),

Hi, I have a scatter plot of the variables GNI and Lifeexp (Gross National Income and Life Expectancy, both metric). So I plotted them and I want to add a regression line and a lowess line. I use lowess and not loess because I have missing values. My code: plot(GNI,Lifeexp) abline(lm(Lifeexp~GNI), col="red") y.loess<-loess(Lifeexp~GNI,na. action = na.exclude)

I have installed Photoshop 7 via winetools. However, I cannot start it: dotancohen@ubuntu:~$ photoshop detecting Wine version... done. Drive C: is /home/dotancohen/.wine/drive_c Wine 0.9.7 wine is executed as wine Parameters are dotancohen@ubuntu:~$ cd ~/.wine/drive_c/Program\ Files/Adobe/Photoshop\ 7.0/ dotancohen@ubuntu:~/.wine/drive_c/Program Files/Adobe/Photoshop 7.0$ ./Photoshop.exe

sorry for cross-posting Dear all, I have tow (several) bivariate distributions with a known mean and variance-covariance structure (hence a known density function) that I would like to compare in order to get an intersect that tells me something about "how different" these distributions are (as t-statistics for univariate distributions). In order to visualize what I mean hear a little

hi. i''ve got a couple migrations where i''ve added foreign keys and tried to update them via the migration but my items don''t save. i''ve used save! in the hopes that an exception gets raised but no luck. the migrations run without errors but the data isn''t in the db. can anyone point out what i might be doing wrong? thanks! code below. class

Hi, I'm testing affyPara on Debian lenny with R 2.7 from unstable, and the corresponding bioconductor packages from bioconductor.org downloaded using biocLite as per usual. The command cl = makeMPIcluster(k) succeeds, but spawns k R slaves which promptly soak up all available CPU, which is odd since they are not doing anything yet. This looks like a bug to me. Similar behaviour is shown

I have been attempting to estimate the periodic contribution of an effect to some data but have not been able to fit a sine wave within R. It would be nice to start by being able to fit a sine wave with an amplitude and frequency. x<-seq(0,20,by=0.5) y<-2*sin(2*pi*.5*x) #amplitude =2, frequency=0.5 # This failed to converge r<-nls(y ~ A*sin(2*pi*F*x), start=list(A = 1, F = 1),

On Wed, 2005-12-21 at 15:00 +0100, edd at debian.org wrote: > --CE+1k2dSO48ffgeK > Content-Type: text/plain; charset=us-ascii > Content-Disposition: inline > > Resend after type in mailheader. I have the bug on my system with > yesterday's R 2.2.1. A library mismatch is still a likely cause. > > Dirk <SNIP> > On 21 December 2005 at 01:50, Ethan

--CE+1k2dSO48ffgeK Content-Type: text/plain; charset=us-ascii Content-Disposition: inline Resend after type in mailheader. I have the bug on my system with yesterday's R 2.2.1. A library mismatch is still a likely cause. Dirk -- Hell, there are no rules here - we're trying to accomplish something. -- Thomas A. Edison

Hello i was wondering how i can change the IP address information for my Asterisk box, IP addy, Gateway, DNS. I have a smoothwall router that i am using and i am tring to put the Asterisk box on the orange interface so if anyone can help me please i can use it. Thanks alot William Ray -------------- next part -------------- An HTML attachment was scrubbed... URL:

On Wed, 21 Dec 2005 edd at debian.org wrote: > > --CE+1k2dSO48ffgeK > Content-Type: text/plain; charset=us-ascii > Content-Disposition: inline > > Resend after type in mailheader. I have the bug on my system with > yesterday's R 2.2.1. A library mismatch is still a likely cause. This seems to be a readline 5.1 error. I was unable to reproduce it with 4.3 or 5.0. gdb)

Hi, stepAIC generic plot function creates useful graphics for the diagnosis of multiple regressions. To create predicted versus observed plots, I use to look for the coefficients, copy them by hand, calculate R?, then plot. Is there a more automated way to plot predicted versus observed with its associated R? output using stepAIC, or another function? Kind regards, S.-?. Parent Universit?

Dear useR's, I was comparing results for a logistic regression model between different library's. themodel formula is arranged as follows: response ~ (intercept) + value + group OR: glm( response ~ (intercept) + value + group , family=binomial(link='logit')) lrm( response ~ (intercept) + value + group ) ROC( from = response ~ (intercept) + value + group ,

Dear , This is Elaine. I am computing moran's I using moran.test for a generalized linear model (multiregression). The following contents are the results, and I cannot find the observed Moran's I mentioned as estimate in the manual. Please kindly help indicate if there is observed Moran's I did not notice or other method for calculation. Thanks Elaine data: residuals(modg1)

The following code library(Hmisc) x = 1:10 y = x latex(summary(x~y),table.env=FALSE) latex(summary(cbind(x,y)),table.env=FALSE) should produce latex output that is not a table. The second one produces just a tabular, as it should. However, the first one produces a tabular embedded in a table. (This is the effect if you leave table.env=FALSE out).

Hi, I've two fitted models, one binomial model with presence-absence data that predicts probability of presence and one gaussian model (normal or log-normal abundances). I would like to evaluate these models not on their capability of adjustment but on their capability of prediction by calculating the (log)likelihood between predicted and observed values for each type of model. I found

Hi All, Not sure exactly a memory leak or not. I was porting my nagios from Redhat 7.3 to CentOS 5.4 and I observed the memory usage was gradually increasing on the new centos box. When I ran all my perl plugins with Valgrind -3.2.1, all the plugins complained about a memory leak. Not sure if it's a leak or if its fault with valgrind way of evaluating things. I have attached

(If I am so lucky, please include me in reply, I'm not on the list) Tony Breeds wrote (2004-10-13): > On Wed, Oct 13, 2004 at 04:10:23PM +1300, Michael Woodhams wrote: >> Background: Linux, Debian (Sarge). I want to auto-mount an smbfs at >> boot. smbmount version is 3.0.7-Debian. > > <snip> > >> username=<user>/<domain> >>

Dear friends. I believe this problem has been discussed in various forms now and then, so I hope you will forgive me I ask how to do a truncated model like this, where the observed y is recorded as 10 whenever it is higher or equal to that value x <- rnorm(1000) y <- 10*x + rnorm(1000) y[which(y>10)] <- 10 and recover the "true" model ? Best wishes Troels Troels Ring,

Hello R Users, I have two data sets i) TOMS aerosol optical depth(AOD) and ii) TOMS ozone(O3). AOD data is on 1x1 grid and O3 data is on 5x5 grid. First I want to read AOD and O3 as it is and then I want to regrid AOD on 5x5 grid as O3. Reading is first problem. FIRST PROBLEM READING AOD: AOD data is in following format: ######### Latitute: 89.5 167 0 0 0 0 0 182 0 0 0 0 0 0 0 0 0 0 0 0 200