similar to: nonlinear fitting when both x and y having measurement e

Dear r-help, I want to conduct nonlinear fitting to a data frame having x and y variables. Because both x and y have measurement error, I want to include error term of x variable in the model. I'm not sure but I think ordinary nls model only consider error term of y variable. How can I do this kind of nonlinear fitting in R. Is there any examples in nls package? Thanks in advance, --

Hi, I recently tried to conduct non-parametric multiple comparison with npmc package, and encountered some problem; sometimes it stopped with following error. Error in uniroot(f = function(arg) p - z.dist(arg, corr = corr, df = df, : f() values at end points not of opposite sign Although I'm not quite sure about calculation of z.quaintile in npmc.R, it seems increment of

I would like to install the package "growth" as it contains the function "corgram" and some other presumably useful stuff for time series analysis. I can see it is in R standard library list: http://hosho.ees.hokudai.ac.jp/~kubo/Rdoc/doc/html/packages.html<http://hosho.ees.hokudai.ac.jp/%7Ekubo/Rdoc/doc/html/packages.html>

I am using this code to calculate the moving average mean.It worked fine but when I wanted to also calculate based on sd(stander deviation) I got the error shown below. I read this documentation of R movingFun and found that sd was mentioned at http://hosho.ees.hokudai.ac.jp/~kubo/Rdoc/library/raster/html/movingFun.html I wonder if needs any other things to conseder or sd is not at all supplied

Hi R users, I'm a quite extensive user of the heatmap() function, and as many others i'm frequently frustrated by its fixed square width / height ratio. Here is a typical example : dev.new() heatmap(matrix(rnorm(1200), nrow=10)) I have a non square matrix to plot with heatmap, as there are many more columns than rows, row labels are larger than necessary but column labels

Hi Can anyone tell me how to convert an image file (a jpg file for instance) into a matrix file? The command image(blabla) transforms a matrix into an image file, I am searching exactly the inverse. Thanks... [[alternative HTML version deleted]]

Dear researchers I wish to plot a box plot without the mean line (the black line) and the i wish a full line for the standard deviation This is an example mytest <- c(2.1,2.6,2.7,3.2,4.1,4.3,5.2,5.1,4.8,1.8,1.4,2.5,2.7,3.1,2.6,2.8) boxplot(mytest) really thanks Gianni [[alternative HTML version deleted]]

Dear editor: I'm writing this letter to enquire where can I download the package of "lmtest". Can you send me this package? THanks a lot. Best regards, Shu-Fei Wu

Could someone please help me by giving me a reference to how one computes standard errors for the coefficients in an orthogonal linear regression, or perhaps someone has some R code? (I would accept a derivation or formula, but as a former teacher, I know how that can rankle.) I tried to imitate what's done in the code for lm() but went astray somewhere and got nonsense. (This type of

I just switched to Mongrel, and it''s been working much better than my previous lighttpd/fastcgi setup. So thanks for the awesomeness. My current problem: once or twice an hour, I get following error in production Mongrel timed out this thread: too many open files I never get it in testing or on our staging server. Any ideas what would cause that? It doesn''t *appear*

Hello, I am trying to save my work space but I recieve the following message error in save.image(name), how can I overcome this problem? Ahmad A. Al-Mallahi Laboratory of Bioproduction Engineering Graduate School of Agriculture Hokkaido University Sapporo-Japan

I am experiencing similar issues as noted in this bug id: http://bugzilla.mindrot.org/show_bug.cgi?id=510 I am ssh'ing from a dchp'd address to a nat'd address (tried both hostname & ip). after a successful login, I launch an X app. Shortly thereafter I get: "Disconnecting: Corrupted MAC on input." I was not experiencing this problem w/3.7, but I cannot place full

Dear R helpers, For an allometric analysis (allometric equation y = a*x^b) I would like to apply a non-linear regression instead of using log-log transformations of the measured parameters x and y and a Model II linear regression. Since both of the variables x and y are random, I would like to apply a Model II non-linear analog of either Reduced Major Axis or Major Axis Regression. The

Hi, I'm using (or trying to at least) a combination of PAM/passwd and MySQL authentication: auth vmail { mechanisms = plain userdb = sql /etc/dovecot/dovecot-sql.conf passdb = sql /etc/dovecot/dovecot-sql.conf user = root } auth system { mechanisms = plain userdb = passwd passdb = pam imap user = root } MySQL

Hi all, As part of my research I want to investigate the relation between the grid's geometry and the memory accesses of a kernel in common gpu benchmarks (e.g Rodinia, Polybench etc). As a first step i want to answer the following question: - Given a kernel function with M possible memory accesses. For how many of those M accesses we can statically infer its location given concrete values

@Ees, Oh, I see what you mean now. Doing such analysis would be useful for a thread block and not just a single thread but as you say you are onto something bigger than just a thread. We had published a short paper in ICS around this which uses polyhedral techniques to do such analysis and reason about uncoalesced access patterns in Cuda programs. You can find paper at

I hope somebody can help. A student of mine is doing a study on Measurement Error models (errors-in-variables, total least squares, etc.). I have an old reference to a "multi archive" that contains leiv3: Programs for best line fitting with errors in both coordinates. (The date is October 1989, by B.D. Ripley et al.) I have done a search for something similar in R withour success. Has

Dear R users,    my question concerns my interest in comparing the beta coefficients between two identical regression models in two (actually 3) groups. Disclaimer: I am quite new to R, so I might be missing some terminology that I have not come across.   I am trying to find out if I can easily implement a Wald test in R for this, but the only relevant thing that I came across is this link

I have on ubuntu a digital library application Kete which gives in the log/production.log: Please help, thnx. /!\ FAILSAFE /!\ Fri Feb 05 19:05:08 +0200 2010 Status: 500 Internal Server Error No connection to server (localhost:11211 DEAD (Errno::ECONNREFUSED: Connection refused - connect(2)), will retry at Fri Feb 05 19:05:38 +0200 2010)

Michael makes a great point about aliasing here and different indexing that accesses the same element! Another note: x = A[0][2] is fundamentally different depending on the type of `A`. If e.g. A was declared: int A[10][20], there's only _one_ load. A is a (and is treated as) a linear buffer, and GEPs only pinpoint the specific position of A[0][2] in this buffer (i.e. 0*10 + 2). But if A was