similar to: Network from package functions

Dear R-helpers, I´m in the context of writing a general function for error propagation in R. There are somehow a few questions I would like to ask (discuss), as my statistical knowledge is somewhat restricted. Below is the function I wrote, the questions are marked. Many thanks in advance. propagate <- function(expr, varList, type = c("stat", "raw"), cov = TRUE) {

Dear all, I encounter some discrepancies when comparing the deviance of a weighted and unweigthed model with the AIC values. A general example (from 'nls'): DNase1 <- subset(DNase, Run == 1) fm1DNase1 <- nls(density ~ SSlogis(log(conc), Asym, xmid, scal), DNase1) This is the unweighted fit, in the code of 'nls' one can see that 'nls' generates a vector

Dear All, I encounter some discrepancies when comparing the deviance of a weighted and unweigthed model with the AIC values. A general example (from 'nls'): DNase1 <- subset(DNase, Run == 1) fm1DNase1 <- nls(density ~ SSlogis(log(conc), Asym, xmid, scal), DNase1) Now for a weighted fit: fm2DNase1 <- nls(density ~ SSlogis(log(conc), Asym, xmid, scal),

I am using R on a Windows XP professional platform. The following code is part of a bigger one CODE press=function(y,x){ library(qpcR) models.press=numeric(0) cat("\n") dep=y print(dep) indep=log(x) print(indep) yfit=dep-PRESS(lm(dep~indep))[[2]] cat("\n yfit\n") print(yfit) yfit.orig=yfit presid=y-yfit.orig press=sum(presid^2)

Hi all, I am wondering what people are doing for security when registering IP phone's remotely if you do not have the equipment to do a VPN tunnel at the remote site. The phone I would be working with mainly is the Polycom lineup. Thanks, Connor Spiess -------------- next part -------------- An HTML attachment was scrubbed... URL:

Ich werde au?er Haus sein von 18-06-99 Bis 28-06-99. Ich werde Ihre Nachrichten nach meiner R?ckkehr beantworten. Bitte wenden Sie sich in dringenden F?llen an: Frank R?seler, Tel: MON 3855 Uwe Spiess, Tel: MON 3743 Guido Zimmermann, Tel: MON 3962 Ferdinand Gro?, Tel: MON 2368 lieber Gru? Axel Hellmig

Hello Bernardo, --------- If I understood your problem this script solve your problem: q<-0.15 + c(-.1,0,.1) h<-10 + c(-.1,0,.1) 5*q*h [1] 2.475 7.500 12.625 --------- OK, this solves the simple example. But what if the example is not that simple. E.g. P = 5 * q/h Here, to get the maximum tolerances for P, we need to divide the maximum value for q by the minimum value for h, and

Hi, I got this error every time when I put upper and lower bound on my parameters in nls. I'm not quite sure what this error is about. I also wanted to find the content of nls_port_fit. But typing this directly into R doesn't seem to give any result out. I also tried to see the content of If anybody could shed some light I'd greatly appreciate. Thank you. rc

Full_Name: Manuel A. Morales Version: 2.8.1 OS: Linux Submission from: (NULL) (137.165.199.246) When fitting a model in nls using the algorithm port with constraints and the shorthand parameter[factor] in the model, I get the following error message: "Error in nls_port_fit(m, start, lower, upper, control, trace) : (list) object cannot be coerced to type 'double' In addition:

Dear All, I was trying to call a self-defined function that performs nonlinear regression and gets the corresponding prediction upper limit using nls2 package. However, weird thing happened. When I called the function in the main program, an error message "fitted(nlsmodel): object 'nlsmodel' not found" came up. But when I directly ran the codes inside the function, no error came

I apologize for posting this in essence the second time (no light at the end of the tunnel yet..): is there a way to enforce that "nls" takes both, the data *and* the model definition from the parent environment? the following fragment shows the problem. #======== cut here========== wrapper <- function (choose=0) { x <- seq(0,2*pi,len=100) y <- sin(1.5*x); y <-

Dear R useRs, I have a vector with positive and negative numbers: A=c(0,1,2,3,0,4,5) Now if i-th element in vector A is > 0, then i-th element in vector B is a+1 else i-th element in vector b=a (or 0) vector A: 0 1 2 3 0 4 5 vector B: 0 2 3 4 0 5 6 What's the right way to do this. I still have some problems with for and if statements... Cheers, Andrej

Dear useRs, I have two different length vectors: one column (1...m) and one row vector (1...n): 20 40 20 60 5 4 2 Now I have to calculate ratios between column vector elements and each row vector elements: 4 5 10 8 10 20 4 5 20 15 12 30 Thank's in advance for any suggestions, Andrej

Hello I was trying to estimate the weibull model using nls after putting OLS values as the initial inputs to NLS. I tried multiple times but still i m getting the same error of Error in nlsModel(formula, mf, start, wts) : singular gradient matrix at initial parameter estimates. The Program is as below > vel <- c(1,2,3,4,5,6,7,8,9,10,11,12,13,14) > df <- data.frame(conc, vel) >

Dear listers, I have two questions: (1) Is there a way in R to change the base-n of the calculations. I wnat to run some calculations either in binary (base-2) or base-4. Is there a way to specify that in R - to chnage from the decimal? (2) I also want to extract the digits from a larger number and store them as a vector. I could do it through converting top string, parsing the string and

Hai Andrej, Oeps, Sorry, my mistake 4.6.x needs at least 1.1.29. But you can upgrade to 4.6.x ;-) Check again and if that problem then still exists, then i think you can better ask this one in the samba-technical list. Or if Rowland knows something more about this. One tip, wait untill tomorrow or the weekend, there wil be a new release due to the regression bug. There planning is to be

hello, I try to use samba as PDC on NetBSD, but I 've had problem with that combination.. samba work perfectly..but if someone want to write into homedir, they couldn't, and w2k workstation appear message , that disk is full.. disk isn't full, and permission are set correct I don't know, where should be a mistake, I check configuration problably 100x, but I haven't found any

HI all, I might add some more information in order to possibly solve my problem. I'm really stuck and no obvious solutions do the trick. I'm using R 1.7.1 on Windows 2000 with the packages regurarly updated. I'm using hypothetical data constructed as a pseudo population conforming to a certain Var-Cov structure. I might add that just > predict(level2) works. But when I add the

Hello Rowland, of course it will be started by samba, I saw this output if I run "samba -i". But I can trigger this output also by starting samba_kcc manually. Andrej -----Ursprüngliche Nachricht----- Von: Rowland Penny [mailto:rpenny at samba.org] Gesendet: Donnerstag, 9. November 2017 14:04 An: samba at lists.samba.org Cc: Andrej Gessel <Andrej.Gessel at janztec.com>

Dear R useRs, I have big (23000 rows), vertical bar delimited file: e.g. A00001|Text a,Text b, Text c|345 A00002|Text bla|456 ... .. . Try using A <- read.table('filename.txt', header=FALSE,sep='\|') process stop at line 11975 with warning message: number of items read is not a multiple of the number of columns I have no problems with processing similar file, which is