similar to: fit a gamma pdf using Residual Sum-of-Squares

Dear all, in one part of my code I want to sum two vectors element-wise the problem is that either the 1st vector or the 2nd vector always have one or two less elements example of my problem In TotalVector + (datalist2[[1]]$dataset$Results[[j]]$Results[[time]]$Sweep) :   Länge des längeren Objektes          ist kein Vielfaches der Länge des kürzeren Objektes Browse[1]> str(TotalVector)  int

Hi, In R, I know the method table(), which builds a contingency table of the counts y at each level for the factor x. But I would like to know what is the inverse method of table(), if it exists, to obtain the vector z, from the two vectors x and y? x <- (86, 90, 94, 98, 102, 106, 110, 114, 118, 122, 126, 130, 134, 138, 142, 146, 150, 154, 158, 162, 166, 170, 174) y <- c(2, 5, 10, 17,

I am looking for a way to get a warning message immediately after an evaluation within a function. To get error messages you can use geterrmessage(). But I found no function that allows me to check for warnings. Five years ago this questions has been posted but I haven't found any answer. Thanks for any help. Peter Wolf ------------ For illustration purpose a simple example follows: you

R 2.7.2 PPC Mac OS X 10.4.11 library mice 1.13.1 I try to use mice for multivariate data imputation. My variables are numeric, factors, count data, ordered factors. First I created a vector for the methods to use with each variable ImpMethMice<-c(rep("logreg", 62), rep("polyreg",1), rep("norm",12), rep("polyreg",12)) next step was

Dear list I tried to read an xml file using the xml package. Unfortunately, some encoding problems occure. E.g. german Umlaut will be red correctly. I assume that the occurs due to (internal?) conversion to utf-8. To illustrate the problem, I have wrote to xml files. File Test 1 ----------- <?xml version="1.0" encoding="ISO-8859-1"?> <Daten> <ITEM>

Hello everybody, I have a problem with the commando of boxplot -bclust. http://127.0.0.1:13155/library/e1071/html/boxplot.bclust.html > data(iris) > bc1 <- bclust(iris[,1:4], 3, base.centers=5) Committee Member: 1(1) 2(1) 3(1) 4(1) 5(1) 6(1) 7(1) 8(1) 9(1) 10(1) Computing Hierarchical Clustering > boxplot(bc1) Warnmeldungen: 1: In if (x$datamean) { : Bedingung hat Länge > 1 und

There are two misspellings in the german Error message for friedman test: Fehler in friedman.test.default(cont$score, group = cont$goup, blocks = cont$cont) : y, Gruppen und blöcke müssen die sekbe Länge haben The correct spelling would be: Fehler in friedman.test.default(cont$score, group = cont$goup, blocks = cont$cont) : y, Gruppen und *Blöcke* müssen die *selbe* Länge haben Instead of

Hi all, I am trying to use in RStudio the latest code given in https://github.com/systematicinvestor/SIT/blob/master/R/bt.test.r, which seems to work fine but with the following warning for download limits (one for each of the tickers). I searched in options() something which could be related to this setting, w/o success. Any hint for me in order to raise or remove these limits? Where is this

Folks; I have looked for this subject extensively using all the resources I have available and could not find a solution; therefore I am sending you this email seeking for help. I need to solve a discontinuous nonlinear function, but R 2.6 keeps giving me an error. Below is a simple example: f <- function (t,k,l,u) { r0<-346.1 x<-t-l if (x<0) r0 else (r0-u)*exp(-k*x)+u

Greetings R users, maybe there is someone who can help me with this problem: I define a function "optim.fun" and want as output the sum of squared errors between predicted and measured values, as follows: optim.fun <- function (ST04, SM08b, ch2no, a, b, d, E) { predR <- (a*SM08b^I(2)+b*SM08b+d)*exp(E*((1/(283.15-227.13))-(1/(ST04+273.15-227.13)))) abserr

Dear folks, We are running SAMBA (2.0.5.a) in a test-enviroment using an IBM Risc 6000 / 370 (AIX 4.3.2) as server. The frontends are pc's running NT 4.0 (SP6). Is it possible to connect from one frontend to the same share as two different users (with different passwords and privileges)? The connection works out perfectly for every single user but when we try to connect one user after the

Dear all, I am new R user and I am sure that this question has been asked quite often and I have also googled it and read about it! I understood that in order to read excel sheet into R you need to open it and saved it as csv or text, is this true? or you can use read.delim2 and read.csv2 to do this without the following error > dat <- read.csv2(file="C:\\Dokumente und

Hello together, this is ma first post, so please aplogize me if post this in the wrong section. I have problem concerning ma two matrizes. After a regressione and so on, I got two matrizes Matrixres contains the results of ma calculation. Matrixr contains my detiene, which where Aldo used for the regression. Please ser the following code: #Datei einlesen residual =

Hi, I want to calculate the t- and p-values for a linear model using the Newey West estimator. I tried this Code and it usually worked just fine: > oberlm <- lm(DYH ~ BIP + Infl + EOil, data=HU_H) > coeftest(oberlm, NeweyWest(oberlm, lag=2)) t test of coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 0.1509950 0.0743832 2.0300 0.179486 BIP

Hi Everyone, I'm trying to write selfStart non-linear models for use with nls. In these models some combinations of parameter values are illegal; the function value is undefined. That's OK when calling the function directly [e.g. SSmodel(x, pars...)]; I return an appropriate non-value such as NA or Inf. However, when called from nls [e.g. nls(y~SSmodel(x, pars...), ...)] those

Hi all, I am not a mathematician and I am trying to fit a function which could fit my observed data. Which function should I use and how could I fit it to data in R? Below are the data: x <- c(0, 9, 17, 24, 28, 30) y <- c(500, 480, 420, 300, 160, 5) I use R for Mac OS, version 2.10-1 2009-08-24 Thank you for your help. Vincent. [[alternative HTML version deleted]]

Hi, I have some data I want to fit with a non-linear function using nls, but it won't solve. > regresjon<-nls(lcfu~lN0+log10(1-(1-10^(k*t))^m), data=cfu_data, > start=(list(lN0 = 7.6, k = -0.08, m = 2))) Error in nls(lcfu ~ lN0 + log10(1 - (1 - 10^(k * t))^m), data = cfu_data, : step factor 0.000488281 reduced below 'minFactor' of 0.000976562 Tried to increase minFactor

Hello, I am in fly-by mode as I am still more than busy, so I just dumb this one out for anyone interested. Adapt it for your needs if you like it. This effect is similar to ScrollTo, but scrolls the viewport vertically the minimal amount necessary to bring the element into view. Usage: new Effect.ScrollToFullView(element, {duration: 0.6, offset: 10}); Offset is the distance in pixel to the

Hi Everyone, I posted this a couple of weeks ago with no responses. My interface (via gmane) seemed a bit flakey at the time, so I'm venturing to repost with some additional information. I'm trying to write selfStart non-linear models for use with nls. In these models some combinations of parameter values are illegal; the function value is undefined. That's OK when calling the

Dear all, I have a huge problem after downloading and exporting data from Reuters3000 XTra: I downloaded many many monthly, quarterly and yearly data. I do not know why, but after saving, I have mixed-data sets, i.e. files which can be imported as “read.csv” and others that are in the format of “read.csv2”. Sure I could change them, but normally it should be possible to mix them… For