search for: datax

...=2,sd=1) x3 <- rnorm(n=50,mean=2,sd=1) x4 <- rnorm(n=50,mean=2,sd=1) y <- 1+ 2*x1+4*x2+3*x3+2*x4+e x2[1] = 10 #influential observarion y[1] = 10 #influential observarion data.x <- matrix(c(x0,x1,x2,x3,x4),ncol=5) data.y <- matrix(y,ncol=1) data.k <- cbind(data.x,data.y) dataX <- data.k[,1:5] dataY <- data.k[,6] theta <- function(data) { B.cap <- solve(crossprod(dataX)) %*% crossprod(dataX,dataY) P <- dataX %*% solve(crossprod(dataX)) %*% t(dataX) Y.cap <- P %*% dataY e <- dataY - Y.cap dX <- nrow(dataX) - ncol(dataX) var.cap <- crossprod(e) /...

...#influential observarion X <- matrix(c(x0,x1,x2,x3,x4),ncol=5) Y <- matrix(y,ncol=1) Design.data <- cbind(X, Y) for (j in 1:nrow(X)) { result <- vector("list", ) for( i in 1: 3100) { data <- Design.data[sample(50,50,replace=TRUE),] ##### and unique(Design.data.....) dataX <- data[,1:5] dataY <- data[,6] B.cap.simulation <- solve(crossprod(dataX)) %*% crossprod(dataX, dataY) P.simulation <- dataX %*% solve(crossprod(dataX)) %*% t(dataX) Y.cap.simulation <- P.simulation %*% dataY e.simulation <- dataY - Y.cap.simulation dX.simulation <- nrow(data...

Just found because an old (written for S+) function of mine did label plots wrongly with R. Example --- inspired from example(qqnorm) --- data(precip) qqnorm(precip, ylab = "Precipitation [in/yr] ...") qqnorm(precip, ylab = "Precipitation [in/yr] ...", datax = TRUE) this is all "fine" -- `datax' is a switch that just switches the axes but keeps the meaning of "ylab labels the data axis" (i.e, when changing datax from FALSE to TRUE, you don't have to switch xlab and ylab in your call). Unfortunately, S+ (6.1) does this di...

....txt file). I wish to go for filtering this 50 rows for any one of the six coloumns satisfying the value >= 64. I need to have a final table with rows having >= 64 value in any one of the six coloumns and the rest could be <=64. For this purpose I use the following R code; --------------- datax<-read.table("filter_test.txt",row.names=1,sep="\t",header=TRUE,dec = ".",as.is =TRUE,na.strings = "NA", colClasses = NA,check.names = FALSE,strip.white = FALSE, blank.lines.skip = TRUE, allowEscapes = FALSE, flush = FALSE,encoding = "unknown")...

I am a beginner in R, so please don't step on me if this is too simple. I have two data sets datax and datay for which I created a qqplot qqplot(datax,datay) but now I want a line that indicates the perfect match so that I can see how much the plot diverts from the ideal. This ideal however is not normal, so I think qqnorm and qqline cannot be applied. Perhaps you can help? Ralf

...ges the types to symbol and if there is a way to get the actual data pointed to by the symbol? (I have tried eval but that does not work and most functions just treat a symbol as a string). ( An example which shows the type conversion is given below with the output - the key is the following "dataX=data" which makes the object data passed as a symbol and not the actual data). Many thanks Tom ====EXAMPLE CODE==== data=c(1,2,3,4,5,6,7,8,9); x <- function( ... ) { args <- list(); extras <- match.call(expand.dots = FALSE)$...; for( i in names(extras) ) {...

Buenas tardes, quedarme con las muestras de una BD (data) que están presentes en otra (datax), cuando se tiene una variable que nunca se repite (Key) es fácil: data <- subset(data,data$Key %in% datax$Key). Mi problema es cuando la exclusividad viene dada por dos variables. P.e., las coordenadas de un mapa: lon y lat. ¿Como puedo quedarme con las muestras de una df cuya lon y lat son igu...

...;- data[data$x3 < 5, ] lm(formula = formula, data = data, subset = counter == 1, weights = data$x2) # lm(formula = formula, data = data, subset = counter == 1, weights = data$x2)$weights } my.lm2(formula, datai) ################## ## The following all work in both R-1.2.0 and R-1.1.1 datax <- datai[datai$x3 < 5 & datai$counter==1,] lm(formula, data = datax, weights = x2) lm(formula, data = datax, weights = datax$x2) lm(formula, data = datai, subset = (counter == 1 & x3 < 5), weight = x2) lm(formula, data = datai, subset = (datai$counter == 1 & datai$x3 < 5), w...

Dear all, I am trying to my dataframe for the PLS analysis using the PLS package. However I have some trouble generating the correct dataframe. The main problem is how to use one name to represent several columns in the dataframe. The example dataframe in PLS package is called "sensory". I cannot directly read the data file since it's a binary file. If I use

...1 1 0 1 1 1 . . . . . . 1 0 1 1 1 1 0 1 1 1 . . . . . . 1 1 1 1 On this dataframe (here named datax), i perform LCA as follows: datax.int=datax+1 ### poLCA can only analyze 1,2,... f<-as.formula(paste("cbind(", paste(colnames(datax.int), collapse = ","), ")~1")) #all items are dependent variables fit<-list() #collect fits Kmax=5 #maximum nr of classes bic=rep...

Hi everybody I'm rather new to running wine so this might be a problem that is easy to correct. Anyways here goes: I have been trying to install a program called datax using Setup.exe on an original cd-rom. I have done some debugging on my own without any luck. My installation is able to run programs like Textpad 4 and notepad.exe. It's an fake windows installation on my linux partition. I'm using wine 20030911 on a AMD Duron 1,3 GHz, 256 MB ram, a Tri...

Hi, there. Is there any function in R to plot the probability-probability plot (PP plot)? Suppose I am testing some data against normal. Thanks. Yulei $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ Yulei He 1586 Murfin Ave. Apt 37 Ann Arbor, MI 48105-3135 yuleih at umich.edu 734-647-0305(H) 734-763-0421(O) 734-763-0427(O) 734-764-8263(fax) $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

Hi all, I have just started using R to produce qqnorm plots. I am trying to switch the x and y axes so that the theoretical values are plotted on the y axis and my data on the x axis. Can anyone help me with this? Thanks -- View this message in context: http://www.nabble.com/Flipping-axes-of-qqnorm-tp23248007p23248007.html Sent from the R help mailing list archive at Nabble.com.

two ARIMA models, both have several bars signicant in ACF and PACF plots of their residuals, but when run Ljung Box tests, both don't show any significant correlations... however, one model has p-value that is larger than the other model, based on the p-values, can I say the model with larger p-values should be better than the model with smaller p-values? [[alternative HTML version

Hi, I'm trying to run an accounting program called Mamut DaTax. Running the installation part is not a problem, but when entering the program itself all keys are dead but the following: ?????? I'm using a Norwegian keyboard, and wine version 1.3 under Kubuntu 10.04. Any suggestions for how to approach? Thanks :-) Best regards, Erik

Hi, I am trying to fit a multivariate time series model using DCC GARCH model and forecast it. The data looks like this: > head(datax) x vibration_x Speed 1 2017-05-16 17:53:00 -0.132 421.4189 2 2017-05-16 17:54:00 -0.296 1296.8882 3 2017-05-16 17:55:00 -0.572 0.0000 4 2017-05-16 17:56:00 -0.736 1254.2695 5 2017-05-16 17:57:00 0.000 0.0000 6 2017-05-16 17:58:00 0.00...

Hi, I am completely new to GARCH models and trying to fit a multivariate time series model using DCC GARCH model and forecast it. The data looks like this: > head(datax) x vibration_x Speed 1 2017-05-16 17:53:00 -0.132 421.4189 2 2017-05-16 17:54:00 -0.296 1296.8882 3 2017-05-16 17:55:00 -0.572 0.0000 4 2017-05-16 17:56:00 -0.736 1254.2695 5 2017-05-16 17:57:00 0.000 0.0000 6 2017-05-16 17:58:00 0.00...

...dataframe #######x1 <- c(1,NA,NA,NA,NA,1,2,2)x2 <- c(1,NA,NA,NA,NA,1,2,1)x3 <- c(1,NA,1,1,1,1,1,5)x4 <- c(1,NA,NA,NA,NA,NA,NA,5)x <- rbind(x1,x2,x3,x4)test <- rowSums(is.na(x)) ## count numer of NA's in rowx <- cbind(x, test) ## add row NA count to datax <- data.frame(x) ## make dataframe # FOR LOOP to apply across all rows of dataframe -- for(i in 1:nrow(x)) {if(x[i,9] > 4) {x[i,1:7] <- NA} else { x[i,1:7] <- x[i,1:7]i = i+1}} ##################################################### [[alternative HTML version deleted...

Does R do cumulative frequency distribution plots? -- Tommy E. Cathey, Senior Scientific Application Consultant High Performance Computing & Scientific Visualization SAIC, Supporting the EPA Research Triangle Park, NC 919-541-1500 EMail: cathey.tommy at epa.gov My e-mail does not reflect the opinion of SAIC or the EPA. Federal Contact - John B. Smith 919-541-1087 - smith.johnb at epa.gov

...ing that this little stunt be rather be hidden and automatically effected within "qqline" proper, whenever par('ylog') is TRUE. I thought about providing a patch, as "qqline" is so small. Yet it would be more noise than useful, as I'm not familiar with the "datax" argument usage, which should probably be addressed as well. Here is the data, in case useful: freq <- as.integer(c(33, 79, 21, 436, 58, 18, 1106, 498, 1567, 393, 2, 104, 50, 67, 113, 76, 327, 331, 196, 145, 86, 59, 12, 215, 293, 154, 500, 314, 246, 587, 85, 23, 323, 3, 13, 576, 29,...