nsure the above single variable table will
plot correctly? and what can I do to suppress errors on the datasets which
do not have data so the loop continues?*
full code: http://pastebin.com/LB88hpfM
Thank you in advance.
--
all the best,
Lee Mueller
ISA Certified Arborist MI-4148A
Registered
it to plot the barplot based on the table which is
created above. However, R doesn't like the lack of several columns.
Further guidance would be appreciated.
--
all the best,
Lee Mueller
ISA Certified Arborist MI-4148A
Registered Forester #46043
On Wed, May 9, 2012 at 4:30 AM, Jim Holt
m effects with a
>> common variance (to be estimated). In principle this works with any number
>> of arguments to s(...,bs="re").
>>
>> See ?random.effects (and its linked help files) in mgcv for more.
>>
>> There are mechanisms for allowing random smooth
tructed a data frame with all the time
points for all subjects and fed it into the predict.gam() function. If I
only want to get the CI for the fixed effect only for the derivatives, what
should I change here?
Many thanks in advance!
L
On Thu, Apr 6, 2017 at 2:16 PM, Leon Lee wrote:
> Hi, Simon
&
Hello everyone,
I have a question and I need your precious kind help.
I am working on matching two string. However, the length of the two strings
is different. For example, one is "example" (nchar=7), the other one is
"example " (nchar=10). The R considers them as different strings but t
Hi Jeff,
Thank you a lot!
Best,
Yen
-Original Message-
From: Jeff Newmiller [mailto:jdnew...@dcn.davis.ca.us]
Sent: Thursday, May 18, 2017 11:11 PM
To: r-help@r-project.org; Yen Lee
Subject: Re: [R] Question about change the length of a string.
http://stackoverflow.com/questions
Hi Dear Rusers,
I am working on a survey data with the "survey" package. The logistic
regression and multinomial regression would be the main statistic method I
want to use. I found that the svyglm function could be used to conduct the
logistic regression with the complex design but not the mul
Hi! All.
I'm not much familiar with R.
So I tried to find a R function or packages that could work with my problems.
What I wonder is,
Whether there is any R function or package that includes the cluster analysis
considering with the weighted attribute.
I saw several papers that dealt w
Hello forum members,
I am taking this R course which I have to admit that it is beyond my
capability. I am asked to analyze the following data(file attached) according
to the question provided below:
1. The file [data_13-9.txt] contains a data of length 225, which seems to
have some cycle pattern
from the list produced by a call of the form:
data(mtcars)
lmTest = lm(cbind(mpg, hp) ~ disp + wt + carb, data=mtcars)
Best,
--Lee
--
Lee Hachadoorian
Assistant Professor in Geography, Dartmouth College
http://geospatial.commons.gc.cuny.edu
http://freecity.commons.gc.cuny.edu
earch.
Thanks !
All the best,
Ahreum Lee.
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://ww
Hi,
I am looking for an implementation of the fused lasso that allows the
predictor matrix to be an "abstract" linear operator, namely the cumulative
sum (that is, (X.b)_i = sum(b_k, k=1..i)) (due to the size of the problem,
forming the entire matrix is unlikely to be a good approach). Any pointer
Hi, R experts
I am new to R & GAM toolbox and would like to get inputs from you all on my
models. The question I have is as follows:
I have 30 subjects with each subject being scanned from one to three times
in the first year of life. The brain volume from each scan was measured.
The scan time was
dels list, do you mind sending me a link?
Thank you!
Longchuan
On Tue, Mar 28, 2017 at 4:28 PM, David Winsemius
wrote:
>
> > On Mar 28, 2017, at 9:32 AM, Leon Lee wrote:
> >
> > Hi, R experts
> >
> > I am new to R & GAM toolbox and would like to get inputs f
Hi R-listers,
>
> I would like to know how can i extract (ncomps) component no. when the
> RMSEP is lowest for PLS package?
>
> Currently, I only plot it manually, validationplot(pls) and then only feed
> the ncomp to the jack knife command. However, I would like to automate this
> step.
>
> Pleas
Hi Bjorn,
Greetings. I would like to know how can i extract component no. when the
RMSEP is lowest?
>
> Currently, I only plot it manually and then only feed the ncomp to the jack
> knife command. However, I would like to automate this step.
>
> Please let me know. Many thanks.
>
> Rgrds,
>
Hi,
I have issue a jacktest command and it produces me a list of numbers. I
would like to map them to a dataframe so I could save it to a file.
I tried to convert the mvr.obj to vector but when I check the class, it is
still jacktest.
converted<-as.vector(mvr.obj)
Or any otehr way to map this t
mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
--
Sincerely
- Derrick
Derrick Lee, MSc Candidate
Department of Statistics
The University of British Columbia
LSK-314A |
Hello,
I'm running R 2.7.2 in a windows XP environment and I run the
following in an R console:
library(ncdf)
nc <-open.ncdf('c:/file.nc')
aa <- get.var.ncdf(nc,'var1')
This works fine, but 'aa' takes up about 100mb and I want to release
the memory after using it. I try:
rm(aa)
close.ncdf(nc)
I created variables automatically like this way
for(i in 1:5){
nam <- paste("a",i,sep="")
assign(nam,1:i)
}
and then, i want to insert a new data into "a2" variable. so, i did next
sentence
paste("a",2,sep="") <- 4
so, i got this error message
Error in get(paste("a", 2, sep =
I'm new to post in R-help and my native language is not English.
I apologize if my sentence is not fluent to read.
I am doing a simulation study and I need to execute SAS and read a SAS code
in R.
I try the following code but it doesn't work.
system('"c:\\program files\\SAS\\SAS 9.1\\sas.exe
-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Ben Bolker
Sent: Monday, March 01, 2010 12:32 AM
To: r-h...@stat.math.ethz.ch
Subject: Re: [R] Calling SAS from R
Yen Lee ntu.edu.tw> writes:
> I apologize if my sentence is not fluent to read
Thanks for all insightful replies.
I've solved the problem by using the code as following,
system('"C:\\Program Files\\SAS\\SAS 9.1\\sas.exe"
C:/entropy/output7/scale/syntax.sas')
Thanks for you all.
Yen Lee
-Original Message-
Fro
Hi guys,
I have a very simple question.
I'm trying to make multiple columns to a single column.
For example,
*ttx1* is a 46*72 matrix.
so, I tried this.
*d1=ttx1[,1]
d2=ttx1[,2]
...
d72=ttx1[,72]*
now, d1, d2, ...,d72 become a 46*1 matrix.
And then.. I tried..
*dd=rbind(d1, d2, ..., d72)*
I
Hi. I have a dataset of unique coordinates (no missing coordinates). When I did
the following:
lonlat<-mydata[, c("lon, "lat")]
sp1<-corSpatial(1,form=~lon+lat,type="g")
scor<-Initialize(sp1, lonlat , nugget=FALSE)
I got a value for sp1 (2.67). But got an error message for scor--Error in
getC
ent is much too small")
:
the condition has length > 1 and only the first element will be used
3: In if (dd < 100 * .Machine$double.eps) return(from) :
the condition has length > 1 and only the first element will be used
4: In 0L:n : numerical expression has 2 elements: only th
I have two data frames ( x and y -- sample values below). The rows have
HH:MM:SS and columns have names of devices.
I am trying to find a list of 5 least used devices during recorded time period.
When apply function is used to sum on the column, I get the correct answer on
data frame called x
# remove time from column
return(as.data.frame(pivot))
}
On Apr 25, 2010, at 11:07 PM, David Winsemius wrote:
>
> On Apr 25, 2010, at 10:15 PM, robert lee wrote:
>
>> I have two data frames ( x and y -- sample values below). The rows have
>> HH:MM:SS and columns have names
Hi guys,
I need your help!!
My goal is to make a csv file from ncdf file.
This is the code i've used :
> hyo=open.ncdf("C:/CRUTEM3.nc")
> hyo
[1] "file C:/CRUTEM3.nc has 4 dimensions:"
[1] "longitude Size: 72"
[1] "latitude Size: 36"
[1] "unspecified Size: 1"
[1] "t Size: 1916"
[1] "-
Hi guys,
I need your help.
I would like to select a subset from a dataset.
This is the dimension of the dataset.
> dim(data1)
[1] 72 36 1916
so, it's like.. there are 1916 of 72 * 36 matrix. ==> looks like 72 * (
36*1916 )
**
*1)*
And I would like to select the first 72*36 matrix. This is h
Hi guys,
I need your help.
I'm trying to sort the data by the variable "obs".
This is how I tried to sort the data below.
The problem is, I have a variable name "obs"; this is.. a counter variable.
something like _n_ in SAS.
I do not know why it is not working.
I even tried a similar example in UC
Thanks.. but not working...
> sortedx1=ttx1[order(-ttx1$obs),]
Error in ttx1$obs : $ operator is invalid for atomic vectors
On Fri, Oct 2, 2009 at 8:47 PM, P Ehlers wrote:
> Try
>
> sortedx1=ttx1[order(-ttx1$obs),]
>
> (and ask yourself where obs lives)
>
> -Peter
ttx1,"c:/i.csv")
> m=read.table("c:/i.csv", header=T, sep=",")
> attach(m)
> sortedx1=m[order(-obs),]
Do you know what I should do?
Thanks.
-Hyo
On Fri, Oct 2, 2009 at 8:50 PM, jim holtman wrote:
> It is clear that 'obs' does not exist.
>
> O
if anyone in the community has used
this library before, and if so how?
Thankyou in advance,
Lee Kelvin
--
School of Physics & Astronomy
University of St Andrews
North Haugh
St Andrews KY16 9SS
United Kingdom
Phone (+44) [0]1334461668
Email l...@st-andrews.ac.uk
(i in 1:4217){b[i,]=colnames(max(cor_expN[i,]))}
And got the following error:
Error in b[i, ] = colnames(max(cor_expN[i, ])) : number of items to replace
is not a multiple of replacement length
Any thoughts?
Lee
[[alternative HTML version deleted]]
___
Thanks! guys for the help.
cheers!
Lee
On Wed, Oct 28, 2009 at 7:03 PM, Tony Plate wrote:
> Here's a simple example that might help you get what you want:
>
> set.seed(1)
>> x <- matrix(rnorm(30), ncol=3, dimnames=list(NULL, letters[1:3]))
>> (xc <- co
---
- - -- ---
- - -- ---
I am running out of ideas. So I decided to post this on mailing list. Please
Help!
Best
Lee
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r
e matrix M accordingly.
Any thoughts are welcome!
Best
Lee
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.h
files.
Any suggestions, thanks,
Wenjie Lee
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and pr
Please read this:
On Fri, Dec 11, 2009 at 6:28 PM, Wenjie Lee wrote:
> Hi R Experts,
>
> I'm aware of pdf(), jpeg(),... functions. But,
>
> 1. Is it also possible to export graphs directly to word or RTF? I use to
> copy and paste graphs but resolutions are not so gr
ut is
> directed to Word including figures and tables. It is still a beta-version
> but works fine.
>
> Frank Bloos
>
> >>> Wenjie Lee 12.12.2009 00:28 >>>
>
> Hi R Experts,
>
> I'm aware of pdf(), jpeg(),... functions. But,
>
> 1. Is it
entre for Plants and the Environment
> University of Western Sydney
> Hawkesbury Campus
> Richmond NSW 2753
>
> Dept of Biological Science
> Macquarie University
> North Ryde NSW 2109
> Australia
>
> Mobile: +61 (0)422 096908
> www.remkoduursma.com
>
>
>
> On
Hi R Users,
I'm using following sample code to save png plots.
png(file="sample.png", width=8, height=6)
Then I copy or use "insert pictures" function to get the image into MS Word.
After copying, if I try to stretch or re-size the image it starts getting
distorted. Is there a way to standardize
images into word (instead of latex).
Thanks,
Wenjie
On Sat, Dec 19, 2009 at 1:49 PM, Duncan Murdoch wrote:
> On 19/12/2009 1:28 PM, Wenjie Lee wrote:
>
>> Hi R Users,
>>
>> I'm using following sample code to save png plots.
>>
>> png(file="sample.png
.off()
>
>
> On Sat, Dec 19, 2009 at 1:28 PM, Wenjie Lee
> wrote:
> > Hi R Users,
> >
> > I'm using following sample code to save png plots.
> >
> > png(file="sample.png", width=8, height=6)
> >
> > Then I copy or use "
Hi,
I'm running R 2.7.2 on windows XP. I'd like to find the maximum of a
3-d array over it's third index to create a 2-d array. For example:
> x <- array(c(1,2,3,10,11,12,3:8),c(2,3,2))
> x
, , 1
[,1] [,2] [,3]
[1,]13 11
[2,]2 10 12
, , 2
[,1] [,2] [,3]
[1,]3
Hi,
I have a question regarding the R command options under Linux
Is there any options that I can save history only?
The only save option is --save, which saves both data and history.
I don't want to save temporary data, because sometimes the it takes much
time if the temp data is too big.
Howeve
I put savehistory() in my profile file, and it works nicely.
Thank you very much, jdeisenberg.
alex
jdeisenberg wrote:
>
>
> alex lee wrote:
>>
>> Is there any options that I can save history only?
>> The only save option is --save, which saves both data and hi
ot' function from library boot
but I am not able to figure out how? I am really stuck. Please help!
Best
Lee
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE
Hi, all
I have a question on mgcv and ns. Now I want to compare the results from
glm, gam and ns. Take a simple model y~x for example.
glm1 = glm(y~x, data=data1)
gam1 = gam(y~s(x), data=data1)
ns1 = glm(y~ns(x),data=data1)
In order to confirm the result from glm1 is consistent to those from gam
Hi, all
I have a question on mgcv and ns. Now I want to compare the results from
glm, gam and ns. Take a simple model y~x for example.
glm1 = glm(y~x, data=data1)
gam1 = gam(y~s(x), data=data1)
ns1 = glm(y~ns(x),data=data1)
In order to confirm the result from glm1 is consistent to those from gam
Hello,
Can you tell me what R function to use to do a two-sample chi-squared
test? I want to see if two distributions are significantly different
from each other, and I don't specify the theoretical distribution of
either. For example, I have the following fake count data:
x <- sample(1:10,50,r
Dear Sir/Madam,
May I seek for your kind help for the use of R software for conducting
two-proportional z-test which is test for non-inferiority? I have found one
saying using function (TwoSampleProportion.NIS(alpha, beta, p1, p2, k,
delta, margin) but the software said there is not such function.
Dear Sir/Madam,
May I seek for your kind help for the use of R software for conducting
two-proportional z-test which is test for non-inferiority? What function
will the R software provide for doing such test with 95%CI and P-value
calculated. Thank you very much!
With kindest regards,
Simon
packages("paran")
library(paran)
IS<- read.table ("C:/Users/Deborah Lee/Documents/R/IS.csv")
paran(IS, cfa=TRUE, graph=TRUE, color=TRUE, col=c("black", "red", "blue"))
When I run the code above, I get the error message below.
in cor(x) : 'x
Thank you very much. I appreciate your help!
Deborah D. Lee, PhD
Associate Director of Student Affairs Research and Assessment
The Pennsylvania State University
105 White Building
University Park, PA 16802
(814) 863-9609
-Original Message-
From: PIKAL Petr
Sent: Wednesday, August 26
Hello, I have a problem... I tried using a function and I get the error:
"ERROR: missing value where TRUE/FALSE needed"
Here is my code:
a<-mctp(Ecoli~Fecha, data=Datos, type = "Tukey", alternative = "two.sided",
asy.method = "mult.t", plot.simci = TRUE)summary(a)
Ecoli, Fecha and Datos are well
Hi,
The function locator() allows me to add points by mouse on function plot()
I have draft. But the points can not be edit once I have made the click.
How can I remove or just move the points I made on the draft?
I use the locator function like this:
> plot(1,1)
> locator(type="o")
Thanks.
ding on the complexity of your graphic, you could also
> just re-draw
> the entire thing every time you want to make a change.
>
> Daniel
>
>
> On Fri, Aug 2, 2013 at 10:17 AM, David Winsemius
> wrote:
>
> >
> > On Aug 2, 2013, at 2:26 AM, Klot Lee wrote:
&g
R code) to do what you need.
>
>
> On Fri, Aug 2, 2013 at 3:26 AM, Klot Lee wrote:
>
>> Hi,
>> The function locator() allows me to add points by mouse on function plot()
>> I have draft. But the points can not be edit once I have made the click.
>> How can I r
hi,
when I am doing factor analysis, there is continous and Non continuous
variables(classified variables). Does "psych" package handle this? Is
there any difference between continous variables and non ones or mixed ones
when doing factor analysis? If there is, which packages should I use then?
M
Hi!
I am just learning the 'svyglm' package and have run into an error for which I
have not found a solution. My data have been collected from a stratified random
survey. Here is the code:
NB3<-glm.nb(COLNUM~Year+Depth+MESH+offset(LogEffort),data=data)
dstrat <- svydesign(id=~1, strata=~STRATA
Thanks! That worked. However, what error distribution is assumed? I was hoping
to use a negative binomial.
-Original Message-
From: Milan Bouchet-Valat [mailto:nalimi...@club.fr]
Sent: Friday, April 11, 2014 11:37 AM
To: Lee, Laura
Cc: r-help@r-project.org
Subject: Re: [R] weights error
hi,
I have got some data and draw a line graph with ggplot2 geom_line. The
output seems too sharp to me, is there any way to smooth it?
I mean smooth the line, NOT ADDING a smooth line. Or there is other
functions in ggplot2 can help me out with a smooth line?
Thank you.
--
Stand Alone Complex
hi,
I'v got some data attached as apx.csv. It shows the relationship about a
kind of enzyme and the temperature.
"yhf" and "xbt" are 2 kind of cowpea. "apx" is the kind of enzyme and "low"
means the lowest temperature intraday. "date" means 1 to 121 days.
What I want to show with it is how is "a
rocess it so that I can read it into R?
Best,
--Lee
R version 2.14.1 (2011-12-22)
Platform: x86_64-pc-linux-gnu (64-bit)
Linux Mint 13
--
Lee Hachadoorian
Assistant Professor in Geography, Dartmouth College
http://freecity.commons.gc.cuny.edu
__
R-help@r
On 11/08/2012 02:51 AM, Prof Brian Ripley wrote:
On 08/11/2012 07:11, Lee Hachadoorian wrote:
I have a large (105MB) data file, tab-delimited with a header. There are
some odd characters at the beginning of the file that are preventing it
from being read by R.
That is a BOM make in UCS-2
Dear all,
I am using the zeroinfl() function from the pscl package to develop a
zero-inflated Poisson GLM. I would like to calculate the standard errors of
predicted values. I've tried code posted in a previous discussion on this topic
(https://stat.ethz.ch/pipermail/r-help/2008-December/182806
o I don't understand why the package is not available for the version of R I
am running. Anyone has any idea?
---
Shirley Lee
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinf
Hi, I am a student in trouble with R.
please help me!
I have a 16 huge data(.csv).
data1<-read.csv(D://data1.csv, header=FALSE)
.
.
.
data16<-read.csv(D://data16.csv, header=FALSE)
and then, I tried to seperated data w/ column.
a1<-cbind(data1[1:36568,1])
b1<-cbind(data1[1:36568,2])
.
.
.
ac16<
Hi all R-users,
I'm trying to using apply function to input a range of values into a
function I wrote.
I wrote a function with 4 information needed. I would like to make 2 of
them fixed and the other 2 random (but with specified values).
I would like to replicate the function 1 times.
I wa
Hello, everyone,
There's a problem about zero in R and I really need your help.
I have a vector shown as x=c(0.1819711,0.4811463,0.1935151,0.1433675),
The sum of this vector is shown as 1 in R, but when I type 1-sum(x), the
value is not zero, but -2.220446e-16.
I can accept that this value i
Hello, everyone,
Thank you for all your kindness.
I've solved the problem through your help with the function all.equal.
Thank you very much!
Yen
-Original Message-
From: Gabor Grothendieck [mailto:ggrothendi...@gmail.com]
Sent: Wednesday, June 16, 2010 6:40 AM
To: Yen Lee
Cc: r
Hello,
I have two series, one with stream stage measurements every 5 minutes, and
the other with barometric pressure measurements every hour. I want to
subtract each barometric pressure measurement from the 12 stage measurements
closest in time to it (6 stage measurements on either side of the ho
Liviu,
Try
sh -c 'R_DEFAULT_PACKAGES="$R_DEFAULT_PACKAGES Rcmdr" R "$@"'
This is pulled directly from the launcher that Ubuntu creates when it
installs Rcmdr from the repository.
--Lee
On 10/23/2010 11:37 AM, Liviu Andronic wrote:
> On Sat, Oct 23, 2010
l all be restored. Of course, if you're not on Linux, you
might not want to go to the trouble of running a VM. Criteria 1 Y 2 Y 3 Y
4 Y
Hope this helps,
--Lee
--
Lee Hachadoorian
PhD Student, Geography
Program in Earth & Environmental Sciences
CUNY Graduate Center
ill a room when its empty? Does the room,
> the thing itself have purpose? Or do we, what's the word... imbue it."
> - Jubal Early, Firefly
>
>
> From: Aaron Lee To: r-help@r-project.org Date:
> 10/27/2010
> 11:06 AM Subject: [R] Generate random percentages a
questions about specific packages, I would like to
know I'm not trying to do something impossible.
Thanks,
--Lee
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/po
Thanks very much for the information. I also see that
http://cran.r-project.org/bin/windows/contrib/2.12/ReadMe has
information on some of the packages I'm interested in.
Thanks for the help,
--Lee
On Mon, Nov 1, 2010 at 1:48 PM, Duncan Murdoch wrote:
> On 01/11/2010 1:40 PM, Lee Hach
Dear all,
I've written a function and repeated it for 5000 times with loops with
different value, and the messages returned are the output I set and 15
warnings.
I would like to trace the warnings by stopping the loop when warning came
out.
Does anyone know how to make it?
Thanks a lot fo
On Mon, 15 Nov 2010 13:28:40 -0500, Anthony Damico wrote:
> Do you have any other ideas as to how I might diagnose what's going on
> here? Or, alternatively, is there some workaround that would get this giant
> CSV into a database? If you think there's a reasonable way to use the
> IMPORT comman
the upper left and you will see a list of all
objects (the second part of your question), including variables, data
frame, and open packages. Then if you right click on the name of a data
frame and select Edit from the context menu, you will see a
spreadsheet-style editor.
Best,
Lee Hachadoorian
PhD
ebody help me with this? Thanks!
Hana Lee
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and
Hello,
I have a list of data, such that:
[[1]]
[1] 0.00 0.00 0.03 0.01 0.00 0.00 0.00 0.00 0.01 0.01 0.00 0.00 0.01 0.00
0.00 0.03 0.01 0.00 0.01 0.00 0.03 0.16 0.14 0.02 0.17 0.01 0.01 0.00 0.00
0.03 0.00 0.00 0.00 0.00 0.00 0.01 0.00 0.00 0.00 0.00 0.00
[42] 0.00 0.00 0.00 0.00 0.01 0.00 0.00
Hi, all
I would like to convert xls files to xlsx files with R commands in R console
instead of saving xls files as xlsx files after opening xls files.
Please show me how.
Thanks.
Wonjae
--
View this message in context:
http://r.789695.n4.nabble.com/Conversion-to-xlsx-file-tp3487118p3487118.h
Hello,
I am doing an Engle Granger test on the residuals of two I(1) processes.
I would like to get the MacKinnon (1996) critical value, say at 10%. I
have 273 observations with 5 integrated explanatory variables , so that
k=4. Could someone help me with the procedure in R?
Thank you
Hello,
I'm running version R x64 v2.12.2 on a 64bit windows 7 PC. I have two data
vectors, x and y, and try to run archmCopulaFit. Most of the copulas produce
errors. Can you tell me what the errors mean and if possible, how I can set
archmCopulaFit options to make them run? I see in the do
Hello,
I'm running version R x64 v2.12.2 on a 64bit windows 7 PC. I'm trying to read
a text file using read.table where the values have a format like "1,234,567".
What I want is "1234567". Is there a quick way to strip out the commas? I can
use strsplit and paste, but the file is quite larg
Thanks, Peter and to Joshua Wiley and Peter Alspach for similar comments. Gsub
did the trick.
-Original Message-
From: Peter Langfelder [mailto:peter.langfel...@gmail.com]
Sent: Monday, June 13, 2011 4:35 PM
To: Lee, Eric
Cc: r-help@R-project.org
Subject: Re: [R] remove commas in a
im trying to write a for loop so that every leap year, the number of days
becomes to 366 instead of 365. could someone help me out?
and also, this set of data has 99.99s I set all the 99.99 ==NA.
however, when im doing the order function to find the max value of that
year, it still reads 99.99 as t
Hi, I'm new to R. I'm trying to plot my data into a circle.
my data sort of looks like 12,12,4,5,6,5,11,10,3,9,9,9,12,12,2 total of 15
numbers.
I'm trying to add all the same numbers, such that, there are 4 of 12s,1 of
11, 1 of 10, 3 of 9s, and such
so the circle plot would have 4 parts of 12,
x<-circular(dayofmax[,2]*2*pi/365)
res25 <- density(x, bw=25, control.circular=list(units="degrees"))
circularp(res25$x)
plot(res25, col=4, points.plot=TRUE, xlim=c(-1.5,1.5))
res50 <- density(x, bw=50, adjust=2)
lines(res50, col=2)
lines(res50, col=3, shrink=0.9)
this is my code, basically, it
Could you please take my email off the help center? Thank you very much.
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project
Hey List,
I did a multiple regression and my final model looks as follows:
model9<-lm(calP ~ nsP + I(st^2) + distPr + I(distPr^2))
Now I tried to predict the values for calP from this model using the
following function:
xv<-seq(0,89,by=1)
yv<-predict(model9,list(distPr=xv,st=xv,nsP=xv))
The pr
data
and I thought with the predict function the programm would calculate a
value from the model function for every value of calP... ?
2011/3/20 Ista Zahn :
> Hi Anna,
>
> On Sun, Mar 20, 2011 at 2:54 PM, Anna Lee wrote:
>> Hey List,
>>
>> I did a multiple regressio
to newd
> cbind(newd, predict(mod, newdata = newd)) # predictions at new
> data points
>
> # To just get predictions at the observed points, all you need is
> predict(mod)
>
> HTH,
> Dennis
>
> On Sun, Mar 20, 2011 at 11:54 AM, Anna Lee wrote:
>>
t; # Append predictions to newd
> cbind(newd, predict(mod, newdata = newd)) # predictions at new
> data points
>
> # To just get predictions at the observed points, all you need is
> predict(mod)
>
> HTH,
> Dennis
>
> On Sun, Mar 20, 2011 at 11:54 AM, An
Dear Marc,
Your answer on post
https://stat.ethz.ch/pipermail/r-help/2005-March/067634.html
has a broken link.
Thank you,
W
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.
Dear List!
I want to compare medians of non normal distributed data. Is it
possible and usefull to calculate 95% confidence intervals for
medians? And if so - how can this be achieved in R?
Thanks a lot!
Anna
--
Der Inhalt dieser E-Mail ist vertraulich. Sollte Ihnen die E-Mail
irrtümlich zu
Dear List!
I have a unverse (basic population) which is not normally distributed.
Now from this universe I take some subsets. Each subset is normally
distributed within itself. I now want to compare the subsets and see
if they differ significantly. So what is my assumption - normal
distributed dat
1 - 100 of 341 matches
Mail list logo