I am new to both sem and lavaan package ...
I dint exactly get the difference between sem from sem package and sem from
lavaan package... ,
--
View this message in context:
http://r.789695.n4.nabble.com/Structural-equation-modeling-in-R-lavaan-sem-tp3409642p3997527.html
Sent from the R help ma
On Nov 7, 2011, at 12:09 AM, Jeffrey Joh wrote:
I have a table that looks like this:
structure(list(speed = c(3,9,14,8,7,6), result = c(0.697, 0.011,
0.015, 0.012, 0.018, 0.019), house = c(1,
1, 1, 1, 1, 1), date = c(719, 1027, 1027, 1027, 1030, 1030),
id = c("1000", "1",
"10001"
On Nov 7, 2011, at 12:10 AM, David Winsemius wrote:
On Nov 6, 2011, at 10:08 PM, Lanna Jin wrote:
Hi guys,
I have a matrix with values varying from approximately -0.7 to 0.33
that I
want to create a heatmap/levelplot with.
When I execute the levelplot function for my matrix, I end up
Hi
I'm getting the intercepts of the Random effects as 0. Please help me to
understand why this is coming Zero
This is my R code
Data<- read.csv("C:/FE and RE.csv")
Formula="Y~X2+X3+X4 + (1|State) + (0+X5|State)"
fit=lmer(formula=Formula,data=Data)
ranef(fit).
My sample Data
State Year Y X2 X
On Nov 6, 2011, at 10:08 PM, Lanna Jin wrote:
Hi guys,
I have a matrix with values varying from approximately -0.7 to 0.33
that I
want to create a heatmap/levelplot with.
When I execute the levelplot function for my matrix, I end up
getting colors
that are adjusted to the max and min rat
I have a table that looks like this:
structure(list(speed = c(3,9,14,8,7,6), result = c(0.697, 0.011, 0.015, 0.012,
0.018, 0.019), house = c(1,
1, 1, 1, 1, 1), date = c(719, 1027, 1027, 1027, 1030, 1030),
id = c("1000", "1",
"10001", "10002", "10003", "10004")), .Names = c("speed"
Hi,
Try specifying explicit break points together with their corresponding
colors using at and col.regions,
levelplot(m, at= unique(c(seq(-2, 0, length=100), seq(0, 10,
length=100))), col.regions = colorRampPalette(c("blue", "white",
"red"))(1e3))
HTH,
baptiste
On 7 November 2011 16:08, Lanna
Hi guys,
I have a matrix with values varying from approximately -0.7 to 0.33 that I
want to create a heatmap/levelplot with.
When I execute the levelplot function for my matrix, I end up getting colors
that are adjusted to the max and min rather than around 0. In other words,
ideally I would like
In fact, the same thing happens (window freezes) when a long computation
generated by the tcltk window is stopped through the R console.
library(tcltk)
library(svMisc)
win<-tktoplevel()
n<-1
ff<-function(){
for(i in 1:n){progress(i,n)}
}
button<-tkbutton(win,text="test",co
Hello,
I can't seem to figure out how to generate partial dependence plots
for random forest models generated with the 'party' package. Is there
a function for this that I just haven't found yet?
Thanks
-Philip Dilts
__
R-help@r-project.org mailing li
Mucha gracias!! as.vector worked like a charm and, in this case,
produced the same results as c():
c(pdf)*v
as.vector(pdf)*v
At 07:02 PM 11/6/2011, R. Michael Weylandt wrote:
>It looks like pdf is not a "scalar" (that term actually has no
>meaning in R but I know what you mean) but is
Hello useRs
Using the following code:
library(tcltk)
win<-tktoplevel()
ff<-function(){
plot(1:10,1:10)
pol<-locator(1)
print(pol)
}
button<-tkbutton(win,text="test",command=ff)
tkpack(button)
makes the win panel stop responding when the plot is closed befo
Hi
The problem is that BOTH mosaic() and labeling_cells() are calling
seekViewport() to find the right viewport to draw into and for BOTH
plots they are finding the same viewports (on the left side of the
page). The following code solves the problem (for me anyway) by
specifying a different
Hi,
R may not have a special "scalar", but it is common, if informal, in
linear algebra to refer to a 1 x 1 matrix as a scalar. Indeed,
something like:
1:10 * matrix(2)
or
matrix(2) * 1:10
are both valid. Even
matrix(2) %*% 1:10
and
1:10 %*% matrix(2)
work, where the vector seems to be silen
It looks like pdf is not a "scalar" (that term actually has no meaning in R but
I know what you mean) but is rather a 1x1 matrix, as attested by the fact it
has dimensions. If you give dnorm() a matrix it will return one, as it did
here.
Perhaps you should look at the is.matrix() and as.vector
You might find useful tools if you look at Bioconductor as well.
M
On Nov 6, 2011, at 4:58 PM, Minh Bui wrote:
> Hi everyone,
>
> I am new to R-project. I did search through the list for my problem but i
> can't find it. I am sorry if this question has been asked.
>
> I would like to perfo
I would start by reading one or more of the introduction manuals available
here:
http://mirrors.ibiblio.org/pub/mirrors/CRAN/
wizi wrote:
>
> Hi everyone,
>
> I am new to R-project. I did search through the list for my problem but i
> can't find it. I am sorry if this question has been
Hi everyone,
I am new to R-project. I did search through the list for my problem but i
can't find it. I am sorry if this question has been asked.
I would like to perform a correlation analysis between a hiv data and gene
expression.
Basically, i have a file that contains: hiv_name, start_positio
On 11-11-06 3:19 PM, Chee Chen wrote:
Dear All,
I would like to seek for help on this issue:
1. I set par(mfrow=c(2,2)), hoping to plot 4 subgraphs in a whole graph
2. Each subgraph has its own x,y axes and each has x-axis label and y-axis
label
3. moreover, subgraphs in the left column of the
Hello to all,
I am having trouble with intregrating a complicated uni-dimensional function
of the following form
Phi(x-a_1)*Phi(x-a_2)*...*Phi(x-a_{n-1})*phi(x-a_n).
Here n is about 5000, Phi is the cumulative distribution function of
standard normal,
phi is the density function of standard nor
>
> On Nov 5, 2011, at 7:20 PM, ScottDaniel wrote:
>
> > So I have a text file that looks like this:
> > "Label" "X" "Y" "Slice"
> > 1 "Field_1_R3D_D3D_PRJ_w617.tif" 348 506 1
> > 2 "Field_1_R3D_D3D_PRJ_w617.tif" 359 505 1
> > 3 "Field_1_R3D_D3D_PRJ_w617.tif" 35
Hello all,
I have an outstanding request to documents units [1] or other variable-specific
facts in data files. I am thinking in particular of tab-delimited files that
are suitable to be read into a data frame. One suggestion/request is to add a
line of text under the column headings.
My fi
Dear All,
I would like to seek for help on this issue:
1. I set par(mfrow=c(2,2)), hoping to plot 4 subgraphs in a whole graph
2. Each subgraph has its own x,y axes and each has x-axis label and y-axis
label
3. moreover, subgraphs in the left column of the whole graph are all 3D, and
have z axes
Since you did not supply a subset of data, or indicate if you wanted
all the other values transformed in some ways, here is a simple use of
tapply to get your data with respect to 'catch'. If this not what you
wanted, you need to be more clear in your request and also give some
hints as to what yo
There do exist packages for multi-variate integration in R, but sticking to
base functions, what you've described should work but the inner integral will
need to be vectorized before it's passes to the outer integral: Vectorize() can
do this directly, but it won't be particularly fast since it's
Dear list,
I have this dataframe:
> names(events)
[1] "EID""X" "Y" "trip" "tow""catch" "effort" "depth"
[9] "season"
Where some of my unique ID "EID" appears more than once in 162 cases.
> length(events$EID)-length(unique(events$EID))
[1] 162
I would like to combined each
Hello to everyone!
I am working on my final year project about multivariate time series. There
are three variables in the multivariate time series model.
I have a few questions:
1. I used acf and pacf plot and find my variables are nonstationary. But in
adf.test() and pp.test(), the data are s
Hi,
I have a function that I need to do double integration:
\int^T_0 \int^t_0 N(\delta / \sigma \sqrt(u)) (1-N(\delta / \sigma
\sqrt(u))) du dt
where N(x) is a standard normal probability of x.
I start off by writing an inner integral into a function. Meaning
\int^t_0 N(\delta,\sigma \sqrt(u))
Not too difficult.
Rgames> bar<-c(3,5,7)
Rgames> foo<-c(1,3,5,6,8,9,7)
Rgames> match(bar,foo)
[1] 2 3 7 # these are the matching positions
Rgames> foo[-(match(bar,foo))]
[1] 1 6 8 9
Dear list,
I have been struggling for some time now with this code... I have this
vector of unique ID "EI
Well, you could simply use everything from the old library and just apply
update.packages(checkBuilt=TRUE)
in order to get the packages updated for the new release.
Uwe Ligges
On 05.11.2011 19:00, Erich Neuwirth wrote:
Running
rownames(installed.packages())
will tell you the names of all
Hi:
On Sat, Nov 5, 2011 at 11:06 PM, Kaiyin Zhong wrote:
> Thank you Dennis, your tips are really helpful.
> I don't quite understand the lm(y~mouse) part; my intention was -- in
> pseudo code -- lm(y(Enzyme) ~ y(each elem)).
As I said in my first response, I didn't quite understand what you
wer
On 06-Nov-11 14:50:18, David Winsemius wrote:
>
> On Nov 6, 2011, at 9:07 AM, Chee Chen wrote:
>
>> Dear All,
>> I would like to know how to do the following:
>> 1. suppose I have x values from the ordered from 0, 0.5, 1, and
>> would like to label these three points on the x-axis.
>> 2. Howeve
On Nov 6, 2011, at 9:07 AM, Chee Chen wrote:
Dear All,
I would like to know how to do the following:
1. suppose I have x values from the ordered from 0, 0.5, 1, and
would like to label these three points on the x-axis.
2. However, R labels them as 0.0, 0.5, 1.0. But I wan5 them to be
0, .5,
On Nov 6, 2011, at 12:21 AM, R. Michael Weylandt wrote:
There are a few (nasty?) side-effects to c(), one of which is
stripping a matrix of its dimensionality. E.g.,
x <- matrix(1:4, 2)
c(x)
[1] 1 2 3 4
So that's probably what happened to you. R has a somewhat odd feature
of not really consid
Roger,
It's nice to see a reply from the leader in quantile regression. I wonder
if I might ask a somewhat unrelated question. A few recent papers have
developed ways to force quantile regression curves not to cross. Do you
have plans to implement this capability in quantreg?
Thanks very much
On Nov 5, 2011, at 7:20 PM, ScottDaniel wrote:
So I have a text file that looks like this:
"Label" "X" "Y" "Slice"
1 "Field_1_R3D_D3D_PRJ_w617.tif"348 506 1
2 "Field_1_R3D_D3D_PRJ_w617.tif"359 505 1
3 "Field_1_R3D_D3D_PRJ_w617.
Dear list,
I have been struggling for some time now with this code... I have this vector
of unique ID "EID" of length 821 extracted from one of my dataframe (skate). It
looks like this:
> head(skate$EID)
[1] "896-19" "895-8" "899-1" "899-5" "899-8" "895-7"
I would like to remove the compl
Dear All,
I would like to know how to do the following:
1. suppose I have x values from the ordered from 0, 0.5, 1, and would like to
label these three points on the x-axis.
2. However, R labels them as 0.0, 0.5, 1.0. But I wan5 them to be 0, .5, 1,
since the former way uses limited space of a mu
Roger Koenker
rkoen...@illinois.edu
On Nov 5, 2011, at 1:02 PM, Julia Lira wrote:
Dear David,
Indeed rq() accepts a vector fo tau. I used the example given by
Frank to run
fitspl4 <- summary(rq(b1 ~ rcs(x,4), tau=c(a1,a2,a3,a4)))
and it works.
I even can use anova() to test equality
Hello,
The function garchFit in the package fGarch allows for choosing a
conditional distribution, one of which is the t-distribution. The function
allows specification of the shape parameter of the distribution (equal to
the degrees of freedom for the t-distribution), for which the default is set
On Nov 05, 2011 at 11:01pm Francisco Mora Ardila wrote:
> But a problem arised with the predict function: it doesn´t seem to work
> with an object
> from dudi.mix and I dont understand why.
Francisco,
There is no predict() method for dudi.mix() or for any of the dudi objects
in ade4. I don't se
Thank you Dennis, your tips are really helpful.
I don't quite understand the lm(y~mouse) part; my intention was -- in
pseudo code -- lm(y(Enzyme) ~ y(each elem)).
In addition, attach(d) seems necessary before using lm(y~mouse), and
since d$mouse has a length 125, while each elem for each region ha
I sent a post yesterday that I found out why my function didn't work. It's
ok now it works.
Thank you all.
--
View this message in context:
http://r.789695.n4.nabble.com/nested-for-loops-tp3992089p3994996.html
Sent from the R help mailing list archive at Nabble.com.
__
43 matches
Mail list logo
