Re: [R] Highly significant intercept and large standard error

2010-10-06 Thread Daniel Nordlund
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

Re: [R] same random numbers in different sessions

2010-10-09 Thread Daniel Nordlund
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@

Re: [R] MATLAB vrs. R

2010-10-11 Thread Daniel Nordlund
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

Re: [R] MATLAB vrs. R

2010-10-11 Thread Daniel Nordlund
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,

Re: [R] Memory limit problem

2010-10-11 Thread Daniel Nordlund
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

Re: [R] How to read only ten rows from a SAS dataset (read.ssd)?

2010-10-19 Thread Daniel Nordlund
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

[R] is get() really what I want here?

2010-10-19 Thread Daniel Weitzenfeld
# 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("

Re: [R] is get() really what I want here?

2010-10-19 Thread Daniel Weitzenfeld
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

[R] cv.lm only bivariate; other options for prediction intervals

2010-10-22 Thread Daniel Weitzenfeld
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

[R] cv.lm() broken; cross validation vs. predict(interval="prediction")

2010-10-22 Thread Daniel Weitzenfeld
<< 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

[R] trouble opening connection 64bit R with win7 64bit

2010-10-27 Thread Hayes, Daniel
> > 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 _

Re: [R] doubt in climate variability analysis in R!

2010-10-29 Thread Daniel Nordlund
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

Re: [R] About 5.1 Arrays

2010-11-05 Thread Daniel Nordlund
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

Re: [R] connecting to remote database using RMySQL

2010-11-05 Thread Daniel Nordlund
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

Re: [R] About 5.1 Arrays

2010-11-06 Thread Daniel Nordlund
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

Re: [R] Using changing names in loop in R

2010-11-06 Thread Daniel Nordlund
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

Re: [R] Help with Iterator

2010-11-09 Thread Daniel Nordlund
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

[R] Clark method enhancement to ChainLadder package

2010-11-13 Thread Daniel Murphy
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

[R] Converting numbers to percentages for boxplots

2010-11-15 Thread daniel carawan
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

[R] Question about GLMER

2010-11-16 Thread Daniel Jeske
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

Re: [R] Question about GLMER

2010-11-16 Thread Daniel Jeske
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

[R] Bug in agrep computing edit distance?

2010-11-16 Thread Dickison, Daniel
; 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

[R] expand comma separated field vertically in data.frame

2010-11-19 Thread Daniel Brewer
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

[R] Changing the Significance Level in R for Hypothesis Tests and Regression Analysis

2010-11-20 Thread Daniel Carawan
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

Re: [R] Counting things in a time series.

2010-11-21 Thread Daniel Malter
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.

[R] weighted Spearman correlation coefficient

2010-11-29 Thread Daniel Rabczenko
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

[R] Suitable test for ordinal variable vs continuous variable trend

2010-12-02 Thread Daniel Brewer
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

Re: [R] evaluating NAs in a dataframe

2010-12-08 Thread Daniel Nordlund
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[

Re: [R] How can i select a set of element in a list ?

2010-12-08 Thread Daniel Nordlund
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]

[R] melt causes errors when characters and values are used

2010-12-10 Thread Daniel Brewer
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

[R] Remove 100 years from a date object

2010-12-10 Thread Daniel Brewer
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

Re: [R] Remove 100 years from a date object

2010-12-10 Thread Daniel Brewer
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="

[R] quantile linear ridge regression using qrnn

2010-12-11 Thread Daniel Gordon
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

[R] SAS "/slice" and "/diff" equivalent? - Automated multiple comparisons in nlme package?

2010-12-12 Thread Daniel Pereira
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 ___

[R] excessive (?) memory utilization by package foreign when reading SAS xport file

2010-12-16 Thread Daniel Nordlund
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

Re: [R] test whether all elements of a vector are identical

2010-12-17 Thread Daniel Murphy
> 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

Re: [R] A question to get all possible combinations

2010-12-22 Thread Daniel Nordlund
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

Re: [R] Reconcile Random Samples

2010-12-23 Thread Daniel Nordlund
^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 __

Re: [R] linear regression with dates

2010-12-27 Thread Daniel Nordlund
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

Re: [R] Reading sas7bdat files into R

2010-12-27 Thread Daniel Nordlund
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

Re: [R] Importing Sas files to R

2010-12-27 Thread Daniel Nordlund
/", "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

Re: [R] Changing column names

2010-12-30 Thread Daniel Nordlund
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

Re: [R] Help with "For" instruction

2011-01-03 Thread Daniel Nordlund
-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

[R] S4 plot function (Package)...

2011-01-04 Thread Daniel Kosztyla
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

[R] aggregate.formula implicitly removes rows containing NA

2011-01-11 Thread Dickison, Daniel
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

Re: [R] aggregate.formula implicitly removes rows containing NA

2011-01-11 Thread Dickison, Daniel
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

Re: [R] R not recognized in command line

2011-01-12 Thread Daniel Nordlund
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 ___

Re: [R] R not recognized in command line

2011-01-12 Thread Daniel Nordlund
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

[R] to append a column to a data frame, has I use loop/if in my case?

2011-01-17 Thread Daniel Wu
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

Re: [R] [R-sig-Debian] Debian ?Ubuntu version of latest R using synaptic inUbuntu 10.10

2011-01-22 Thread Daniel Nordlund
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

Re: [R] (no subject)

2010-07-30 Thread Daniel Nordlund
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,

[R] grep with search terms defined by a variable

2010-08-02 Thread Daniel Malter
. 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

Re: [R] grep with search terms defined by a variable

2010-08-02 Thread Daniel Malter
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

[R] Tables in Packages

2010-08-17 Thread Daniel Murphy
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

Re: [R] how to implement string pattern extraction in R

2010-08-22 Thread Daniel Nordlund
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

Re: [R] lme4 + R 2.11.0 + mac unavailable

2010-08-24 Thread Daniel Malter
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

[R] Time and space considerations in using predict.glm()

2010-08-24 Thread Daniel Yarlett
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

[R] Time and space considerations in using predict.glm.

2010-08-24 Thread Daniel Yarlett
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

[R] Find classes for each column of a data.frame

2010-08-26 Thread Daniel Brewer
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

[R] documentation to upgrade R-package from 32 to 64bit

2010-09-01 Thread Hayes, Daniel
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

Re: [R] create arrays

2011-05-06 Thread Daniel Nordlund
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

[R] Email out of R (code)

2011-05-17 Thread Daniel Malter
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())", "

Re: [R] Email out of R (code)

2011-05-18 Thread Daniel Malter
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

Re: [R] Email out of R (code)

2011-05-19 Thread Daniel Malter
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

Re: [R] lmer with 2 random effects with only two levels

2011-05-19 Thread Daniel Malter
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

[R] Latent class analysis, selection of the number of classes

2011-05-23 Thread Daniel Malter
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

Re: [R] Latent class analysis, selection of the number of classes

2011-05-26 Thread Daniel Malter
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

[R] Any tests/exams freely available with answers to test basic R knowledge

2011-06-09 Thread Daniel Brewer
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

Re: [R] Linear multivariate regression with Robust error

2011-06-10 Thread Daniel Malter
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

Re: [R] Regressions with fixed-effect in R

2010-05-10 Thread Daniel Malter
. 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

Re: [R] Advice needed on awkward tables

2010-05-10 Thread Daniel Malter
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: [R] Panel data with binary dependent variable

2010-05-11 Thread Daniel Malter
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. __

Re: [R] Regressions with fixed-effect in R

2010-05-11 Thread Daniel Malter
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

Re: [R] function

2010-05-12 Thread Daniel Malter
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

Re: [R] what the problem could be if i am suddenly unable to add abline to the scatter plot?

2010-05-12 Thread Daniel Malter
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

Re: [R] what the problem could be if i am suddenly unable to add abline to the scatter plot?

2010-05-12 Thread Daniel Malter
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

Re: [R] function

2010-05-12 Thread Daniel Malter
. 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

Re: [R] Split data frame by conditional and column at the same time

2010-05-13 Thread Daniel Malter
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

Re: [R] More complex historgram

2010-05-13 Thread Daniel Malter
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

Re: [R] How to draw a graphic using data with coordinates and production rate?

2010-05-13 Thread Daniel Malter
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

Re: [R] How to draw a graphic using data with coordinates and production rate?

2010-05-14 Thread Daniel Malter
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. ___

Re: [R] Discretize factors?

2010-05-15 Thread Daniel Nordlund
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

Re: [R] Generating all possible models from full model

2010-05-19 Thread Daniel Malter
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

Re: [R] Posting an 'S4-creating Package Problem'...

2010-05-22 Thread Daniel Kosztyla
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

Re: [R] difference in sort order linux/Windows (R.2.11.0)

2010-05-28 Thread Daniel Nordlund
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

Re: [R] Very OT: World Cup Statistics

2010-06-07 Thread Daniel Malter
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

[R] R in Linux: problem with special characters

2010-06-11 Thread daniel fernandes
; 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. ___

[R] Latent Class Logit Models in discrete choice experiments

2011-01-31 Thread Daniel Vecchiato
-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

Re: [R] Latent Class Logit Models in discrete choice experiments

2011-02-01 Thread Daniel Vecchiato
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

Re: [R] Function to locate points in 3d octants or points on twoaxes

2011-02-02 Thread Daniel Nordlund
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

Re: [R] Loop to find dates whithin dates

2011-02-07 Thread Daniel Nordlund
?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

Re: [R] [Rd] manipulating the Date & Time classes

2011-02-08 Thread Daniel Cegiełka
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

Re: [R] Problem with long number (from character to numeric class)

2011-02-09 Thread Daniel Nordlund
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

[R] Revolution Analytics reading SAS datasets

2011-02-09 Thread Daniel Nordlund
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

[R] creating NAs for some values only

2011-02-13 Thread Daniel M.
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]<

[R] Optimal Y>=q cutoff after logistic regression

2011-02-13 Thread Daniel Weitzenfeld
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

[R] question regarding basic stat question

2011-02-22 Thread Daniel Harris
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

Re: [R] question regarding basic stat question

2011-02-22 Thread Daniel Harris
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")

Re: [R] question regarding basic stat question

2011-02-22 Thread Daniel Harris
>> 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

[R] clog-log based heckit?

2011-02-24 Thread Daniel McDowell
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

<    2   3   4   5   6   7   8   9   10   11   >