search for: livia

Hi, I would like to generate a series in the following form (0.8^1, 0.8^2, ..., 0.8^600) Could anyone tell me how can I achieve that? I am really new to R. -- View this message in context: http://www.nabble.com/sequences-tf4019146.html#a11414836 Sent from the R help mailing list archive at Nabble.com.

Hello, I would like to fit a linear regression and when I use summary(), I got the following result: Call: lm(formula = weight ~ group - 1) Residuals: Min 1Q Median 3Q Max -1.0710 -0.4938 0.0685 0.2462 1.3690 Coefficients: Estimate Std. Error t value Pr(>|t|) groupCtl 5.0320 0.2202 22.85 9.55e-15 *** groupTrt 4.6610 0.2202 21.16 3.62e-14

Hi, I would like to plot a histogram with the following codes, and I would like to make the tick intervals smaller. I tried to add "lab=c(1,1,12)", but nothing changes. par(mfrow=c(1,1),font=1, cex=0.8) hist (data1, seq(-0.01,0.3,0.01),freq = FALSE, main="Figure1.",xlab="return",lab=c(1,1,12)) lines(density(data1), col="blue") Could anyone give me sone

Hi All, I am trying to make a loop for a function and I am using the following codes. "p" and "var" are some matrix obtained before. I would like to apply the function "gpdlow" for i in 1:12 and get the "returnlow" for i in 1:12. But when I ask for "returnlow" there are warnings and it turns out some strange result. for (i in 1:12){ gpdlow

Hi, I would like to use the following codes to plot the CDF for pareto distribution. Before doing this, I have plot the emperical one. x <- seq(1.6, 3, 0.1) lines(x,pgpd(x, 1.544,0.4477557,), col="red") Could anyone give me some advice whether the above codes are correct? Many thanks. -- View this message in context:

Hello, I would like to use the cbind() function to construct a matrix used in the middle of a function as following for (i in 1:1000) { b[i] <- function(cbind(a[[1]][[i]], a[[2]][[i]],a[[3]][[i]],...a[[67]][[i]])) } Is there an easy way of achieving this rather than "cbind" every column? -- View this message in context: http://www.nabble.com/cbind%28%29-tf4724534.html#a13507633

Hi Livia and everyone, Did you ever get a response on this question from last year (Jan 2008)? I am also looking for more explanatory documentation on the ui and ci parameters for the function constrOptim(). The examples provided in the R help and the full reference man...

I would like to fit a Pareto Distribution and I am using the following codes. I thought the fitted (fit1) should be the fitted value for the data, is it correct? As the result of the "fitted" turns out to be a single value for all. fit=vglm(ycf1 ~ 1, pareto1(location=alpha), trace=TRUE, crit="c") fitted(fit) The result is fitted(fit) [,1] [1,] 0.07752694

Hi, I have a series of return data, a and b are factors. I would like to build a matrix which contains each vector of "returns". I am thinking about something as following, but I guess there should be a sensible way of doing this. returns <- split(return, list(regimef, assetf)) cbind(returns[[1]], returns[[2]],...,returns[[n]]) Could anyone give me some advice? Many thanks. --

I have got a list named "filtered", I would like to construct alist named "fdata" as following: fdata <- cbind(matrix(unlist(filtered),ncol=28), myregime) If I try names(filtered), it gives all the correct name for each vector, but if I try names(fdata), it appears "filtered[[1]]" "filtered[[2]]" ..., How can I keep the name in "fdata"?

Hi, I would like to minimize the value of x1-x2, x2 is a fixed value of 0.01, x1 is the quantile of normal distribution (0.0032,x) with probability of 0.7, and the changing value should be x. Initial value for x is 0.0207. I am using the following codes, but it does not work. fr <- function(x) { x1<-qnorm(0.7,0.0032,x) x2=0.01 x1-x2 } xsd <- optim(0.0207, fr,

Hello all, I would like to plot the emperical CDF, normal CDF and pareto CDF in the same graph and I amusing the following codes. "z" is a vector and I just need the part when z between 1.6 and 3. plot(ecdf(z), do.points=FALSE, verticals=TRUE, xlim=c(1.6,3),ylim=c(1-sum(z>1.6)/length(z), 1)) x <- seq(1.6, 3, 0.1) lines(x,pgpd(x, 1.544,0.4373,-0.2398), col="red") y

Hello, I am using the following codes to plot a histogram and density line for x. For the density line, I just want it to show the two tails, eg, for x larger than 0.05 ans smaller than -0.05 hist (x, seq(-0.1,0.1,0.01),freq = FALSE) lines (density(x,bw="SJ"), x> 0.05 & x< (-0.05), col = "red") But is does not work, can anyone give me some advice? -- View this

Hello, I would like to perform a linear regression and the data is a list.e.g lm(list$abc~., data=list) or lm(abc~., data=list), which would give the same result. The problem is I would like to call the response variable in a more general form. What I try to achieve is sth like lm(list$(paste("a","b","c"))~., data=list), but it does not work. Could anyone give me

Hello everyone, I would like to a linear regression with the following code. lm(a[,"fquamsci"]~., data=a) a is a list with class "mts" "ts" , and "fquamsci" is the name of the response variable in a. I would like to do a linear regression of "fquamsci" to the rest of the variables. But it turns out the "fquamsci" is also included in

Dear members, I would like to pass the histogram settings to each subset of the dataframe, and generate a multiple figures graph. First, can anyone tell me how to generate a multiple figures environment? I am trying mfrow=c(2,4) and nothing appears. Secondly, I want to pass the following function in tapply() hist(x, freq=FALSE) lines(density(x), col="red") rug(x) how can I manage

Hi, I want to import some data from Access and I am using the following codes: testdb <- file.path("c/../db1") channel <- odbcConnect("testdb") sqlFetch(channel,"tbl",colnames = TRUE, rownames = FALSE) It comes out the error message: 1: [RODBC] ERROR: state IM002, code 0, message [Microsoft][ODBC Driver Manager] Data source name not found and no default

Hi, I am trying to write a function with the following codes and I would like it to return the values for "alpha beta para parab " seperately. Then I would like to use this funstion for "variable" with factor "a" and "b". But the result turns out to be a matrix with element like "Numeric,2" ... I guess they are just the values for

Hello, I would like to use the function EWMA() in the fMultivar Package and I have a series of data x, which is the returns series. Basically, I would like to get the variance estimation using EWMA. I am trying something like EWMA(x, lambda) and I have a couple of questions: Should x be the returns series or price series in my case? When I get the result, there are the same numbers of data

Hello everyone, I get some data in the following format and I would like to combine them to form a dataframe. The data is like: cbcname1 = 0.1, cbcname2= 0.2, cbcname3=0.3,... name1, name2, name2 are just some random names. I would like to achieve sth like: (cbcname1=0.1, cbcname2=0.2, cbcname3=0.3,......) I am using the following codes do.call(cbind,