Re: [R] make hyperlink in R

On May 7, 2013, at 11:15 PM, Bill Hyman wrote: > Dear all, > > Does anybody know how to make hyperlink in R? If I want to output > "http://www.r-project.org/"; with hyperlink, how can I do with plain text > 'http://www.r-project.org/'? I don't see any characters that need escaping. What are

Re: [R] How can I find negative items from a vector with a short command?

Thanks, but why does f[f=4] yield 16 instead of 4? > f [1] -24 -8 16 -32 64 -128 > f[f<0] [1] -2 -8 -32 -128 > f[f>0] [1] 4 16 64 > f[f=4] [1] 16 2013/5/8 Jorge I Velez > f [ f < 0 ] > > > On Wed, May 8, 2013 at 11:54 AM, jpm miao wrote: > >> Hi, >> >>I have a vector f

Re: [R] How can I find negative items from a vector with a short command?

Sorry. I just got it > f[f==4] [1] 4 2013/5/8 jpm miao > Thanks, but why does f[f=4] yield 16 instead of 4? > > > f > [1] -24 -8 16 -32 64 -128 > > f[f<0] > [1] -2 -8 -32 -128 > > f[f>0] > [1] 4 16 64 > > f[f=4] > [1] 16 > > > 2013/5/8 Jorge I Velez > >> f [ f < 0 ] >>

Re: [R] How can I find negative items from a vector with a short command?

Check ?"==" and try f[f == 4] More in section 2.4 at http://cran.r-project.org/doc/manuals/R-intro.pdf On Wed, May 8, 2013 at 5:14 PM, jpm miao wrote: > Thanks, but why does f[f=4] yield 16 instead of 4? > > > f > [1] -24 -8 16 -32 64 -128 > > f[f<0] > [1] -2 -8 -32 -128

Re: [R] Calculates the mean/median from grouped data in R?

see if this helps: library(actuar) data(gdental)# grouped data quantile(gdental) mean(gdental) Nuri On Wed, 8 May 2013, jpm miao wrote: Is there a function in R that calculate the mean and median for a grouped data? For example, a survey shows the oil price outlook in the future. How can

[R] Parsing XML to tree.

Hi All, I am struggling to parse a XML file that describes a tree. The XML file is present here: http://www.emouseatlas.org/emap/ema/theiler_stages/StageDefinition/Stage_xml/TS23.xml I want is simple list of parents id of each component. The output will look like Component = [7148 7149 7150 7

[R] Fitting AR(p) model

Hi all, I am trying to fit an AR(p) model to my variable res ar=ar(res, method="yule-walker") The output is : Order selected 0 sigma^2 estimated as 5.87e-06 Does this mean that the optimum value of p is 0? But when I go for summary(ar), I get the following: order 1 -none- numer

Re: [R] Fitting AR(p) model

Hello, You forgot the column names. The output is: Length Class Mode order 1 -none- numeric ar 0 -none- numeric var.pred1 -none- numeric x.mean 1 -none- numeric aic13 -none- numeric n.used 1 -none- numeric order.max 1

[R] R prints empty group on my figure!!

Hello list I am trying to solve a problem when drawing a figure related to the dataframe below. I draw V6 in Y axes and V3 as X axes. The data are grouped by V1 AND V2. I use this commande to make it: xyplot(cm[,6]~cm[,3],type="b",group=interaction(cm[,1],cm[,2],sep="/"), auto.key =list( title="Ma

Re: [R] R help for creating expression data of Differentially expressed genes

Hi Arun, I am still facing trouble as I can see the data output is identical for all rows when I am using this merge function. It seems that since in my data2 which I have provided I have not given you the exact genes I have. There are likely to be repeatations of ID in both the files but the pro

[R] Fuzzy rules definition (package "sets") from data.frame

Hi, I just discovered the package "sets" and its ability to perform fuzzy systems calculations. The example in the manual of fuzzyinference() gives an overview how to develop rules. However my "rules" are already available as data.frame and I'd like to convert them into the format that is needed

[R] gsub regex simplification

Dear all, I want to use gsub to change a vector of strings. Basically, I want to replace any dot by a space, remove the possibly appended ".f" and I want to capitalize each word. I did that by chaining multiple gsubs together, but I was wondering (for the sake of learning - maybe the current ve

Re: [R] R prints empty group on my figure!!

hello, I attach the figure generated by R: 2013/5/8 Adel ESSAFI > Hello list > I am trying to solve a problem when drawing a figure related to the > dataframe below. > I draw V6 in Y axes and V3 as X axes. The data are grouped by V1 AND V2. > I use this commande to make it: > > xyplot(cm[,6]~cm

Re: [R] gsub regex simplification

On Wed, May 8, 2013 at 5:08 AM, Thaler,Thorn,LAUSANNE,Applied Mathematics wrote: > Dear all, > > I want to use gsub to change a vector of strings. Basically, I want to > replace any dot by a space, remove the possibly appended ".f" and I want to > capitalize each word. I did that by chaining mul

Re: [R] Fuzzy rules definition (package "sets") from data.frame

Following my last mail, I found a simple solution using eval(parse)): df <- data.frame(depth=c("low","medium","high"),velocity=c("medium","medium","low"),suitability=c("low","medium","low")) df$rule <- paste("fuzzy_rule(depth %is% ",df[,"depth"]," && ","velocity %is% ",df[,"velocity"],", ","suitab

Re: [R] is.numeric () FALSE

Jeff, The dataset are 3D landmarks taken from a skull (Homo being the dataset of 3 species of the genus Homo) I am trying to perform Generalised Procrustes Distance and superimposition and as you can see, I am told x is not True. If I also perform is.numeric(Homo) I am told that it is FALSE. I hop

Re: [R] Vector allocation problem while trying to plot 6 MB data file

Thanks for your answer Uwe On Fri, 03 May 2013 23:36:24 +0200 Uwe Ligges wrote: > On 02.05.2013 14:37, Ramon Hofer wrote: > > > > I'm trying to analyse the network speed and used iperf to create a > > csv file containing the link test data. It's only about 6 MB big but > > contains about 40'00

Re: [R] R prints empty group on my figure!!

try using: group = factor(paste(cm[, 1], cm[, 2], sep = '/')) instead of group=interaction(cm[,1],cm[,2],sep="/") On Wed, May 8, 2013 at 4:25 AM, Adel ESSAFI wrote: > Hello list > I am trying to solve a problem when drawing a figure related to the > dataframe below. > I draw V6 in Y axes and

Re: [R] Vector allocation problem while trying to plot 6 MB data file

On 08.05.2013 13:38, Ramon Hofer wrote: Thanks for your answer Uwe On Fri, 03 May 2013 23:36:24 +0200 Uwe Ligges wrote: On 02.05.2013 14:37, Ramon Hofer wrote: I'm trying to analyse the network speed and used iperf to create a csv file containing the link test data. It's only about 6 MB

Re: [R] How does one set up logical functions?

Bruce, Try something like this df$limlo <- ifelse(df$BG>=100, df$BG-(df$BG*0.15), df$BG-15) df$limhi <- ifelse(df$BG>=100, df$BG+(df$BG*0.15), df$BG+15) plot(seq(df$BG), df$BG, ylim=range(df$limlo, df$limhi)) arrows(seq(df$BG), df$limlo, seq(df$BG), df$limhi, length=0.1, angle=90, code=3) Jean

[R] NaN-result from fuzzy_inference (package "sets") with certain input

Hi, I am trying to use the fuzzy_inference system on multiple input values. Basically I combine two variables (depth, velocity) into the variable suitability. Both depth and velocity have 3 trapezoid classes (verylow,medium,high) each. The consequent (suitability) has 4 fuzzy levels (none,low,medi

Re: [R] chronological season results assistance

Hi, May be this helps: funFluxSumSeason <- function(DF, FUN = sum){  month <- as.integer(format(DF$Date, format="%m"))  year <- format(DF$Date, format="%Y")  DF$season<- NA  DF$season[month %in% 10:12] <- paste(year[month %in% 10:12], "Fall")  DF$season[month %in% 1:3] <- paste(year[month %in% 1:3

Re: [R] How to calculate Hightest Posterior Density (HPD) of coeficients in a simple regression (lm) in R?

Richard Asturia gmail.com> writes: > > Hi! > > I am trying to calculate HPD for the coeficients of regression models > fitted with lm or lmrob in R, pretty much in the same way that can be > accomplished by the association of mcmcsamp and HPDinterval functions for > multilevel models fitted wit

Re: [R] How to calculate Hightest Posterior Density (HPD) of coeficients in a simple regression (lm) in R?

Thanks a lot for your reply! That is interesting function. And you are completely right. Of course it may be an overkill to use hpd with lm... But anyway, I am precisely interested in comparing CI and HPD for a model fitted with lm and then do the same comparison for the same model fitted with lmr

Re: [R] State space models with regime switching

On 05/06/2013 04:28 PM, David Hoppe wrote: Hello everyone, I'm new to this mailing list, but i hope this is the right place to post my question. I'm trying to do some time series analysis with state space models in R. So far I used the packages dse and dlm. I was wondering if there is a package,

[R] gam (mgcv), multiple imputation, f-stats/p-values, and summary(gam)

Dear All, I'm using gam for a project that involves multiple imputation, and it has led me to a question about how f-statistics/p-values work in gam. Specifically, how do the values in summary(gam) get generated? As is made clear by the dumb example below, I'm manipul;ating gam objects to r

Re: [R] is.numeric () FALSE

You have still left out important information. You created Homo as a data.frame but the next command arrayspecs() assumes Homo is a matrix. There appears to be an initial column of characters that indicates the landmark? What does str(Homo) produce after read.table()? We need the steps you used to

Re: [R] gam (mgcv), multiple imputation, f-stats/p-values, and summary(gam)

?summary.gam ## The Help page Since the Help page is presumably not enough, why don't you look at the code?? R is open source. summary.gam ## at the prompt -- Bert On Wed, May 8, 2013 at 7:12 AM, Andrew Crane-Droesch wrote: > Dear All, > > I'm using gam for a project that involves multiple i

Re: [R] gam (mgcv), multiple imputation, f-stats/p-values, and summary(gam)

For smooths the method is described in Wood 2013 On p-values for smooth components of an extended generalized additive model, Biometrika 100(1),221-228 http://biomet.oxfordjournals.org/content/early/2012/10/18/biomet.ass048.full.pdf+html best, Simon On 08/05/13 15:12, Andrew Crane-Droesch wr

Re: [R] Parsing XML to tree.

Hi, On May 8, 2013, at 3:43 AM, avinash sahu wrote: > Hi All, > > I am struggling to parse a XML file that describes a tree. The XML file is > present here: > http://www.emouseatlas.org/emap/ema/theiler_stages/StageDefinition/Stage_xml/TS23.xml > I want is simple list of parents id of each compo

Re: [R] rowspan and readHTMLTable

Sorry to answer my own question - I guess here's one way to read this table. Other suggestions are still welcome. Chris -- x<-htmlParse(" abX YZ c ") # split by rows z <- getNodeSet(x, "//tr") # create empty data.frame - probably not the best solution... t1<- data.frame(matrix(NA, nrow

Re: [R] Some unrelated questions.

R gives you facilities for doing this (using '...' and 'missing') without passing in a data.frame. For example, > foo <- function(arg1,arg2,...) { if (missing(arg1)) cat('missing arg1\n'); > print(list(...)) } > foo(1,2,arg3=3,arg4=4,arg5=5) $arg3 [1] 3 $arg4 [1] 4 $arg5 [1] 5 > foo(arg2=2

Re: [R] chronological season results assistance

Hi Arun, thank you. I am well. Yes, that worked perfectly. Irucka <-Original Message-> >From: arun [smartpink...@yahoo.com] >Sent: 5/8/2013 11:39:54 AM >To: iruc...@mail2world.com >Cc: r-help@r-project.org >Subject: Re: chronological season results assistance > > > >Hi, >I am fine. Ho

[R] Estimating confidence intervals from several simulations

Hi, I have 1000 simulations each arranged beneath another. The attached data is on monthly basis. First column is number of simulation, second is year and the rest, months. Problem:  1) Take sim 1 for example to be my original data from 1961-2005 with twelve months. You can delete the years. I

Re: [R] Tranfer data from exel into R

Hi, how are you? You are welcome for my assistance. I'm sorry, but I made a mistake in my previous post. Is your .csv file in your working R directory? If so, then H<-read.csv("Book.csv", header=TRUE, sep=",") should work. I am assuming that your data is comma separated and that the names of

Re: [R] chronological season results assistance

Hi Arun, how are you? Thank you for your assistance. I'm sorry, but I made a mistake in what I was seeking. For this example, this is how I would like the results to look: 2000 Fall 2001 Winter 2001 Spring 2001 Summer 2001 Fall 2002 Winter 2002 Spring 2002 Summer 2002 Fall 2003 Winter 2003 Sprin

[R] Using the Sun Grid Engine with R-3.0

Dear all In previous R versions I ran R jobs on our cluster managed by SUN Grid Engine using the Rsge package. With the R-3.0 release, the package Rsge has been discontinued. What would be the equivalent way of running jobs on the cluster using current packages? I understand that the paralle

[R] scatter plot matrix with different x-y variables

Dear list-members, I wonder if there is a way of creating a scatter plot table/grid with different variables on the y-axis compared to the x-axis? Something like this: A*** B*** C*** -XYZ I know that you can create scatter plot matrix with the same variables on the y-axis as on the x-axis, lik

[R] generate gaussian

Dear Sir, My name is Iut Tri Utami. i am beginning user. I have a problem about generate data in R. It consists of one disk generated by a Gaussian N(0, 0.167) and one ring generated by a Gaussian N(R, 0.1). The mean R was generated from its polar coordinates. The angle was drawn from a uniform d

[R] Correctly Setting New Seed

Could someone please suggest a method to store the current random seed. I'm having trouble understanding how to correctly use set.seed and .Random.seed. Specifically, I have the following code that crashes: set.seed(seed) for(i in 1:10){ print( runif(1)) } To get around this I need to

Re: [R] Correctly Setting New Seed

"Crashes" is not a clear description. When I execute the statement: set.seed(seed) I get the error Error in set.seed(seed) : object 'seed' not found which makes sense because there is no variable defined that is called "seed". Choose a literal value like 42: set.seed(42) or define the variab

Re: [R] Very basic statistics in R

On 5/6/2013 7:02 AM, arun wrote: stErr<- sd(vec1)/sqrt(length(vec1)) Or possibly, stErr<- sd(vec1)/sqrt(!is.na(vec1)) -- Robert W. Baer, Ph.D. Professor of Physiology Kirksille College of Osteopathic Medicine A. T. Still University of Health Sciences Kirksville, MO 63501 USA

Re: [R] generate gaussian

Someone probably could help but you did not tell us what the problem is. http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example John Kane Kingston ON Canada > -Original Message- > From: triutami@gmail.com > Sent: Thu, 9 May 2013 00:03:25 +0700 > To:

Re: [R] Correctly Setting New Seed

By crash, I mean that it shuts the program and no error message is provided. However, my question is more general - how to get the two versions of code below to return the same string of random numbers. The code provided will run without incident and is only provided to make the problem clear.

Re: [R] scatter plot matrix with different x-y variables

Several ways but I find that grid.arrange in the gridExtra package combined with ggplot2 to draw the actual graphs works nicely. You might also want to look mfcol or mfrow in ?par ### quick and dirty example of ggplot2/gridExtra library(ggplot2) library(gridExtra) mdat <- data.frame(a= 1:5

[R] PPCOR: Semipartial Correlation & Regression weights

In linear regression, regression weights of x1 on Y given x2 and x3 should be mathematically identical to the semipartial correlations between x1 and Y, given x2 and x3. However, I do not obtain identical results, so apparently I'm doing something wrong in R. Data preparation: data<-read.csv("fil

Re: [R] Correctly Setting New Seed

I was wondering about that code. It seemed to run alright for me. We really need to see the code and hopefully the data, or a,sample thereof, that is causing the actual crash. Prefereably cut out anything that is not causing the problem and ship us the bare minimum of code. http://stackover

Re: [R] how to read numeric vector as factors using read.table.ffdf

Please read the posting guide. Your question is far from clear. First you're apparently unhappy because character vectors are being converted to factors. Then later you ask how to tell a function that some numeric datas should be considered factors. Which is it that you want? A simple short exampl

Re: [R] Very basic statistics in R

On May 8, 2013, at 20:38 , Robert Baer wrote: > On 5/6/2013 7:02 AM, arun wrote: > stErr<- sd(vec1)/sqrt(length(vec1)) > > Or possibly, > > stErr<- sd(vec1)/sqrt(!is.na(vec1)) > You probably intended sqrt(sum(!is.na(...))). -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen B

Re: [R] Correctly Setting New Seed

> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > project.org] On Behalf Of Ellerbe, Caitlyn Nicole > Sent: Wednesday, May 08, 2013 11:51 AM > To: Ellerbe, Caitlyn Nicole; r-help@R-project.org > Subject: Re: [R] Correctly Setting New Seed > > By crash, I

[R] ecdf --- title suggestion and question

dear R-experts---first, a suggestion to martin: the ecdf() function could have an optional parameter to set the title. by looking at str(), I see the plot title is set in an attr named "call". i.e., I can reset it as ee <- ecdf( rnorm(25 ) ) attr(ee,"call") <- "my own title" plot(ee) alas, I ca

Re: [R] ecdf --- title suggestion and question

Hi Ivo, On Wed, May 8, 2013 at 1:37 PM, ivo welch wrote: > dear R-experts---first, a suggestion to martin: the ecdf() function > could have an optional parameter to set the title. by looking at > str(), I see the plot title is set in an attr named "call". i.e., I > can reset it as > > ee <- ecd

Re: [R] ecdf --- title suggestion and question

ooops...never mind. I mixed up "title" and "main" as options. Ivo Welch (ivo.we...@gmail.com) On Wed, May 8, 2013 at 1:54 PM, Steve Lianoglou wrote: > Hi Ivo, > > On Wed, May 8, 2013 at 1:37 PM, ivo welch wrote: >> dear R-experts---first, a suggestion to martin: the ecdf() function >> cou

Re: [R] Parsing XML to tree.

Hi All, I am using XML package. Even this type of simple parse is not giving intended result. tissue.tree <- xmlTreeParse(" http://www.emouseatlas.org/emap/ema/theiler_stages/StageDefinition/Stage_xml/TS23.xml";, handlers=list( anatomy=function(x,attr) {x},

Re: [R] Very basic statistics in R

Hello, Or, to make it complete, use the na.rm argument. sd(vec1, na.rm = TRUE)/sqrt(sum(!is.na(vec1))) Rui Barradas Em 08-05-2013 20:55, peter dalgaard escreveu: On May 8, 2013, at 20:38 , Robert Baer wrote: On 5/6/2013 7:02 AM, arun wrote: stErr<- sd(vec1)/sqrt(length(vec1)) Or possib

Re: [R] Parsing XML to tree.

On Wed, May 8, 2013 at 3:43 AM, avinash sahu wrote: > Hi All, > > I am struggling to parse a XML file that describes a tree. The XML file is > present here: > http://www.emouseatlas.org/emap/ema/theiler_stages/StageDefinition/Stage_xml/TS23.xml > I want is simple list of parents id of each compone

Re: [R] Correctly Setting New Seed

This approach gives the same results: > set.seed(1) > for (i in 1:10) print(runif(1)) [1] 0.2655087 [1] 0.3721239 [1] 0.5728534 [1] 0.9082078 [1] 0.2016819 [1] 0.8983897 [1] 0.9446753 [1] 0.6607978 [1] 0.629114 [1] 0.06178627 > > set.seed(1) > for (i in 1:5) print(runif(1)) [1] 0.2655087 [1] 0.372

Re: [R] Very basic statistics in R

On May 8, 2013, at 12:55 PM, peter dalgaard wrote: > > On May 8, 2013, at 20:38 , Robert Baer wrote: > >> On 5/6/2013 7:02 AM, arun wrote: >> stErr<- sd(vec1)/sqrt(length(vec1)) >> >> Or possibly, >> >> stErr<- sd(vec1)/sqrt(!is.na(vec1)) >> > > You probably intended sqrt(sum(!is.na(...)))

[R] NMDS with missing data?

Hi, I'm trying to run NMDS (non-metric multidimensional scaling) with R vegan (metaMDS) but I have a few NAs in my data set. I've tried to run it 2 ways. The first way with my entire data set which includes variables such as ID, sex, exposure, treatment, sodium, potassium, chloride mydata.mds

[R] R unable to access internet to install packages

Hi, I am having some trouble trying to install a number of packages in R on my work computer. An example of the code I am using is as follows: install.packages("ggplot2",dependencies=TRUE) An example of the error messages I am getting is: Warning: unable to access index for repository http://

[R] [newbie] *apply to matching elements of n arrays?

How to apply a function to all elements with the same indices in multiple arrays? E.g.: I have two spatial grids defined the same way (i.e., same number of rows and columns--and dimensions, both 2D). Wherever both * the value of an element i,j in the first grid is NA * the value of element i,j i

Re: [R] R unable to access internet to install packages

a) You really should upgrade your version of R. b) You have not indicated clearly what your operating system is as the Posting Guide requests of you (e.g. sessionInfo()), but if it is Windows then you should look at the R for Windows FAQ (section 2.19). --

Re: [R] R unable to access internet to install packages

Hello, About problem for R to download files: http://cran.r-project.org/bin/windows/rw-FAQ.html#The-Internet-download-functions-fail_002e The repertory "2.9" for the second mirror doesn't exist. Hope this helps, Pascal On 05/09/2013 09:30 AM, Alex Collier wrote: Hi, I am having some trouble

Re: [R] [newbie] *apply to matching elements of n arrays?

Not particularly interested in "points" from you. Would like reproducibility from you [1], including dput of before and after data, and (inefficient) code of course. For loops are part of R.. If your data structure is not set up to take advantage of vectorization, then for loops are generally

[R] Choice of statistical test (in R) of two apparently different distributions

I have the following list of data each has 10 samples. The values indicate binding strength of a particular molecule. What I want so show is that 'x' is statistically different from 'y', 'z' and 'w'. Which it does if you look at X it has more values greater than zero (2.8,1.00,5.4, etc) than othe

Re: [R] Choice of statistical test (in R) of two apparently different distributions

Hello, 1) Why 'x' should be statistically different from others? 2) 'y' looks to be bimodal. The mean is not an appropriate measurement for this kind of distribution. Regards, Pascal On 05/09/2013 10:29 AM, Gundala Viswanath wrote: I have the following list of data each has 10 samples. The

[R] Inputting white noise values for ARMA prediction

Hi All, I have the following time series model for prediction purposes *Loss_t = b1* Loss_(t-1) + b2*GDP_t + b3*W_(t-1)* where W_t is the usual white noise variable. So this is similar to ARMA(1,1) except that it also contains an extra predictor, GDP at time t. I have only 20 observat

[R] error with lhoat() function from the package hydroPSO

Dear all, I just installed the package hydroPSO and when I try to run the example of the lhoat() function I got this error: nparam <- 5 lhoat( fn=sphere, lower=rep(-100,nparam), upper=rep(100,nparam), control=list(N=10, f=0.1, write2disk=TRUE) ) *Error in fn.n

Re: [R] Bootstrapped 1-sided confidence intervals

Great! Thanks to all for your assistance. Regards Janh On Tue, May 7, 2013 at 11:37 PM, Pascal Oettli wrote: > Hello, > > You already asked that question on May 7, 2013. And David Winsemius > already responded to you: > https://stat.ethz.ch/**pipermail/r-help/2013-May/**353044.html

Re: [R] [newbie] *apply to matching elements of n arrays?

Is the following what you want? It works for vectors or arrays of any number of dimensions. It assumes that the dimensions of the grids are the same. > grid1 <- c(1,2,NA,3,NA) > grid2 <- c(101,102,103,104,NA) > shouldCopy <- is.na(grid1) & !is.na(grid2) > grid1[shouldCopy] <- grid2[shouldCopy]

Re: [R] R unable to access internet to install packages

Alex For using a proxy server configuration with R please ensure that you use the internet2 ddl to connect to the internet. It can be set in an R session by using setInternet2(TRUE). Besides all of these you should really set up a new version of R. Regards The woods are lovely, dark and deep But

Re: [R] R unable to access internet to install packages

also try this command to set the proxy Sys.setenv(http_proxy="http://proxy_server_ip :port_for_proxy") The woods are lovely, dark and deep But I have promises to keep And miles to go before I sleep And miles to go before I sleep - On Thu, May 9, 2013 at 7:50 AM, Nilesh Gupta wrote: > Alex >

[R] Re-start R issue on win7 [SEC=UNOFFICIAL]

Hi all, It is really annoying. The first time starting R, it works, but I cannot re-start R unless re-start my PC. Operation system is Win7. Anyone knows how to deal with issue? Cheers, Xiaoping ** IMPORTANT: This e-mail is

[R] ARMA(p,q) prediction with pre-determined coefficients

I have the following time series model for prediction purposes *Loss_t = b1* Loss_(t-1) + b2*GDP_t + b3*W_(t-1)* where W_t is the usual white noise variable. So this is similar to ARMA(1,1) except that it also contains an extra predictor, GDP at time t. I have only 20 observations on eac

Re: [R] Re-start R issue on win7 [SEC=UNOFFICIAL]

Have not experienced that. --- Jeff NewmillerThe . . Go Live... DCN:Basics: ##.#. ##.#. Live Go... Live: OO#.. Dead: OO#.. Playi

Re: [R] Re-start R issue on win7 [SEC=UNOFFICIAL]

I have not experienced that, but have you tried to see if killing the processes helps? On Windows CTRL+ALT+DEL usually brings up the system monitor or whatever it's called, where you can see programs and processes currently running. If your R gets stuck in a loop then you may see some refer