On 03/06/2012 03:49 PM, R_beginner_starter wrote:
My input file:
...
http://r.789695.n4.nabble.com/file/n4448762/ScreenHunter_01_Mar._06_12.42.jpg
Hi R_beginner_starter,
If you don't really need the 3D effect, you can get more or less that
plot with this (br is your data):
library(plotrix)
many thanks for the response.
Rstudio seems to be a very clever robust tool and i suspect that if R
can
do it (ie send text output to its own window) then Rstudio would probably
handle it.
Can it be done in 'vanilla' R ?
cheers Bob
Uwe Ligges
On 06.03.2012 10:37, Robert Kinley wrote:
many thanks for the response.
Rstudio seems to be a very clever robust tool and i suspect that if R
can
do it (ie send text output to its own window) then Rstudio would probably
handle it.
Can it be done in 'vanilla' R ?
No, but you can sink()
On 03/06/2012 06:42 PM, Ajay Askoolum wrote:
I expected the row names to be unique but a data frame appears to be able to
hold duplicate row names. This makes me thing that there must be circumstances
when it is necessary. However, I cannot think of any. Please enlighten me.
Hi Ajay,
An
HI, I would like to update my DESeq package version on R-2-14 using
bioclite() and get this message, could somebody help please?
biocLite(DESeq)
BioC_mirror: 'http://www.bioconductor.org'
Using R version 2.14, BiocInstaller version 1.2.1.
Installing package(s) 'DESeq'
Installing package(s)
Jim, have a look here:
http://www.r-bloggers.com/select-operations-on-r-data-frames/. Specifically at
the Duplicate Row Names section towards the end of the page.
From: Jim Lemon j...@bitwrit.com.au
Cc: R General Forum r-help@r-project.org
Sent: Tuesday, 6
On 03/06/2012 02:20 AM, nathalie wrote:
HI, I would like to update my DESeq package version on R-2-14 using
bioclite() and get this message, could somebody help please?
biocLite(DESeq)
BioC_mirror: 'http://www.bioconductor.org'
Using R version 2.14, BiocInstaller version 1.2.1.
Installing
On 06-03-2012, at 11:23, Ajay Askoolum wrote:
Jim, have a look here:
http://www.r-bloggers.com/select-operations-on-r-data-frames/. Specifically
at the Duplicate Row Names section towards the end of the page.
Look carefully at that example.
It is NOT providing rownames for a dataframe.
Currently, the GtkDrawingArea object has no real knowledge that it is being
used as a graphics device. You could do something like: stick the device ID
on da1/da2 as an attribute, and then have a function that does dev.set with
that attribute.
Michael
On Mon, Mar 5, 2012 at 1:14 PM, Mark
Hello,
I am trying to reshape a data.frame in wide format into long format.
Although in the reshape R documentation
they programmer list some examples I am struggling to bring my data.frame
into long and then transform it back into wide format. The data.frame I look
at is:
df - data.frame(ID1
Thanks a lot, Jim.
Do you mind to post out the complete content of your R script suggested?
I got add in your R code in my existing R script.
It shown the following error message:
Error in rbind(br[, 2] + 1, br[, 4]) : object 'br' not found
Thanks for your further reply.
--
View this message in
Hi,
I work with hyperspectral remote sensing data and I try to built a pls
model with this data. I already built the model but if I try to
calculate the RMSEP and R2 with a test data set I get the following
error message:
Error: variable 'subX' was fitted with type nmatrix.501 but type
Dear List,
I managed to solve a part of my problem. Anyway, I still do not
understand how to use the lattice package the right way.
An example:
I do have multiple stations (say, with climatological data) and I want
to plot each station on its on panel. I managed to do this of course
with one
Hi sir,
When is used search option as ?lm, it showed following thing:
?lm
starting httpd help server ... done
Error in shell.exec(url) :
access to 'http://127.0.0.1:16499/library/stats/html/lm.html' denied
how to get rid of it
..plz help
waiting 4 response
Naveen Verma
Hi,
I have a question about understanding PLS. If I use the predict function
of R than it seems to me the function only uses the last latent variable
to model new Y values. But should the function not use all latent
variables to model new Y´s?
Thank you very much for your help.
/Thomas
--
i run like ,but it giving other than what i expect
i also want to capture output print it
e-system(date)
e
[1] 127
e-system(dir)
e
[1] 127
--
View this message in context:
Hi
I have a matrix
a = 1 2 3
1 2 3
1 2 3
1 2 3
1 2 3
1 2 3
b = 5 6 7
8 9 10
I want output as
c = 5 2 3
8 2 3
1 6 3
1 9 3
1 2 7
1 2 10
Please any on help to avoid reduncy of my code
-
Please read the help for system, paying particular attention to the
intern argument.
e - system(date, intern=TRUE)
e
[1] Tue Mar 6 08:33:48 EST 2012
Sarah
On Tue, Mar 6, 2012 at 6:56 AM, sagarnikam123 sagarnikam...@gmail.com wrote:
i run like ,but it giving other than what i expect
i also
Those are error messages from the shell, not real values for what
you think you are doing. Are you looking for the getwd() and
Sys.date() functions instead?
Michael
On Tue, Mar 6, 2012 at 6:56 AM, sagarnikam123 sagarnikam...@gmail.com wrote:
i run like ,but it giving other than what i expect
i
On Mar 6, 2012, at 13:02 , arunkumar wrote:
Hi
I have a matrix
a = 1 2 3
1 2 3
1 2 3
1 2 3
1 2 3
1 2 3
b = 5 6 7
8 9 10
I want output as
c = 5 2 3
8 2 3
1 6 3
1 9 3
1 2 7
1 2 10
library(reshape2)
melt(df, id.vars = c(ID1, ID2, ID3))[, -4]
# To drop an extraneous column (but you should take a look and see
what it is for future reference)
Michael
On Tue, Mar 6, 2012 at 6:17 AM, mails mails00...@gmail.com wrote:
Hello,
I am trying to reshape a data.frame in wide
On 3/5/2012 6:35 AM, Robert Kinley wrote:
platform: R 2.14.2 -Windows XP - Rstudio
When I use cat() or print() in a script, the output text is of course
mixed up with the lines of the script.
I can sink() the output to a textfile and then play the file back at the
end of the
On Tue, Mar 6, 2012 at 4:20 AM, Stefan Luedtke slued...@gfz-potsdam.de wrote:
I managed to solve a part of my problem. Anyway, I still do not
understand how to use the lattice package the right way.
An example:
I do have multiple stations (say, with climatological data) and I want
to plot
On 12-03-05 10:03 PM, eigenvalet wrote:
I get this error before the prompt every time I start R version 2.14.2.
Don't know if it is causing any problems, but it bugs me. Any one have any
ideas? Running R 64 bit on Windows 7 Ultimate SP1.
Presumably you have a startup file containing
Just looking at this, but it looks like ix doesn't exist:
sapply(1:length(inxlist), function(i) if(length(ix[[i]]))
fin1[ix[[i]], tkr
+ 1] - ua[i, tkr])
Trying to sort it out now.
Ben
On Mon, Mar 5, 2012 at 7:48 PM, Rui Barradas rui1...@sapo.pt wrote:
Hello,
Mar 05, 2012; 8:53pm
I think this is what you meant:
z.dates =
c(2007-03-31,2007-06-30,2007-09-30,2007-12-31,2008-03-31,2008-06-30,2008-09-30,2008-12-31)
nms = c(A,B,C,D)
# these are the report dates that are the real days the data was available
rd1 = matrix(c(20070514,20070814,20071115, 20080213,
20080514,
Thomas Möckel thomas.moc...@nateko.lu.se writes:
I work with hyperspectral remote sensing data and I try to built a pls
model with this data. I already built the model but if I try to
calculate the RMSEP and R2 with a test data set I get the following
error message:
Error: variable
Thomas Möckel thomas.moc...@nateko.lu.se writes:
I have a question about understanding PLS. If I use the predict function of R
than it seems to me the function only uses the last latent variable to model
new Y values. But should the function not use all latent variables to model
new Y´s?
It
Unfortunately, your solution is does not scale well. (Tough for you to
test this without my real data.) If ua is my data and rd1 are my report
dates (same as the code below) and I use more columns, it appears that your
solution slows considerably. Remember I have ~11k columns in my real data,
so
Hi list:
I work with a lot of laboratory analytical data and I often have
inconsistent names of files and variables within those files so I wrote
this sourcearg function to facilitate handling file and variable names
as both character and R names. The source of the function is given below
Just using the reshape() function in base R:
df.long = reshape(df, varying=list(names(df)[4:7]), direction=long)
This also gives two extra columns (time and id) can can be dropped.
Andy
-Original Message-
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On
On Mar 6, 2012, at 10:24 AM, Lucas wrote:
PERFECT
Thank you so much.
If by 'execute' you meant translate a name for a function then you
should also look at do.call.
--
David.
Lucas.
2012/3/6 R. Michael Weylandt michael.weyla...@gmail.com
eval(parse(text = STRING))
but you often
Hello to everyone.
In MATLAB you can find a function that executes a text string.
I have not find it in R.
Does anyone knows a similiar function in R, i think a valuable tool.
Thank you.
Lucas.
[[alternative HTML version deleted]]
__
The basic question is what do you want the function to do. It seems
convoluted and I am not sure exactly what you are giving it as input
and what you expect as output. Not exactly sure of what the 'sink' and
'source' are supposed to do. So tell me what you want to do, not how
you want to do it.
Hi,
Is there a simple way of doing the following in R?
a - data.frame(name = c(rep(A, 3), rep(B, 5), rep(C, 10)))
a
name
1 A
2 A
3 A
4 B
5 B
6 B
7 B
8 B
9 C
10C
11C
12C
13C
14C
15C
16C
17C
18C
Do some fancy R code
S+ does allow duplicates in the row names if you ask for them
to be allowed (with data.frame(..., dup.row.names=TRUE) or
attr(df, dup.row.names) - TRUE).
They were allowed because bootstrappers were sampling the
rows of data.frames with replacement and it turned out that
a fair bit of time was
Here's one possible approach. It assumes that a$name is a factor, as
it is in your example, but does not require that each sequence has a
unique value (see second example).
a - data.frame(name = c(rep(A, 3), rep(B, 5), rep(C, 10)))
data.frame(name=a,
eval(parse(text = STRING))
but you often don't need to do this. A somewhat canonical quote in the R-world:
If the answer is parse() you should usually rethink the question.
-- Thomas Lumley, R-help (February 2005)
Michael
On Tue, Mar 6, 2012 at 10:06 AM, Lucas lpchaparro...@gmail.com wrote:
PERFECT
Thank you so much.
Lucas.
2012/3/6 R. Michael Weylandt michael.weyla...@gmail.com
eval(parse(text = STRING))
but you often don't need to do this. A somewhat canonical quote in the
R-world:
If the answer is parse() you should usually rethink the question.
-- Thomas Lumley, R-help
Thanks so much it worked very well
-
Thanks in Advance
Arun
--
View this message in context:
http://r.789695.n4.nabble.com/Help-in-matrix-tp4449622p4450560.html
Sent from the R help mailing list archive at Nabble.com.
__
Hey,
for a statistical analysis I have to work with some micro dates in SPSS.
I only have txt-files though (any further details concerning the format or
something like that, I don´t have. I can only see a lot of numbers when I
open it). On the same page I found another Download especially for
Dear all,
I'm trying to construct confidence intervals for a LOWESS estimation (by not
using bootstrapping).
I have checked previous posts and other material online and I understand that
the main procedure is:
my.count- seq(...)
fit- loess (y ~ x, data=z)
pred- pred(fit, my.count,
Hi!
Does anybody know about an R function for importing paired
name/value data, where each pair is separated by a fixed
string (e.g., a semi-colon), and within each pair a
different delimeter is used (e.g., an equals sign).
Example:
Thanks a lot. I have changed the calculation method by using optimism defined
by Efron. The results from using boot and rms packages are quite close now.
Jim
--- On Tue, 6/3/12, Frank Harrell f.harr...@vanderbilt.edu wrote:
From: Frank Harrell f.harr...@vanderbilt.edu
Subject: Re: [R] 632
Dear list,
A follow up on this thread - the solution I ended up with does indeed
involve parse().
However, it does the job, and it is more efficient than what I came up with
in my other attempts.
suspicious.vowels - function(data,factors,vowelcolumn,f1,f2) {
for(currfac in
I'd just sort them first, if you want consecutive numbers for each
value. (It would have been nice to specify that in the original
question.)
If you need for some reason to put them back in order, you can always
cbind(1:nrow(x), x) before sorting, to get an index to sort on after
you're done.
As
Hey Michael,
thanks for your help. I think the reshape2 library works much better than
the normal reshape function.
However, I still cannot retransform my data.frame into wide format once used
melt function to transform
it into long format. How do I get it back to wide format?
The documentation,
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
project.org] On Behalf Of Aims
Sent: Tuesday, March 06, 2012 9:00 AM
To: r-help@r-project.org
Subject: [R] TXT-file with R to SPS-file
Hey,
for a statistical analysis I have to work with some micro
Hi Sarah,
thanks a lot for this peace of code.
Is it possible to give the second sequence of B in your second example
(data frame b)
the numbers 6, 7 and 8 instead of 1, 2 and 3 again? In my real data I have
more columns than
only those two and sometimes the are sorted differently so that
Type at your console:
library(lattice)
example(xyplot.zoo)
I think the second and third do what you were asking. Otherwise,
please describe in more detail what you mean about add other
variables to each panel -- how do you want these displayed and what
sort of variables are they?
On Tue, Mar 6,
I feel this is a very easy thing but I've never done it before and it is
getting frustrating.
I have a big data.frame (1445846 rows, 15 col)
that looks like this:
V1 V2 V3 V4 V5
1home sisterbrother chair 0
2cat
You'll need to use the cast() function but I can't say more unless you
can provide more specifics. I believe Hadley's website has more
documentation than the package: http://had.co.nz/reshape/
Michael
On Tue, Mar 6, 2012 at 9:34 AM, mails mails00...@gmail.com wrote:
Hey Michael,
thanks for
Good afternoon,
MilanoR.net is proud to announce the first MilanoR meeting on Tuesday 8th
May 2012 at 18:00.
Venue:
Fiori Oscuri Bistrot
Via Fiori Oscuri 3
IT-20122 Milano
www.fiorioscuri.it
Agenda:
Introduction
Andrea Spano'
Quantdie Srl
Past, present and future of R
Stefano Iacus,
Hi All,
I am writing a function that reads a file in
myfile = file('myfile.raw', 'rb')
.
.
.
.
.
close(myfile)
No matter what, I get the warning
Warning message:
closing unused connection 3 (myfile.raw)
Since the whole thing is in a function, I'd like to avoid unecessary noise for
the user,
A reproducible example would be helpful, but lacking that
here's some untested code. If your data frame has NA values,
those will also need to be dealt with.
apply(mydata, 1, function(x)sum(x != 0))
Sarah
On Tue, Mar 6, 2012 at 12:03 PM, mari681 marianna.bologn...@gmail.com wrote:
I feel this
Or perhaps faster but less general:
rowSums(mydata != 0)
Michael
On Tue, Mar 6, 2012 at 12:56 PM, Sarah Goslee sarah.gos...@gmail.com wrote:
A reproducible example would be helpful, but lacking that
here's some untested code. If your data frame has NA values,
those will also need to be dealt
Hello,
Just looking at this, but it looks like ix doesn't exist:
sapply(1:length(inxlist), function(i) if(length(ix[[i]]))
fin1[ix[[i]], tkr + 1] - ua[i, tkr])
Trying to sort it out now.
Right, sorry.
I've changed the name from 'ix' to 'inxlist' to make it more readable just
before
i have file as in txt format
@ATTRIBUTE f1996 REAL
@ATTRIBUTE f1997 REAL
@ATTRIBUTE f1998 REAL
@ATTRIBUTE f1999 REAL
@ATTRIBUTE f2000 REAL
@ATTRIBUTE class {-1,1}
@DATA
2.080750,1.099070,0.927763,1.029080,-0.130763,1.265460,-0.43
1.109460,0.786453,0.445560,-0.146323,-0.996316,0.555759,0.29
In following a thread on this mailing list, I encounter an apparent issue
with **Edit | Paste commands only** in the Windows R-GUI.
Reproducible steps:
Go to (using IE?):
http://www.r-bloggers.com/select-operations-on-r-data-frames/
Find the section entitled 'Duplicate row names'
Copy the
The help for warning offers some suggestions.
Sarah
On Tue, Mar 6, 2012 at 12:53 PM, Federico Calboli
f.calb...@imperial.ac.uk wrote:
Hi All,
I am writing a function that reads a file in
myfile = file('myfile.raw', 'rb')
.
.
.
.
.
close(myfile)
No matter what, I get the warning
i used it on windows system, but giving error like
e - system(date, intern=TRUE)
Error in system(date, intern = TRUE) : 'date' not found
--
View this message in context:
http://r.789695.n4.nabble.com/how-to-run-system-command-tp4449597p4449906.html
Sent from the R help mailing list archive at
You could perhaps set the comment.char argument of read.table() to @.
More generally, you can use scan() and process the lines within R.
On Tue, Mar 6, 2012 at 8:41 AM, sagarnikam123 sagarnikam...@gmail.com wrote:
i have file as in txt format
@ATTRIBUTE f1996 REAL
@ATTRIBUTE f1997 REAL
Dear useRs,
I am sure this is a fairly simple problem, but I just cannot get my head around
it.
I have a dataframe which contains several factor variables. I can use table()
to tell me how many different combinations there are of these variables. What I
should like to do is to add a column to
On Tue, Mar 6, 2012 at 8:48 AM, sagarnikam123 sagarnikam...@gmail.com wrote:
i used it on windows system, but giving error like
e - system(date, intern=TRUE)
Error in system(date, intern = TRUE) : 'date' not found
If you type date at a windows command line does it work?
You can't use
I spent some more time to understand the examples .. it helped a bit.
But there is still one issue.
Taking example data (my first on, hope it is not confusing ;-) ) like:
library(latticeExtra)
library(zoo)
# example Data
timeVect1=seq(as.Date(1990-01-01), as.Date(2010-12-31), by=day)
On Mar 6, 2012, at 12:39 PM, Nordlund, Dan (DSHS/RDA) wrote:
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
project.org] On Behalf Of Aims
Sent: Tuesday, March 06, 2012 9:00 AM
To: r-help@r-project.org
Subject: [R] TXT-file with R to SPS-file
Hey,
for
On Mar 6, 2012, at 12:45 PM, R. Michael Weylandt wrote:
You'll need to use the cast() function but I can't say more unless you
can provide more specifics. I believe Hadley's website has more
documentation than the package: http://had.co.nz/reshape/
in reshape2 Hadley did away with cast() and
One possible approach is to use unique() to get the list of distinct
combinations, cbind() an identifying variable to that list, then use
merge() to join it to your existing data frame.
But I'm not seeing how you are getting four unique combinations.
Given your sample data (with the missing comma
On 6 Mar 2012, at 18:01, Sarah Goslee wrote:
The help for warning offers some suggestions.
none that seem to work though.
F
Sarah
On Tue, Mar 6, 2012 at 12:53 PM, Federico Calboli
f.calb...@imperial.ac.uk wrote:
Hi All,
I am writing a function that reads a file in
myfile =
Well, if you can get this to run your version of R is markedly\
different than mine.
#Start of code
x1=c(rep(0:1,6))
x2=c(rep(c(1,1,0,0)6))
Error: unexpected numeric constant in x2=c(rep(c(1,1,0,0)6
x3=c(rep(1,6),rep(0,6))
On Tue, Mar 6, 2012 at 1:23 PM, O'Hanlon, Simon J
On Tue, Mar 6, 2012 at 1:20 PM, Federico Calboli
f.calb...@imperial.ac.uk wrote:
On 6 Mar 2012, at 18:01, Sarah Goslee wrote:
The help for warning offers some suggestions.
none that seem to work though.
Where's your reproducible example of this? Doesn't
seem to work is unanswerable.
F
On 06/03/2012 1:02 PM, Robert Baer wrote:
In following a thread on this mailing list, I encounter an apparent issue
with **Edit | Paste commands only** in the Windows R-GUI.
Reproducible steps:
Go to (using IE?):
http://www.r-bloggers.com/select-operations-on-r-data-frames/
Find the section
On Tue, Mar 6, 2012 at 1:32 PM, O'Hanlon, Simon J
simon.ohan...@imperial.ac.uk wrote:
Ah!
Thanks.
I had already made vector x2 previously and then went and changed it for some
reason, which was why I didn't notice the error (because the subsequent code
was able to run regardless). Sorry
On Tue, Mar 06, 2012 at 09:32:37AM -0800, syrvn wrote:
Hi Sarah,
thanks a lot for this peace of code.
Is it possible to give the second sequence of B in your second example
(data frame b)
the numbers 6, 7 and 8 instead of 1, 2 and 3 again? In my real data I have
more columns than
only
On 06/03/2012 17:53, Federico Calboli wrote:
Hi All,
I am writing a function that reads a file in
myfile = file('myfile.raw', 'rb')
.
.
.
.
.
close(myfile)
No matter what, I get the warning
Warning message:
closing unused connection 3 (myfile.raw)
Since the whole thing is in a function, I'd
As I have seen no response to your post
1. I am in no sense an expert on loess, so caveat emptor.
2. AFAICS, you just need to read the seq() Help file carefully. I do
not know how to tell you how to use seq() more clearly than it already
does, and so will not attempt to do so. You might
Hello everyone!
I'm working with Decision tree and I have doubt about one of the arguments
of plot.rpart function:
When we use uniform=F, the vertical spacing of nodes will be proportional
to the error in the fit.
But, I want to build a scale next my classif tree to show it.
So, how could I
Dear R users,
Given a user-defined cumulative distribution function F, I want to compute
F^{-1}(x). How is that possible with R?
Best Regards,
--
Gildas Mazo
PhD student
MISTIS team at INRIA
Grenoble, France
[[alternative HTML version deleted]]
Rui and Dimitris, thank you much, both approaches work great! (although rle
won me, I had a feeling something like that had to exist)
--
View this message in context:
http://r.789695.n4.nabble.com/split-array-into-groups-by-value-tp4448578p4451048.html
Sent from the R help mailing list archive
Check function uniroot().
I hope it helps.
Best,
Dimitris
On 3/6/2012 7:47 PM, Gildas Mazo wrote:
Dear R users,
Given a user-defined cumulative distribution function F, I want to compute
F^{-1}(x). How is that possible with R?
Best Regards,
--
Dimitris Rizopoulos
Assistant Professor
system(cmd.exe /c date /T, intern=TRUE)
[1] Tue 03/06/2012
Bill Dunlap
Spotfire, TIBCO Software
wdunlap tibco.com
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of sagarnikam123
Sent: Tuesday, March 06, 2012 5:49 AM
To:
I have a data frame in wide format. There are six variables that represent two
factors in long format 3x2, Valence and Temperature:
head(dpts)
File Subj Time Group PainNeg.hot PainNeg.warm SociNeg.hot
SociNeg.warm Positiv.hot Positiv.warm Errors
1 WB101_1_1_dp.txt 1011
I know this is something simple that I cannot do because I do not yet think
in R.
I have a data frame has a variable participation (a factor), and several
other factors.
I want a chisq test (no contingency tables) for participation vs all of the
other factors.
In SPSS I would do:
CROSSTABS
List and Jim
-Original Message-
From: jim holtman [mailto:jholt...@gmail.com]
Sent: Tuesday, March 06, 2012 09:47
To: Robert Sandefur
Cc: r-help@r-project.org
Subject: Re: [R] sourcearg function is there a better way already built into R
The basic question is what do you want the
I wish to fit a dynamical model in R and I am running in a problem that
requires some of your wisdom to solve. For SAS users I am searching for
the equivalent of the retain statement.
For people that want to read complicated explanations to help me:
I have a system of two equations written as
List and Jim:
My solution works to my satisfaction so perhaps I should not have not ask. The
problem was basically this:
1. I receive a CSV file containing assays on various chemical elements. The
names of these fields are inconsistent in terms of spelling, length and case.
2. As an example I
?pmatch
?switch
_may_ be of use...
But if the file calls it Fe and you want to label it Iron, then
you'll need some sort of explicit dictionary.
-- Bert
On Tue, Mar 6, 2012 at 12:57 PM, Robert Sandefur rsande...@cam-llc.com wrote:
List and Jim:
My solution works to my satisfaction so
Hello,
In case anyone is interested in a faster solution for lots of columns. This
solution is slower if you only have a few columns. If anyone has anything
faster, I would be interested in seeing it.
### some mockup data
z.dates =
Hello, I am in search of any command or package that can be used for exclusion
of some variable estimates from one model but not from the other. I am using
MNP package for data analysis and want the result in the form of unequal number
of covariates for each model choice.
Saqlain RAZA
PhD
Dear all,
I know that it is a silly question, but I have a zoo object like this:
valor1 valor2
01/02/08 36.7381 17.0097
01/03/08 36.9296 16.8331
01/04/08 35.6934 16.3539
01/07/08 35.3539 16.3791
01/08/08 36.3811 16.9340
01/09/08 37.9220 17.6655
The
I have a problem with Winbugs.
The code is executed perfectly, the only problem is that the stats option
isn't available (grey color) and i don't know how to make it work.
Maybe if i run the winbugs throug R i will be able to get stats results?
--
View this message in context:
On Mar 6, 2012, at 3:00 PM, David Perlman wrote:
I have a data frame in wide format. There are six variables that
represent two factors in long format 3x2, Valence and Temperature:
head(dpts)
File Subj Time Group PainNeg.hot PainNeg.warm
SociNeg.hot SociNeg.warm Positiv.hot
Try this.
set.seed(123)
a - sample(letters[1:4], 250, replace=TRUE)
b - sample(letters[1:4], 250, replace=TRUE)
c - sample(letters[1:4], 250, replace=TRUE)
d - sample(letters[1:4], 250, replace=TRUE)
e - sample(letters[1:4], 250, replace=TRUE)
df - data.frame(a, b, c, d, e)
result -
On Mar 6, 2012, at 1:47 PM, Gildas Mazo wrote:
Dear R users,
Given a user-defined cumulative distribution function F, I want to
compute F^{-1}(x).
Isn't this quantile()? (But pay attention to the type some say
to use type=1.)
Alternatively: If q = CDF(x) then generate a matrix
On Tue, Mar 6, 2012 at 3:06 PM, Mago84 davidquia...@gmail.com wrote:
Dear all,
I know that it is a silly question, but I have a zoo object like this:
valor1 valor2
01/02/08 36.7381 17.0097
01/03/08 36.9296 16.8331
01/04/08 35.6934 16.3539
01/07/08 35.3539 16.3791
Hi all,
I'm trying to generate a Weibull distribution including four covariates in
the model. Here is the code I used:
T = rweibull(200, shape=1.3,
scale=0.004*exp(-(-2.5*b1+2.5*b2+0.9*x1-1.3*x2)/1.3))
C = rweibull(n, shape=1.5, scale=0.008) #censoring time
time = pmin(T,C) #observed time is
On Mar 6, 2012, at 5:53 PM, FU-WEN LIANG wrote:
Hi all,
I'm trying to generate a Weibull distribution including four
covariates in
the model. Here is the code I used:
T = rweibull(200, shape=1.3,
scale=0.004*exp(-(-2.5*b1+2.5*b2+0.9*x1-1.3*x2)/1.3))
C = rweibull(n, shape=1.5, scale=0.008)
Dear List,
I normally try to solve this kind of problems by myself but I am running out of
time and have
no idea where to start. So I hope you might excuse this rather basic question.
I'm trying to fit a Canonical Negative Binomial model using the ml.nbc function
from the
COUNT package. This
I needed to compute a complicated cross tabulation to show weighted means
and standard deviations and the only method I could get that worked uses a
series of nested for next loops. I know that there must be a better way to
do so, but could use some assistance pointing the way.
Here is my
Thanks for your advise, David.
I did read the help for survreg and using the followings to calculate.
survreg's scale =1/(rweibull shape)
survreg's intercept = log(rweibull scale)
However, the scale in rweibull has been transformed by exp(betaX's). In my
case, the baseline hazard for T is
1 - 100 of 124 matches
Mail list logo
