similar to: NEWTON ALGORITHM

Hello, I'm trying to fit some points with a 8-degrees polynom (result of lm is stored in pfit). In most of the case, it is ok but for some others, some coefficients are "NA". I don't really understand the meaning of these "NA". And the problem is that I can't perform a derivation (pderiv<-as.function((deriv(polynomial(pfit$coefficients))))) on pfit due to the

Good morning to all, I am encountering a blocking issue when using the function 'breackpoints' from package 'strucchange'. *Context:* I use a data frame, 248 observations of 5 variables, no NA. I compute a linear model, as y~x1+...+x4 x4 is a dummy variable (0 or 1). I want to check this model for structural changes. *Process & issues:* *First, I used function Fstats.* It

Hi, I want to assign the ar1 , ma 1 and the intercept estimated by the following<http://r.789695.n4.nabble.com/assign-estimated-values-td4658138.html#>code to three variables a, b and c respectively. Can anyone help me with this please? code: a0 = 0.05; a1 = 0.1; b1 = 0.85 nu = rnorm(2500) epsi = rep(0, 2500) h = rep(0, 2500) for (i in 2: 2500) { h[i] = a0 + a1 * epsi[i-1]^2 + b1 * h[i-1]

Hi, I'm still not quite there with my H-F (G-G) correction code. I have it working for the main effects, but I just can't figure out how to do it for the effect interactions. The thing I really don't know (and can't find anything about) is how to calculate the covariance matrix for the interaction between the two (or even n) main factors. I've looked through some books

Hi, I'm still not quite there with my H-F (G-G) correction code. I have it working for the main effects, but I just can't figure out how to do it for the effect interactions. The thing I really don't know (and can't find anything about) is how to calculate the covariance matrix for the interaction between the two (or even n) main factors. I've looked through some books

I'm an experimental psychologist and when I run ANOVA analysis in SPSS, I normally ask for a test of non-sphericity (Box's M-test). I also ask for output of the corrections for non-sphericity, such as Greenhouse-Geisser and Huhn-Feldt. These tests and correction factors are commonly used in the journals for experimental and other psychology reports. I have been switching from SPSS to R

Hi, I'm trying to make a regression of the form : formula <- y ~ Asym_inf + Asym_sup * ( (1 / (1 + (n1 * (exp( (tmid1-x) / scal1) )^(1/n1) ) ) ) - (1 / (1 + (n2 * (exp( (tmid2-x) / scal2) )^(1/n2) ) ) ) ) which is a sum of the generalized logistic model proposed by richards. with data such as these: x <- c(88,113,128,143,157,172,184,198,210,226,240,249,263,284,302,340) y <-

Hi everybody, I try to do an exponential model with decreasing. In my data, there is no data missing but there negative values. I adapted this following code whose the origin code comes from in this forum: Regress<-read.table("C:\\Users\\Regression.csv",sep=";",dec=",",header=TRUE) Regress f <- function(x,a,b) {a * exp(b * x)} x<-Regress$EWT

Hi, I'm using R to treat a table (with a lot of missing values) with Rulefit. The matter is when I use the command rfmod. Actually, I don't know how to deal with the error message. I don't know were "true" or "false" is missing. Someone can help me? Thanks The following part is the script I used with the error at the end. platform = "windows" rfhome

> > On Sat, 25 Apr 1998, Helgi Tomasson wrote: > > > I am creating postscript-graphics in R. When I import them into > > > a LATEX document it is impossible to get the caption close > > > to the figure. I run the same program in Splus, and no problem > > > of getting the caption just below the figure. Another thing, > > > the ps file that S

I am trying to get a SPA-3000 to work behind NAT - for the sake of the exercice. The SPA is on the local network at the address 192.168.0.125 behind a NATted linux router. The machine I am trying to work with is a friend's (let's call it lolo.dyndns.org) and I've installed Asterisk 1.0.3 on it. I can see the SPA register but when I try to make an outbound call I get the message:

Hello evryone, Thank you for helping me to choose a CRAN MIRROR. When Icome to choose a mirror an error message comes: >Message d'avis : In open.connection(con, "r") : connexion ? 'cran.r-project.org' impossible sur le port 80. Again, I managed to download th package ttda.zip but an error message during the installation comes to me: >Erreur dans gzfile(file,

I am creating postscript-graphics in R. When I import them into a LATEX document it is impossible to get the caption close to the figure. I run the same program in Splus, and no problem of getting the caption just below the figure. Another thing, the ps file that S created is only 10% of the size of the file that R created. What is wrong? Are there some parameters to be set? Helgi

Hi, the second step in my exercice is to calculate the sum of the amout for each class et not the frequency i have this vector x           y 1          100 2          1500 3          3250 4          6250 5          2000 6           450 i want to use the function table and cat to calculate the sum of the amount in each class [0-1000], [1000-3000],[> 3000] Thank you for your help??

Hi I would like to change the locale when using R CMD BATCH. Usually, if I want to run it in english, for R in console/GUIs, I edit the .Rprofile file, adding: Sys.setlocale("LC_ALL","en_US.UTF8") Sys.setlocale("LC_MESSAGES","en_US.UTF8") But while this works for interactive R, it does not for R CMD BATCH. The problem is that running tests for a package,

hello! je veux importer un fichier de donnees excel que j'ai au prealable converti en fichier txt avec separateurs tab, fichier de la forme entree simple (suite de colonnes contenant des variables). Voila ma ligne de commande : > poussins <- read.table("poussins.txt", header=T, sep="\t") et sa reponse Erreur dans scan(file = file, what = what, sep = sep, quote =

Dear r-help, Could you help me to find a solution for this error: Il y a eu 50 avis ou plus (utilisez warnings() pour voir les 50 premiers) > warnings() Messages d'avis : 1: In if ((data[pa, k] == df[, j]) & (data[ch, k] == i)) { ... : la condition a une longueur > 1 et seul le premier élément est utilisé 2: In if ((data[pa, k] == df[, j]) & (data[ch, k] == i)) { ... : la

This is day one on R for me, I am trying to figure out how to do simple computations. For example I have a data set with 200 observations. I am trying to compute the mean and variance in r for 1:25 (first 25 observations); 1:50 (first 50 obs) ; 1:100th observation etc. Here is the dataset: Id value 1 2.2338 2 3.13597 3 1.98685 4 4.35593 5 2.43963 6 4.20262 7 3.12131 8 4.79583 9 3.13937 10

From: jypuppy36@hotmail.com To: r-help-bounces@r-project.org Subject: R programing help-newton iterations for the square root Date: Tue, 7 Dec 2010 12:00:01 -0800 NEWTON ITERATIONS FOR THE SQUARE ROOT Newton iterations to find the root of a real valued function f , i.e. a number x for which f (x) = 0, are of the form Example. To find the square root of a positive number y we can use

Hi r-users,   I hope somebody can help me with this code. I would like to solve for z values using newton iteration method.  I 'm not sure which part of the code is wrong since I'm not very good at programming but would like to learn.  There seem to be some output but what I expected is a vector of z values.  Thank you so much for any help given.   newton.inputsingle <-