similar to: Formula variable help

Hello all, Probably my concepts about the data.frame and matrix and array in R are not clear, I need some clarification to help me understand them better. >M <- read.table("test1.csv",sep=",",row.names=NULL,header=T) gives me: M as M1 M2 M3 M4 M5 M6 M7 M8 M9 M10 1 9 11 14 15 18 20 20 20 20 20 2 3 4 8 9 11 12 14 15 15 15 3 4 5 8 8 9 9 9 9 9 9 4 4

Hi, I am trying to find an efficient way of applying a linear regression model to different factor combinations in a data frame. I want to obtain the output with minimal or no use of loops if possible. Please let me know if this query is unclear. Thanks, Murtaza

I am analyzing data from an experiment with two factors: Carbon (+/-) and O3 (+/-), with 4 replicates of each treatment, and 4 harvests over a year. The treatments are assigned in a block design to individual Rings. I have approaches this as a repeated measures design. Fixed factors are Carbon, O3 and Harvest, with Ring assigned as a random variable. I have performed repeated measures analysis

Can someone enlighten me on why the following doesn't work? setwd('C:/Temp/R') d <- rep(1:53,2) (s <- ts(d, frequency=53, start=c(2000,10))) n <- length(s) k <- n%/%3 for(i in (n-k):n) { st <- c(start(s)[1] + (start(s)[2] + i)%/%frequency(s), (start(s)[2] + i) %% frequency(s)) ed <- c(start(s)[1] +

Hi to all the people again, I was writting a simply function in R, and wish to collect the results in a excel file. The work goes as follows, Ciervos<-function(K1, K0, A, R,M,Pi,Hembras) {B<-(K1-K0)/A T1<-(R*Pi*Hembras-M*Pi+B)/(Pi-M*Pi+R*Pi*Hembras) P1<-Pi-B R1<-P1*Hembras*R M1<-P1*M T2<-(R1-M1+B)/(P1-M1+R1) P2<-P1-B R2<-P2*Hembras*R M2<-P2*M

Hi All, How do i calculate KMO for a dataset? *Dataset:---------------------* m1 m2 m3 m4 m5 m6 m7 m8 1 2 20 20 2 1 4 14 12 2 9 16 3 5 2 5 5 15 3 18 18 18 13 17 9 2 4 4 7 7 2 12 2 11 11 11 5 7 8 5 19 5 2 20 18 6 7 4 7 4 7 9 3 3 7 5 5 5 12 5 13 13 12 8 6 6 4 3 5 17 17 16 9 12 12 4 2 4 4 14 14 10 5 14

Hi all! I?ve got some problems with the function gam (library mgcv). For some models I get the error message : Error: no valid set of coefficients has been found:please supply starting values In addition: Warning message: NaNs produced in: log(x) This is a shortened code I used: gam(y ~ M1 + M3 + M4 + M5 + M6 + sex + M1*M3 + s(age), family=Gamma(link ="identity"), weights=days) If

Hi Just upgraded smoothly to long overdue Aptana M7 Beta which was released on 25/9/2007. FYI, the best way to upgrade from July''s M5 (M6 was skipped) is as follows :- 1) Click Help/Software Upgrades/Finad and Install 2) Manually adding the following update site to the update manager: http://update.aptana.com/update/rails/beta/3.2/site.xml

Hi! I am new to looping and R in general; and I have sent waaaay to much time on this one problem and am about a hair away from doing it manually for the next two days. So, there is a package that while calculating the statistic creates lists (that look like matrices) in the background. Each item (there are 10 items) has one of these ‘matrix looking list’ that I need to extract data from. The

At 10:37 PM 5/1/00 -0400, E. S. Venkatraman wrote: >I have the following problem. I have survival data (time, status) along >with several covariates (X1, X2,..., Xn). I want to fit a Cox model for >each of the covariate (univariately) and obtain the fitted probability of >survival at a fixed time point t0 and covariate value Xi0. I tried to do >this in a for loop where the index

Dear all, I am new to plotting Taylor Diagram using plotrix package within R, hence this post. I have written a script which plots Taylor Diagram with one reference and 7 model values. However the font size, line width and line type are not clear when saving the diagram as a jpeg file. I tried the functions lty, lwd and font but no apparent change. I am attaching the script here. Any help would

Dear R helpers,   Following is a R script I am using to run the Fast Fourier Transform. The csv files has 10 columns with titles m1, m2, m3 .....m10.     When I use the following commands, I am getting the required results. The probelm is if there are 100 columns, it is not wise to define 100 commands as fk <- ONS$mk and so on. Thus, I need some guidance to write the loop for the STEP A and

I am used to java (well, i dont remember it really well, but anyway) I have having a really difficult time making simple loops to work. I got the following to work: ## ##Creates objects Ux1, Ux2, Ux2 etc. that all contain n numbers in a random distribution ## m<-c(m1,m2,m3,m4,m5,m6,m7,m8,m9,m10)#these are defined as numbers (means) v<-c(v1,v2,v3,v4,v5,v6,v7,v8,v9,v10)#these

I ran the following R script under both Linux and Windows, and got 2 different results. Linux R version 2.7.1 and Windows R version 2.7.2. > library(FactoMineR) >x1=read.table("freqtest.txt",header=TRUE) >xrcc2=x1[,1:8] >p1=PCA(xrcc2, graph=FALSE) >p1$var freqtest.txt file lines of text : M1 M2 M3 M4 M5 M6 M7 M8 -1 -1 -1 -1 -1 -1 -1 -1 0 0 0 0 -1 -1 1 1 -1 -1 -1 -1 -1 -1

Dear All I am reposting because I my problem is real issue and I have been working on this. I know this might be simple to those who know it ! Anyway I need help ! Let me clear my point. I have huge number of datapoints plotted using either base plot function or xyplot in lattice (I have preference to use lattice). name xvar p 1 M1 1 0.107983837 2 M2 11

Dear All: I am interested in multinomial logit models (function multinon, library nnet) but I'm having troubles in choose whether to define the predictors as factors or not. I had posted earlier this example (thanks for the reply ronggui): worms<- data.frame(year= rep(2000:2004, c(3,3,3,3,3)),age=rep(1:3,5),

This email is also sent to cluster-devel at redhat.com. Since this issue is about both dlm and ocfs2, I send the email here to look for help from upstream. This is an already known issue. on ocfs2 with user space cluster stack, run the test script from http://people.redhat.com/~teigland/make_panic on the mounted ocfs2 volume from 2 nodes simultaneously, the access to ocfs2 volume on both nodes

HI, I have a scenario where a filesystem i.e /fs13/m6.uecommon.1/wmt8/utils is residing on a fileserver i.e. fs13. Contents of this filesystem are as follows: # ls -l /fs13/m6.uecommon.1/wmt8/utils lrwxrwxrwx ........ include-> /fs12/m6.uecommon.1/utils/include lrwxrwxrwx ........ tool_utils-> /fs12/m6.uecommon.1/utils/tool_utils The filesystems being pointed to i.e

Hey everybody, i am using the rugarch-package and its great! I have a pretty easy problem, but i just dont get it, so thanks if you can help me. Normally i use: / data(DATANAME) spec = ugarchspec() fit = ugarchfit(data = x[,1], spec = spec) fit slotNames(fit) names(fit at fit) coef(fit) infocriteria(fit) likelihood(fit) nyblom(fit) signbias(fit) head(as.data.frame(fit)) head(sigma(fit))

I've got a complicated semi-parametric model that I'm fitting with mgcv. I start with a model based on theory. Its got lots of interaction terms. I want to winnow it down: removing each interaction term or un-interacted main effect one by one, checking the AIC, and retaining the model that gives me the lowest AIC. I then want to repeat the procedure on the retained model. Here is