Displaying 20 results from an estimated 4000 matches similar to: "Formula variable help"
2009 Apr 22
4
read.table or read.csv without row index?
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
2008 Dec 26
3
lm() with same formula but different column/factor combinations in data frame
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
2006 Nov 14
2
Repeated measures by lme and aov give different results
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
2011 Nov 08
3
window?
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] +
2011 Aug 13
1
Own R function doubt
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
2009 Sep 30
1
How to calculate KMO?
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
2005 Feb 14
1
gam(mgcv) starting values
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
2007 Sep 28
4
Aptana Radrails M7 Beta - Best Way to Upgrade
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
2012 Dec 13
1
How do I make a loop to extract a column from multiple lists and then bind them together to make a new matrix?
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
2000 May 02
2
Variable names in model formula
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
2012 Jan 30
1
about changing line type and line width in Taylor Diagram
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
2009 May 27
1
How to write a loop?
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
2007 Apr 04
2
Newbie: Simple loops: complex troubles
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
2008 Sep 16
1
Different PCA results under Windows and Linux
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
2011 Mar 05
2
please help ! label selected data points in huge number of data points potentially as high as 50, 000 !
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
2005 Apr 12
1
factors in multinom function (nnet)
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),
2009 Aug 18
15
dlm stress test hangs OCFS2
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
2002 Jun 09
1
Rsync of filesystems with sym-links to other hosts
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
2011 Sep 20
1
Data
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))
2013 Mar 19
1
How to automate this model selection algorithm?
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