Old Old Old Old
> > [804] Old Old Old Old Old Old Old Old Old Old Old
> > [815] Old Old Old Old Old Old Old Old Old Old Old
> > [826] Old Old Old Old Old Old Old Old Old Old Old
> > [837] Old Old Old Old Old Old Old Old Old Old Old
1234)
>
> What else could be causing this? Regards
> Liviu
>
Could you be reloading a workspace at start-up that is setting the seed? What
happens if you start R using the --vanilla option?
Dan
Daniel Nordlund
Bothell, WA USA
__
R-help@
d be greater than
a') else dx=b/nint
sum(f[-npts]+f[-1])/2*dx
}
#Call my quadrature
x = seq(0,2000,10)
h = 10*(cos(((2*pi)/2000)*(x-mean(x)))+1)
u = (cos(((2*pi)/2000)*(x-mean(x)))+1)
a = x[1]
b = x[length(x)]
plot(x,-h)
a = x[1];
b = x[length(x)]
#call your q
10)
> h = 10*(cos(((2*pi)/2000)*(x-mean(x)))+1)
> u = (cos(((2*pi)/2000)*(x-mean(x)))+1)
> a = x[1]
> b = x[length(x)]
> plot(x,-h)
> a = x[1];
> b = x[length(x)]
>
> #call your quadrature function. Hint, the answer should be 3.
> f = u*h
> result = myquadrature(f,a,
y one process (e.g. R). And even with
those special incantations, at most you will get no more than about 3.2-3.5 GB.
The other thing to remember is that even if you had more than enough free
space, R requires the free space for an object to be contiguous. So if memory
was fragmented and you
t or doesn't want to use sas.get, it is fairly easy to add a parameter
to the read.ssd function to limit the number of records read from the SAS
dataset being exported. If the OP or anyone else is interested, they can
contact me off-list and I will provide an example (no need to clutt
# Let's say I have 5 objects, object_1, object_2, etc.
for (i in 1:5) {
assign(paste("object_",i, sep=""), i+500)
}
# Now, for whatever reason, I don't know the names of the objects I've
created, but I want to operate on them.
list<-ls(pattern="^obj")
#Is get best?
for (l in list) {
cat("
hin a loop.
Thanks,
Dan
On Wed, Oct 20, 2010 at 4:46 PM, Joshua Wiley wrote:
> Hi Daniel,
>
> get() will work for any object, but cat() may not. cat() should work
> for arrays, but it will be messy even for relatively small ones. For
> example, run:
> cat("Hello&quo
Hi Folks,
I have a pretty simple problem: after building a multivariate linear model,
I need to report my 95% confidence interval for predictions based on future
observations.
I tried doing K-fold cross validation using cv.lm() from the DAAG package,
but it currently only uses the first independen
<< repost because previous attempt was not plain text, sorry! >>
Hi Folks,
I have a pretty simple problem: after building a multivariate linear model,
I need to report my 95% confidence interval for predictions based on future
observations.
I know that one option is to use predict(interval="predi
> > file.exists("C:\\Users\\me\\Documents\\My
PhD\\Modelling\\R\\Scripts\\Functions\\functions 20jan10.r")
[1] TRUE
I suspect it something with windows user rights or file properties but
am not sure what to do.
Any help is appreciated.
Thanks in advance,
Daniel Hayes
_
d read the posing guide referenced at the bottom of every email. It
will tell you how to ask questions and what information you need to provide in
order to get useful help. At a minimum, we need a self-contained reproducible
example that demonstrates the problem you are having. Provide say
ems
are still stored consecutively in memory. Notice that layer one in the stack
was "filled" first, and the first layer was "filled" just like the previous
2-dimensional example. But the items are still physically stored linearly, in
consecutive locations in memory.
Hope
ng
> RODBC. Is that how I need to go?
>
> - Jan
This would be your best resource for databases and R.
https://stat.ethz.ch/mailman/listinfo/r-sig-db
Hope this is helpful,
Dan
Daniel Nordlund
Bothell, WA USA
__
R-help@r-project.org maili
I am correcting you. :-) You are using dim() incorrectly, and not accessing
the array correctly. In all of your examples you should be using dim(3,4,2).
Then you need to specify the indexes of the array element you want to look at.
So, to use your example
ucible
example would be very helpful in helping you. You mention you have several
"data files", but I see no reference to data files in your code. Did you mean
data frames? What is Data_i? A data frame or something else?
You said you n
uot;. The error
message says that the function f cannot find object 'itel'. That is because
itel is defined in the function iterate and is local to that function. You
might want to add a third parameter to function f so you can pass in the value
of itel. (I would avoid using itel f
d values, standard
errors and CV's as on page 65 of the paper, and four residual plots.
Type ?ClarkLDF and ?ClarkCapeCod for help on running Clark's two methods.
Please contact the author, Daniel Murphy at chiefmurph _at_ yahoo.com, with
questions, suggestions, ideas or problems.
If anyon
Hi,
I need to convert general numbers to percentages in R to create a boxplot for
purposes of comparison.
Can you please tell me how I can do that? My data is attached as a csv file.
Thanks.
Daniel
__
R-help@r-project.org mailing list
https
Dear R Help,
I believe the glmer() function in lme4 automatically fits an
unstrucruted covariance matirx for the random effects.
Is that true?If so, do I have an option to somehow ask for a
diagonal structured covariance matrix?
Thank you,
Daniel Jeske
Department of Statistics
University
gonal structure.
Dan
-Original Message-
From: dmba...@gmail.com [mailto:dmba...@gmail.com] On Behalf Of Douglas
Bates
Sent: Tuesday, November 16, 2010 1:01 PM
To: Daniel Jeske
Cc: r-help@r-project.org; R-mixed models mailing list
Subject: Re: [R] Question about GLMER
I am cc:ing the r-sig-mixe
; sessionInfo()
R version 2.12.0 (2010-10-15)
Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit)
locale:
[1] en_US.UTF-8/en_US.UTF-8/C/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] stats graphics grDevices utils datasets methods base
loaded via a namespace (and not attached):
[1
t out. ANy help
would be gratefully received.
Dan
--
**
Daniel Brewer, Ph.D.
Institute of Cancer Research
Molecular Carcinogenesis
Email: daniel.bre...@icr.ac.uk
**
The Institu
instructions on how to change the significance level in R for hypothesis
testing and regression analysis I would appreciate it.
Thanks,
--
Daniel Carawan
Master's of Environmental Studies '12
(804) 263-3116
[[alternative HTML versi
you wanted
irregular intervals. You want to know how many occurrences there are between
midnight and 0900, 0900 and noon, noon and 1500, 1500 and 1800, and 1800 and
midnight.
nx=cut(x,c(0,9,12,15,18,24))
table(nx)
for the second question: ?as.ts()
HTH,
Daniel
--
View this message in context:
http://r.
Hello,
I would be grateful if anybody can help me in finding an R function to
compute weighted Spearman correlation coefficient?
Kind regards,
Daniel
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
a chi squared trend might be
appropriate. I don't have any experience dealing with ordinal variables so I
am at a bit of a loss. What is the most appropriate test? and is it
implemented in R?
Many thanks
******
Daniel Brewer
In
ot;
>
> Thanks for any help
>
> Wade
>
Wade,
As you have discovered, you need to test for NA first, and to do that you need
to use is.na(). Something like this should work
for (i in 1:nrow(demo)) {
if (is.na(demo$Area[i])) Class[i] <- "Sna" else
if (demo$Area[
r ( i in 1:length(y) )
> y[[i]] <- x[[ i + length(x) - length(y) ]]
>
> > y
> [[1]]
> [1] 2
>
> [[2]]
> [1] 3
>
> [[3]]
> [1] 4
>
> [[4]]
> [1] 5
>
> Regards,
> Julian
>
Julian,
You use standard indexing to do that.
y <- x[-1]
uld be an appropriate middle step.
Thanks
Dan
--
**
Daniel Brewer, Ph.D.
Institute of Cancer Research
Molecular Carcinogenesis
Email: daniel.bre...@icr.ac.uk
**
The Institute of Cancer Research: Royal
mat="%d.%m.%y") - 100*365
But that doesn't seem to work out correctly. Any ideas how to do this?
Thanks
Dan
--
******
Daniel Brewer, Ph.D.
Institute of Cancer Research
Molecular Carcinogenesis
Email
On 10/12/2010 4:17 PM, Barry Rowlingson wrote:
> On Fri, Dec 10, 2010 at 3:27 PM, Daniel Brewer
> wrote:
>> Hello,
>>
>> I have some data that has dates in the form 27.02.37. I convert them to
>> a date object as follows:
>> as.Date(data$date,format="
only it's dimensions?
Thanks,
Daniel
[[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
a
with "lsmeans" (this is *not* what I am trying to get here!
). I am just trying to get all possible pairwise comparisons from one
simple command. Already checked the "effects package" but it did not
work with a lme() object.
Thanks in advance, Daniel
___
ment of using read.xport(). Or am I doing something wrong.
Thanks,
Dan
Daniel Nordlund
Bothell, WA USA
__
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/posti
> all(vector[1]==vector)
should be quick.
[[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 pro
ves you the 3 Rows,
> the 3 Columns, and the 3 Values for one of the 20 possible
> choices of 3 cells out of the 6 in the matrix. You could
> subsequently rearrange this to suit your purposes.
>
> Ted.
>
I see lots of good ideas and sample code, but what I haven't seen is
^31 - 1) as the random number. SAS also has
a much better generator, rand(). I would instead generate a series of random
number in one program and export the numbers to the other program.
Hope this is helpful,
Dan
Daniel Nordlund
Bothell, WA
__
lp. You haven't shown
us the code you tried to execute, and we don't know where your data came from.
The error message suggests that your data isn't what you think it is.
So, if this data is in a data frame, then show us the results of
str(your data frame)
Also, show us the a
an download the free
SAS viewer from SAS and convert the file to a csv. There are other options
using various commercial products as well. If you can tell us more about your
situation, someone may be able to give you more options. But unfortunately,
SAS doesn't make it easy to use their nati
/", "ex", sasprog= "C:/Program
Files/SAS/SASFoundation/9.2/sas.exe")
If you don't have access to SAS on your system, this will not work.
Hope this is helpful,
Dan
Daniel Nordlund
Bothell, WA USA
__
R-help@r-project.or
or, and you tell us
that the order in which the variables are returned may change, but give us no
clue as to how they might change. So there is no possible way anyone can give
you any assistance without further information. Give us a sample data frame
with a few columns and few rows and show how
-i]
> }
>
> Is there a way to do this?
>
> thank you
>
> A.Dias
Does this do what you want?
for(i in 1:length(a)) assign(paste('a', i, sep=''), a[-i])
Hope this is helpful,
Dan
Daniel Nordlund
Bothell, WA USA
Hello,
I have a problem programming 2 plot methods for 2 classes.
1. I have 2 classes. Each class has its own plot method. !!! How the
plot,xxx-method.R -files have to look like,
so that the plot method for class 1 und plot methods for class 2 can act
for each belonging class. Prob
My example
aggregate(b ~ a, d, mean)
a b
1 1 1.5
2 2 3.0
It's removing whole rows even if just one of the columns is NA, i.e.:
> d <- data.frame(a=rep(1:2, each=2),
+ b=c(1,2,NA,3),
+ c=c(NA,2,3,NA))
> aggregate(cbind(b,c) ~ a, d, mean)
a b c
1 1 2 2
Dani
Oh wow, that would be it. Not sure how I missed that. Thanks for the tip.
Sent from my iPhone
On Jan 11, 2011, at 18:56, "David Winsemius" wrote:
>
> On Jan 11, 2011, at 5:41 PM, Dickison, Daniel wrote:
>
>> The documentation for `aggregate` makes it sound like
&g
rogram Files\Python27; C:\Program Files\SSH
> Communications Security\SSH Secure Shell;C:\Python26;C:\Python26\Scripts
>
For your system, the path to R should probably be
C:\Program Files\R\R-2.12.1\bin\i386
Hope this is helpful,
Dan
Daniel Nordlund
Bothell, WA USA
___
ubject: Re: [R] R not recognized in command line
>
> @ Daniel Nordlund : I've also tried the dir that you recommended. R is
> still
> not recognized as a command in CMD.exe
>
> @ Uwe Ligges : There are no blanks before or after the semicolon. I added
> them after I copied it
days=Sys.Date()-1:70
price=abs(rnorm(70))
regular=rep(c(0,0,0,0,1,0,1,0,0,1),c(7,7,7,7,7,7,7,7,7,7))
y=data.frame(cbind(days,price,regular))
y is like
days price regular
1 14990 0.16149463 0
2 14989 1.69519358 0
3 14988 1.57821998 0
4 14987 0.47614311 0
5 149
CRAN and see the instructions for downloading and installing the
latest version of R for your version of Ubuntu.
Hope this is helpful,
Dan
Daniel Nordlund
Bothell, WA USA
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo
lt;- 999
# predictor variable
ts <- 1:100
tslength <- length(ts)
var.all <- 1
# pre-allocate vector to collect beta coefficients
b <- numeric(n.k)
# generate your e values all at once
# columns contain values for each simuation trial
e <- matrix(rnorm(n.k*tslength,0,
. 1. How do I concatenate terms that I would
normally quote (like "^") with variables that contain search terms and that
normally would not be quoted? 2. How do I run this over indices of the
variable that contains the search terms?
I greatly appreciate any help,
Daniel
--
View this me
A beauty. Works a charm. Many many thanks.
Daniel
--
View this message in context:
http://r.789695.n4.nabble.com/grep-with-search-terms-defined-by-a-variable-tp2311294p2311307.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r
I have a 100x2 matrix to go into a package. I could type it into the R code
but, more generally, what is the *best* way to read a table from disk into,
say, package 'foo' when it is loaded into memory via library(foo)? I thought
of having .onLoad create an environment and assign a matrix read from
de sample?
>
> Thanks much in advance.
>
> --
> Waverley @ Palo Alto
>
Check out the documentation for regexp
?regexp
Hope this is helpful,
Dan
Daniel Nordlund
Bothell, WA USA
__
R-help@r-project.org mailing list
https://stat.ethz.ch/m
Hi, has there been a solution to this issue? I am encountering the same
problem on a Mac with OSX 10.6.4. The problem persists when I try to install
lme4 from the source (see below), and my R version is up to date according
to R's update check.
Thanks for any help,
D
Hello,
I am using R to train a logistic regression model and save the resulting
model to disk. I am then subsequently reloading these saved objects, and
using predict.glm on them in order to make predictions about single-row data
frames that are generated in real-time from requests arriving at an
Hello,
I am using R to train a logistic regression model and save the resulting
model to disk. I am then subsequently reloading these saved objects, and
using predict.glm on them in order to make predictions about single-row data
frames that are generated in real-time from requests arriving at an
hanks
Dan
--
******
Daniel Brewer, Ph.D.
Institute of Cancer Research
Molecular Carcinogenesis
Email: daniel.bre...@icr.ac.uk
**
The Institute of Cancer Research: Royal Cancer Hospital, a charitab
anks you in advance for you help,
Daniel Hayes
[[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
Some good suggestions, just (as always) be aware of floating-point imprecision.
See FAQ 7.31
> s <- seq(1,30,0.1)
> s[8]
[1] 1.7
> s[8] == 1.7
[1] FALSE
Just trying to forestall future questions :-)
Dan
Daniel Nordlund
Bothell, WA USA
> -Original Message-
> From
e.g., google mail SMTP server
"server = smtplib.SMTP('smtp.gmail.com:587')",
"server.ehlo()",
"server.starttls()",
"server.ehlo()",
"server.login(username,password)",
"server.sendmail(fromaddr, toaddrs, msg.as_string())",
"
also do not know how create.post() handles server authentication (though,
my blind guess would be with the settings of your email program or browser
mail). To stop guessing, if you want a solid comparison, I am afraid you
have to do it yourself.
Best,
Daniel
--
View this message in context:
http
As I (thought I) understood from the sendmailR manual, the package does
currently not support server authentication, or does it?
Daniel
--
View this message in context:
http://r.789695.n4.nabble.com/Email-out-of-R-code-tp3530671p3536512.html
Sent from the R help mailing list archive at
on a group's responses
(i.e., it allows for differences in mean responses between groups) and that
this constant group effect is independent of the effect of femaleset and
treatment on a group.
Best,
Daniel
--
View this message in context:
http://r.789695.n4.nabble.com/lmer-with-2-random
wondering whether there is an alternative
implemented in a different package that does exactly that.
Thanks,
Daniel
--
View this message in context:
http://r.789695.n4.nabble.com/Latent-class-analysis-selection-of-the-number-of-classes-tp3545538p3545538.html
Sent from the R help mailing list archive at
Thanks everybody.
--
View this message in context:
http://r.789695.n4.nabble.com/Latent-class-analysis-selection-of-the-number-of-classes-tp3545538p3554357.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing li
wondering if anyone knew if there were any materials on
the web that would be suitable? and is possible had both questions and
answers.
Many thanks
Dan
PS Cross-posted to http://stats.stackexchange.com
--
**
Daniel Brewer, Ph.D.
Institute
other control variables. You should read up on "omitted variable bias." If
that is not the problem, you have to give us more information/reproducible
code.
Hope that helps,
Daniel
--
View this message in context:
http://r.789695.n4.nabble.com/Linear-multivariate-regression-with-Ro
.
Daniel
--
View this message in context:
http://r.789695.n4.nabble.com/Regressions-with-fixed-effect-in-R-tp2173314p2173326.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman
ivial
problems, the list typically likes to see some own effort and where you are
stuck, rather than to solve the whole problem.
Best,
Daniel
--
View this message in context:
http://r.789695.n4.nabble.com/Advice-needed-on-awkward-tables-tp2173289p2173341.html
Sent from the R help mailing list ar
RE models are available in the lme4 and MASS packages in the glmer and
glmmPQL functions, respectively.
--
View this message in context:
http://r.789695.n4.nabble.com/Panel-data-with-binary-dependent-variable-tp2156043p2184223.html
Sent from the R help mailing list archive at Nabble.com.
__
F test that all u_i=0: F(9, 89) =16.05 Prob > F =
0.
HTH,
Daniel
--
View this message in context:
http://r.789695.n4.nabble.com/Regressions-with-fixed-effect-in-R
parts work. Does the "function()" work when applied on just one
column, for example? If not, then you also have a problem with the
definition of the function. If so, you would have to check, whether the
quantile or IQR function fails in one of the instances.
HTH,
Daniel
--
View this
Xin, you plot the scatterplot wrongly.
Note that the lm (your OLS regression) has a wiggle, whereas you plot
command has a comma. The plot command also should have a wiggle so that you
plot y against x and not x against y. See example below:
x=rnorm(100);e=rnorm(100)
y=2*x+e
reg=lm(y~x)
par(mf
ribed above.
Daniel
--
View this message in context:
http://r.789695.n4.nabble.com/what-the-problem-could-be-if-i-am-suddenly-unable-to-add-abline-to-the-scatter-plot-tp2196304p2196489.html
Sent from the R help mailing list archive at Nabble.com.
__
R-h
.
Daniel
--
View this message in context:
http://r.789695.n4.nabble.com/function-tp2196408p2196576.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
Does anything speak against selecting only x where x$B<20 and then splitting
it by x$C?
split(x[!x$B<20,],x$C)
HTH,
Daniel
--
View this message in context:
http://r.789695.n4.nabble.com/Split-data-frame-by-conditional-and-column-at-the-same-time-tp2197908p2216072.html
Sent from the
Are you asking us to do your homework? This is not a homework list.
For the t-test look in any introductory R manual. For the histograms, look
in the lattice library.
HTH
Daniel
--
View this message in context:
http://r.789695.n4.nabble.com/More-complex-historgram-tp2197823p2216076.html
Sent
Hi, assuming your production figures are in a square matrix, where points
with no production take zero and points with production take the production
figure, you could use
image()
This is the most basic approach I would think.
Daniel
--
View this message in context:
http://r.789695.n4
production,maxColorValue=100))
HTH,
Daniel
--
View this message in context:
http://r.789695.n4.nabble.com/How-to-draw-a-graphic-using-data-with-coordinates-and-production-rate-tp2197568p2217136.html
Sent from the R help mailing list archive at Nabble.com.
___
ed to be discretized this way.
>
> Is there a command in R for this, or some easy shortcut? (I tried
> digging into the lm code, but couldn't find where this is being done.)
>
> Thanks!
>
> -N
>
Noah,
You might try something like
model.matrix(~ group -1)
Hope thi
Xs[[1]]=NULL
#Define a function that regresses y on x and shows us the summary
regress=function(x){summary(lm(y~x))}
#Apply regress over all elements of Xs, i.e.,
#regress y on all possible subsets of regressors
lapply(Xs,regress)
HTH,
Daniel
--
View this message in context:
http://r.7896
After trying again your solution finally it worked with the collate field
in the DESCRIPTION file. The Warnings disappeared.
So indeed I have to thank you very much. Best wishes... Dan
On Sun, 7 Feb 2010, Martin Morgan wrote:
On 02/06/2010 03:39 PM, Daniel Kosztyla wrote:
Hello R-Team,
May
indeed get 4 such warnings
> from example(icuSetCollate), indicating that ICU is not being used.
>
> I have now thrown the above experiment straight at Linux, entering
> command-line commands as follows (with the results shown on the
> lines starting with "#"):
>
> sort
helps,
daniel
--
View this message in context:
http://r.789695.n4.nabble.com/Very-OT-World-Cup-Statistics-tp2246551p2246730.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman
;
LC_ADDRESS="C"
LC_TELEPHONE="C"
LC_MEASUREMENT="C"
LC_IDENTIFICATION="C"
LC_ALL=C
Thanks in advance for your help,
Daniel
TRANSFORME SUAS FOTOS EM EMOTICONS PARA O MESSENGER. CLIQUE AQUI E VEJA COMO.
___
-LatentClassModels.ppt
and an overview is given in this paper by Hoyos:
http://dx.doi.org/10.1016/j.ecolecon.2010.04.011
Anybody has any package to suggest for this kind of analysis? (poLCA does not
provide me the same estimates)
Thanks in advance
Daniel
espondent choice depends by the option characteristics (conditional to the
characteristics of the other options in the same set) + the respondent
socioeconomics.
Daniel
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEAS
eneral solution you also need to consider how you want to deal with
"boundary" conditions. For example, if the x and y coordinates are zero for
both points then the points lie in both the xz and yz planes. And if one of
the coordinates of a point is zero, then how do you dec
?subset
Daniel Nordlund
Bothell, WA USA
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
> On Behalf Of patraopedro
> Sent: Monday, February 07, 2011 6:11 AM
> To: r-help@r-project.org
> Subject: [R] Loop to find d
try xts package and functions:
reclass
period.max
period.min
best,
daniel
2011/2/8 Mike Williamson
> Hello,
>
>This is mostly to developers, but in case I missed something in my
> literature search, I am sending this to the broader audience.
>
>
> - Are there any
gt; Is there a way to avoid that ?
> Thanks in advance,
What kind of processing do you need to do with these numbers? There may be
better ways to handle the problem.
Dan
Daniel Nordlund
Bothell, WA USA
__
R-help@r-project.org mailing list
https://s
Has anyone heard whether Revolution Analytics is going to release this
capability to the R community?
http://www.businesswire.com/news/home/20110201005852/en/Revolution-Analytics-Unlocks-SAS-Data
Dan
Daniel Nordlund
Bothell, WA USA
__
R-help@r
hose values that are not "w"s. So far, i was able to do it but
the result replaces the w's also which i do not want to.
Here is my code that i tried to exclude those w's(within my short codes)
ifelse(mydata[sample,var] != 'w',mydata[sm,var]<
Hi,
I understand that dichotimization of the predicted probabilities after
logistic regression is philosophically questionable, throwing out
information, etc.
But I want to do it anyway. I'd like to include as a measure of fit %
of observations correctly classified because it's measured in units
Hello
Is it possible to get summary statistics (inc mean, sd etc) from a
text file that has the following info stored in it?
Height Frequency
123 5
124 8
125 3
126 9
127 7
etc etc
g for but I just
thought their maybe a better way without using the long variable
Thanks
Dan
On Tue, Feb 22, 2011 at 1:32 PM, Bryan Hanson wrote:
> Daniel, how is the data stored? The answer to your question may be as
> simple as
>
>> df <- read.csv("filename.csv")
>> The first column gives the answer that I am looking for but I just
>> thought their maybe a better way without using the long variable
>
> Sounds like you want the frequency-weighted summary statistics.
> For the mean, look at ?weighted.mean.
> For more comprehensive stats, check ?wtd.mean in th
Dear R Help,
I'm using the sampleSelection library to run a heckman model. The default is
to use a probit model for the selection equation. I would like to know if it
is possible to use a clog-log for the selection equation either with the
sampleSelection library or in a different way. Is this pos
601 - 700 of 1319 matches
Mail list logo