calculation for this?
Appreciate your helps in advance.
Cheers! Andre
[[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
give more detail.
Cheers!
On Wed, Oct 1, 2014 at 1:01 AM, Duncan Murdoch
wrote:
> On 30/09/2014 1:31 PM, Andre wrote:
>
>> Dear Sir/Madam,
>>
>> I am trying to use calculation for two-tailed inverse of the student`s
>> t-distribution function presented by Excel
2014, Duncan Murdoch
wrote:
> On 30/09/2014 2:26 PM, Andre wrote:
>
>> Hi Duncan,
>>
>> Actually, I am trying trace the formula for the "Critical value of Z" and
>> manual formula is =(I7-1)/SQRT(I7)*SQRT((TINV(0.
>> 05/I7,I7-2))^2/(I7-2+TINV(0
eers!
On Wed, Oct 1, 2014 at 1:20 AM, Duncan Murdoch
wrote:
> On 30/09/2014 2:11 PM, Andre wrote:
>
>> Hi Duncan,
>>
>> No, that's correct. Actually, I have data set below;
>>
>
> Then it seems Excel is worse than I would have expected. I confirmed R'
incomplete beta function. The incomplete beta function is given by R's
> *pbeta* function. You regularize it with R's *beta* function. Then you
> use R's *uniroot* function to find the inverse. Good homework problem.
>
>
> *Andre >*
> Sent by: r-help-boun.
Dear All,
The manual formula mean that how to calculate that value by hand for
TINV(0.408831, 1221) and the resulted is 4.0891672
Appreciate your help in advance.
Cheers!
On Wed, Oct 1, 2014 at 9:15 PM, wrote:
> What do you mean by a "manual" formula?
>
>
> *Andr
Hello,
I would like to do a dissertation related to the R-Project. This will be the
final work that I have to do to finish my master degree in Information
Management.
Essentially I would like to present a work with the principal
functionalities of R and perhaps I could also do a development of a
Hello,
Anyone knows if there is a way to connect R with the ERP SAP?
Thanks.
--
View this message in context:
http://r.789695.n4.nabble.com/How-to-connect-R-with-SAP-tp2309018p2309018.html
Sent from the R help mailing list archive at Nabble.com.
__
uot; "V" "W" "X" "Y" "Z"
> str( list(letters=letters, LETTERS=LETTERS)[c("letters","LETTERS")] )
List of 2
$ letters: chr [1:26] "a" "b" "c" "d" ...
$ LETTERS: chr [1:26] "A&
Hi Guys,
It's the first time i use R-Help and i hope you can help me.
How can i plot conffidence intervals? with the data bellow:
#Package Austria
library(car)
#head(States)
States1=data.frame(States)
ines=lm(SATM ~ log2(pop) + SATV , data=States1)
summary(ines)
NJ=as.data.frame(States["NJ",c(
Hi Kehl,
First i would like to thank you for the support.
to respond your question i want rather the 4 intervals with the labels, but
another one with the 4 bars with the intervals, it was fantastic!
Thank you again!
2015-05-02 9:24 GMT+01:00 Kehl Dániel :
> Hi Andre,
>
> I think yo
hdquantile step 2 in
hdquantile step 2 in
Tracing function "hdquantile" in package "namespace:Hmisc"
Tracing function "hdquantile" in package "Hmisc"
Also, I can set "RStudio breakpoints" in R Studio on many lines in Misc.s
What is this e
ted levels in factors are deprecated
What am I doing wrong? How do I adjust?
Thanks,
Andre Mikulec
andre_miku...@hotmail.com
DETAILS ( the following R code can be copied and pasted to the R console and
ran )
---
library(nlme)
d.x <- read.table(header=TRUE, text="
o g h i d
A C
, ...) : object 'p' not found
Can someone explain this error and its remedy?
Andre Jackson
jacksonan1...@gmail.com
"You must be the change you wish the world to be"
[[alternative HTML version deleted]]
__
R-help@r-proj
; "G" "H" "G" "H" "G" "H" "G" "H" "G" "H" "G" "H" "G" "H"
$ d: num 1 2 4 8 16 32 64 128 256 1024 2048 4096 8192 16384 32768 65536
- attr(*, "formula")=C
You can't run the Windows GUI in Ubuntu (unless you run it in a Windows
virtual machine). Try installing rkward and using rkward for an alternative
(you will probably be pleasantly surprised):
sudo apt-get install rkward
rkward
http://rkward.sourceforge.net/
-Mose
On Fri, Jul 18, 2008 at 11:09
easy to do it. The help files within GNU R are excellent
quality.
Best regards, André
Andre Cesta, Eng, MSc, SCJP, SCJA
http://geocities.com/aa
l the steps and documented them here:
http://code.google.com/p/jwanalytics/wiki/GnuRe
Thanks for making it easy to do it. The help files within GNU R are excellent
quality.
Best regards, André
Andre Cesta, Eng, MSc, SCJP, SCJA
http://geocities.com/aa
, so that, for example, t1
becomes
0 1 2 3 4
0.4 0.0 0.2 0.0 0.4
Is there a way to do that?
Thanks,
Andre
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R
level after the distributions are created.
Thanks,
Andre
__
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 provide commented, minimal, self-contained, reproducible code.
erely yours
Andre LEBERT
__
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 provide commented, minimal, self-contained, reproducible code.
, and
this would be very useful. I would ideally like to have confidence bands
for the mean function and a single prediction.
I have posted my current Rscript below, with the location of the error
noted, in addition, below the Rscript, is some example data.
Let me know what you think.
Cheers,
|
l +-+ +-+ +-+
label label label
where the margins between each plot should be narrow.
Should I just plot the graphs without axis titles and then use text() to
manually position them?
Thanks in advance,
Andre
__
R-help@r-proj
ylab="")
> plot(1, xaxt="n", yaxt="n", xlab="", ylab="")
> plot(1, xlab="I", ylab="B")
> plot(1, xlab="II", ylab="", yaxt="n")
> plot(1, xlab="III", ylab="", yaxt="n&
should
I use text() and position the title manually?
Thanks in advance,
Andre
__
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 provide commented
text positioning, that is)?
Thanks,
Andre
__
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 provide commented, minimal, self-contained, reproducible code.
r two
plots, something like
++
t ||
||
i ||
++
t
++
l ||
||
e ||
+----+
Thanks,
Andre
__
R-help@r-project.org mailin
plot()? That would make it
easier to position the text.
Thanks again,
Andre
On Mon, 2009-06-01 at 23:02 +0200, baptiste auguie wrote:
> sorry i misread the question -- i thought you wanted the two plots
> were in one row. Here are a few options,
>
>
> 1- use ggplot2 or lattice,
On Mon, 2009-06-01 at 20:46 -0300, Andre Nathan wrote:
> the x and y coordinates seem to always be relative the the axes of the
> first call to plot().
Oops, no they aren't :) It was a mistake in the coordinates. They work
exactly the way I wanted.
sults were the same.
Is there a way to specify a font size that is consistent among calls to
plot()?
Thanks,
Andre
__
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.o
On Thu, 2009-06-11 at 19:09 -0300, Andre Nathan wrote:
> Is there a way to specify a font size that is consistent among calls to
> plot()?
To put it correctly: is there a way to specify a font size that is
consistent among calls to postcript()?
All the fonts in the same file have the sam
in one of them, and par(ps = 18) in the other. Strangely, the generated
files had the following PS code:
/ps 12 def /Font1 findfont 12 s
I don't know PostScript, but I assume this means pointsize = 12.
Best,
Andre
__
R-help@r-project.org mail
Hi All,
I wonder if you can help me with an aparently simple task. I have been
searching examples for this without any luck:
#Assume
x<-1:10 #x ranges from 1 to 10.
y<-x*runif(10)+ 1.5*x #y is a linear function of x with some error. Add
uniform error that is scaled to be larger as x values
Hi All, I wonder if you can help me with an aparently simple task. I have been
searching examples for this without any luck: #Assume
x<-1:10 #x ranges from 1 to 10.
y<-x*runif(10)+ 1.5*x #y is a linear function of x with some error. Add
uniform error that is scaled to be larger as x values a
From: Bert Gunter
To: Andre Cesta
Cc: "r-help@r-project.org"
Sent: Monday, 14 January 2013, 13:34
Subject: Re: [R] How to fit a linear model to data by minimizing the mean
absolute percent error?
Take the logs of both side and minimize the absolute error on the log
scale, then
I started to look at ways to improve times of certain very parallel tasks and
thought that foreach should be a valid candidate to do the job.
So, i opened foreach tutorial by Steve Weston and started timing examples from
it. First example from tutorial is
>system.time(for(i in 1:10) sqrt(i
n 0.4-155
64 bit Win 7 Ultimate
Would greatly appreciate any hints
Andre
[[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.or
quot;
Sent: Tuesday, February 5, 2013 7:22 AM
Subject: Re: [R] failure to connect to Bloomber using Rbbg from batch script on
Windows
Andre,
Does the batch script always fail, or is it sometimes successful? You've
probably tried running in batch and interactively one after the other, bu
I needed to clean someone else code and run it through tidy.source. It
encountered a number of hangups which after some experimentation lead me to
suspect that the culprit is # symbol when it is a part of literal string. For
example if i copy the following to the clipboard
confuse.tidy <-
fu
am using 64-bit R.15.3 version on windows machine
Thanks much
Andre
[[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
Dear all,
I am getting data columnwise that I need to divide by the mean of another
column
If the column is the previous one this code works perfectly well:
fun1 <- function(beginColumn, by, data) { indx <- seq(beginColumn,
ncol(data), by = by) as.data.frame(t(100 - (t(data[, indx])/colMeans(dat
ot seem to find the right option.
If you could help me out, that would be fantastic.
Thanks in advance and have a wonderful day!
--
Sincerely,
Andre Rei Weeks
M.P.H Biostatistics (2014)
B.S. Biology
+1-(850)-443-6592
á§
[[alternative HTML versi
date in my local timezone and avoid this, but if i don't do it, is there an
option in environment or functions to avoid this from happening -- otherwise i
need to watch out for nasty bugs.
Thanks
Andre
[[alternative HTML version deleted]]
__
Hi all,
I got currently some trouble with missvals in netCDF files:
I have a netCDF file, written with an unknown program, which shows via ncdump
no fillval or missval for its variables. Within the variables obviously a
missval of 9.96921e+36 is used.
When I read it into R via the ncdf package
Feb 2014, at 17:45, David W. Pierce wrote:
> > Various questions about missing values in the R ncdf package, and how they
> > are
> > handled if the file lacks the standard "_FillValue" attribute.
>
> Hi Andre,
>
> It sounds like the fundamental pro
(Surv(claimj,censorj==0)~1)
surv.all<-survfit(Surv(claimj,censorj==0)~1)
summary(surv.all)
plot(surv.all)
I would really appreciate any assistance. Thank you.
Regards,
Andre
[[alternative HTML version dele
is code, tmp1 and tmp2 return only one value, where I need
around 520 values. Any ideas would be appreciated.
Regards,
Andre
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing lis
11
[13,] 12
[14,] 13
[15,] 14
[16,] 15
[17,] 16
$series
[1] "x"
$snames
NULL
I need just the $acf values. How do I access them? I would really appreciate
any advice. Thankyou.
Regards,
Andre
[[alt
ere a (quick) way to use reshape without adding a time variable?
2) If the time variable is necessary, is there a quicker way to generate
it? (I know there is, because I did it in Excel...)
Thanks in advance for any advice.
Andre
> foo
Label Value
1 Alpha 0.57911322
2 Alpha 0.02
Hi, everyone. I am using a fair amount of closures in my code. Problem i am
experiencing is i cannot figure out how to mtrace functions defined within a
function. There must be some way to name such function for mtrace to see it
and let me step into it. For example, say i have code
mymodel<-func
After installing RJDBC on RedHat Linux and using it successfully with mySQL
JDBC driver i am trying to use it with Oracle database. I use the JDBC
driver that works fine elsewhere (use it in DBVizualizer). However, when i
try to load the driver,
drv<-JDBC("oracle.jdbc.driver.OracleDriver",
"/hom
I am having trouble loading Oracle JDBC driver into RJDBC.
Basically, i installed RJDBC and used it to connect to MySql. However, when
i try to load JDBC Oracle driver, it gives me an error message
> drv<-JDBC("oracle.jdbc.driver.OracleDriver", "/pathToDriver/ojdbc6.jar")
Error in .jfindClass(as
I asked my local java export and his reply was that the driver was compiled
with the later version of java that i need to install. I did install it in
parallel to existing java and changed JAVA_HOME variable to point to new
java directory. I also pre-pended PATH with path to java. Still, when i tr
Hmm, for the sake of anyone who steps on the same mine. Basically, this class
error always say that you need to install latest java. It's not enough to
get new jre -- you need to get full latest jdk as well. After that you need
to do R reconfiguration, which is done by
R CMD javareconf
Then you ne
ble, one option would be to color each level line
differently and then add a legend. Is it possible to do that?
Finally, how can I remove the tick marks from the top and right axes?
Thanks,
Andre
__
R-help@r-project.org mailing list
https://stat.et
Greetings from Rio de Janeiro, Brazil.
I am looking for advice / references on binary logistic regression
with weighted least squares (using lrm & weights), on the following
context:
1) unbalanced sample (n0=1, n1=700);
2) sampling weights used to rebalance the sample (w0=1, w1=14.29); e
3) a
t; On Apr 27, 2011, at 00:22 , Andre Guimaraes wrote:
>>
>>> Greetings from Rio de Janeiro, Brazil.
>>>
>>> I am looking for advice / references on binary logistic regression
>>> with weighted least squares (using lrm & weights), on the following
o seq() but increasing or
decreasing by factors? I need to do this in one step i.e, not using loops. Any
help would be greatly appreciated.
Regards,
Andre
[[alternative HTML version deleted]]
__
R-h
040th element. I have implemented this with a for loop
an hour ago and it is still loading, can anyone offer any suggestions as to how
I can create this matrix without using loops? I would really appreciate any
suggestions.
Regards,
Andre
[[alterna
Hi,Just a quick one, does anyone know the command for accessing the standard
errors from a survreg object? I can access the coefficients by
model$coefficients, but I cant seem to find a command to access the errors. Any
help would be greatly appreciated.Regards,Andre
hich solve this exact issue.
Cheers,
Andre
**
This email and any attachments are confidential, protect...{{dropped:22}}
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/
s the easiest way to get a list like
(dmin1, dmax1) -> mean(coef such that (dmin, dmax) = (dmin1, dmax1))
(dmin2, dmax2) -> mean(coef such that (dmin, dmax) = (dmin2, dmax2))
...
(dminN, dmaxN) -> mean(coef such that (dmin, dmax) = (dminN, dmaxN))
given a file like
Hi Guys,
Does anyone know how I could constrain my regression coefficients so that they
are positive and add up to one? Any help will be greatly appreciated.
Kind Regards,
Andre
[[alternative HTML version deleted
appreciated.
Kind Regards,
Andre
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
Hi,
Wondering if anyone could help me out with this error.Im trying to fill a
matrix with random numbers taken from an exponential distribution using a loop:
x.3<-matrix(rep(0,3000),nrow=1000,byrow=T)for(i in
1:1000){x[i,]<-rexp(3,rate=2/3)}
I get the error message:
Error in x[i, ] <- rexp(3, ra
uld be appreciated.
Cheers,
Andre
**
This email and any attachments are confidential, protect...{{dropped:22}}
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
he others, or is there an
automatic way to do this?
Thanks in advance,
Andre
__
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 provide commented
rgest value and
call plot() on it, and then points() on the others, or is there an
automatic way to do this?
Thanks in advance,
Andre
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://
's a better way to do that. I've been using
prop.table() to plot these "histograms" for integer numbers, but now I
need binning and I'm not sure what is the best way to proceed.
Thanks in advance,
Andre
__
R-help@r-project.org m
You might try using the Cairo package, which will allow you to bypass using
the Windows graphics device and can give you anti-aliasing. As the manual
says:
"It is ideal for use in server environments (file output) and as a
replacement for other devices that don't have Cairo's capabilities such as
;- cbind(V1,V2,V3);
works. But how could I put together 1200 columns?
I've searched the R mailing help and stumbled upon this entry:
https://stat.ethz.ch/pipermail/r-help/2007-July/137121.html
which doesn't help me.
thanks for your help.
andre
___
uses font size 25 although 28 was requested]
---
system:
R version 2.6.1 (2007-11-26)
on openSuSE 10.3
installed via YaST from
http://download.opensuse.org/repositories/home:/dsteuer/
(provided by Detlef Steuer)
any suggestions?
Thanks!
Hi Joanne
There is a package in Mandriva's repository, called "R-base". You can
isntall it with "urpmi R-base".
HTH,
Andre
On Sat, 2008-09-13 at 16:48 +0100, Joanne Demmler wrote:
> Are there any R rpm's for Mandriva 2008.1? I found a couple of dodgy
&g
Thanks Jim, it worked great!
On Sun, 2008-09-14 at 21:27 -0400, jim holtman wrote:
> try this:
>
> > c(tx,ty)
> 1 2 3 3 4
> 3 2 1 4 1
> > z <- c(tx,ty)
> > tapply(z, names(z), sum)
> 1 2 3 4
> 3 2 5 1
> >
__
R-help@r-project.org mailing list
https://s
ay to "combine" tx and ty in such a way to give me the
distribution below?
1 2 3 4
3 2 5 1
Essentially what I'm looking for is something equivalent to
table(c(x,y)), but x and y are too large and I'd like to avoid the
concatenation.
Thanks in advance,
Andre
This is a text part of the message.
It is shown for the users of old-style e-mail clients
[[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
n searching for
this and couldn't find any info.
Best regards,
Andre
__
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 provide commented, minimal, self-contained, reproducible code.
Thanks Bill and Duncan for your replies, I understand it now. Somehow I
didn't notice the width of the bars was < 1.0.
Andre
On Mon, 2008-02-11 at 11:38 +1000, [EMAIL PROTECTED] wrote:
> Andre,
>
> Regarding your first question, it is by no means clear there is anything
>
ot;
> sprintf("%s.png", name1)
[1] "sample-plot.png"
> paste(name1, "png", sep = ".")
[1] "sample-plot.png"
HTH,
Andre
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
more idiomatic, R-like way to write it:
x <- ...
divs <- ...
for (i in seq(from = 1, to = 100, by = 1000)) {
x[i:(i - 1 + 1000)] <- x[i:(i - 1 + 1000)] / divs[i %/% 1000 + 1]
}
Any suggestions are welcome.
Thanks in advance,
Andre
_
On Mon, 2008-02-25 at 11:52 +1000, [EMAIL PROTECTED] wrote:
> x <- x/rep(divs, each = 1000)
Fantastic :)
Thanks Bill, Andrew and Rolf
Andre
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the p
On Sun, 2008-02-24 at 23:26 -0300, Andre Nathan wrote:
> > gp(1, 2, 10)
> [1]1248 16 32 64 128 256 512 1024
>
Actually,
[1]1248 16 32 64 128 256 512
Andre
__
R-help@r-project.org mailin
On Sun, 2008-02-24 at 21:39 -0500, Charilaos Skiadas wrote:
> 2^(0:9)
I guess it's so simple that I'd never think of that...
Thanks!
Andre
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do rea
8 16 32 64 128 256 512 1024
Any pointers?
Thanks,
Andre
__
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 provide commented, mi
al
data? What I'm doing is setting a small bin length (for example, "breaks
= seq(0, 1, by = 1e-6)", but there's still a chance that points will be
grouped in a single bin.
Is there a better way to do this kind of "raw histogram" plotting?
Thanks,
Andre
_
I know about stem, but the data set has 1 million points, so it's not
very useful here. I want to avoid binning just to have an idea about the
shape of the distribution, before deciding how I'll bin it.
Andre
On Tue, 2008-02-26 at 16:20 -0600, roger koenker wrote:
> t
from a table, so that I can pass them as the x and y arguments to
plot()?
Thanks,
Andre
__
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 p
he prop.table hint. I looked
everywhere but none of my searches hinted at table/table.prop. You guys'
help has been invaluable for me.
Thanks again,
Andre
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
Hey Flo,
I'm not entirely sure what you're problem is -- due to the language thing,
probably -- but I think you want a function that will give you (or computes
in an intermediate step) the "edit distance" between strings. There is the
"sdists" function in the cba package that might help (it can g
Try Sys.sleep()
On Fri, Mar 14, 2008 at 4:33 PM, Rosaura Cordero <[EMAIL PROTECTED]> wrote:
> hola, si me pudieran decir como hago para colocar dentro de un ciclo una
> pausa
>
> -
>
> ¿Con Mascota por primera vez? - Sé un mejor Amigo
>
>
>[[alternative HT
x27; here
Is there a more efficient and/or idiomatic way to do this?
Thanks,
Andre
__
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'm dealing with the following problem:
I have a table with x and y coordinates and corresponding values of
a mineral concentration, let's call it z.
Can someone provide me a short step-by-step manual for the steps
necessary to get a contour map?
How to sort and interpolate my matrix
p/2003-December/043897.html
https://stat.ethz.ch/pipermail/r-sig-geo/2005-October/000609.html
cheers,
andre
S Ellison wrote:
> Doesn't the lattice package's contourplot take arbitrary x and y?
>
> library(lattice)
> ?contourplot
___
Hi,
I would like to use the package netalg made for Stefan Neubauer and Georg
Dorffner. In this package there are functions similars with netlab, a package
for Matlab. I searched for this package and didn't found it. This package had
his name changed??? How I can download this package??
Be
Hi,
I would like estimate a model for function of production's Coob-Douglas using
maximum likelihood. The model is log(Y)= beta[1]+beta[2]*log(L)+beta[3]*log(K).
I tried estimate this model using the tools nlm ( ) and optim ( ) using the
log-likelihood function below:
> mloglik <- fu
95 matches
Mail list logo