[R] Copula in R 2.4.0
Dear R helper, does anyone have an idea on why R.2.4.0 draws the surface for the two command lines below and the next time it renders the error message below for exactly the same command lines: norm.cop - normalCopula(0.5) persp(norm.cop, dcopula) Error in ceiling(length.out) : Non-numeric argument to mathematical function. I will appreciate any help from anyone thanks, Dominique K. __ R-help@stat.math.ethz.ch 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.
Re: [R] copula
Dear R-helpers, is it some sort of R internal fault or am I doing something wrong when typing the following two lines in R, I get the surface ploted and trying it another time it gives the error message below: Error in ceiling(length.out) : Non-numeric argument to mathematical function Why is it that sometimes I try the same commands and i get the copula density plotted? I will appreciate any help. thanks, Dominique K. __ R-help@stat.math.ethz.ch 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.
Re: [R] copula
Dear R-helper, Is there any thing that I am doing wrong in the following codes: norm.cop - normalCopula(0.5) persp(norm.cop, dcopula) The last command produces what follows Error in persp(x, y, z, xlim, ylim, zlim, theta, phi, r, d, scale, expand, : invalid 'x' argument In addition: Warning messages: 1: no non-missing arguments to min; returning Inf 2: no non-missing arguments to max; returning -Inf 3: no non-missing arguments to min; returning Inf 4: no non-missing arguments to max; returning -Inf 5: is.na() applied to non-(list or vector) in: is.na(x) 6: no non-missing arguments to min; returning Inf 7: no non-missing arguments to max; returning -Inf I remember that I tried the same two commands exactly the way they appear above and I had the surface drawn properly. Does anyone have an idea what this could possibly mean? thanks, Dominique __ R-help@stat.math.ethz.ch 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] r 2.4.0
Can someone help interprete the error message below? i was trying to load the package copula from the R command prompt. Error in loadNamespace(package, c(which.lib.loc, lib.loc), keep.source = keep.source) : in 'copula' methods specified for export, but none defined: show, summary, persp, contour Error: package/namespace load failed for 'copula' Thanks, Dominique __ R-help@stat.math.ethz.ch 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.
Re: [R] perspective
Dear r-users, how can I use the persp or a similar command in R to plot two surfaces in the same set of axes? many thanks, Dominique __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] double integral
Dear r-users, Is there any command in R allowing to evaluate a double integral? for instance let say I want to evaluate the following integral: integrate[lo=(0,1),up=(2,3)] f(x,y)=x^2+y^2 where lo is the vector of lower bounds and up that of upper bounds. I thaught the function adapt would work but it did not. Many thanks, Dominique K. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] help
Dear helpeRs,
I seem to be a little bit confused on the result I am getting from the
few codes below:
u=v=seq(0,1,length=30)
u
[1] 0. 0.03448276 0.06896552 0.10344828 0.13793103 0.17241379
[7] 0.20689655 0.24137931 0.27586207 0.31034483 0.34482759 0.37931034
[13] 0.41379310 0.44827586 0.48275862 0.51724138 0.55172414 0.58620690
[19] 0.62068966 0.65517241 0.68965517 0.72413793 0.75862069 0.79310345
[25] 0.82758621 0.86206897 0.89655172 0.93103448 0.96551724 1.
v
[1] 0. 0.03448276 0.06896552 0.10344828 0.13793103 0.17241379
[7] 0.20689655 0.24137931 0.27586207 0.31034483 0.34482759 0.37931034
[13] 0.41379310 0.44827586 0.48275862 0.51724138 0.55172414 0.58620690
[19] 0.62068966 0.65517241 0.68965517 0.72413793 0.75862069 0.79310345
[25] 0.82758621 0.86206897 0.89655172 0.93103448 0.96551724 1.
f=function(x,y){cp=max(x+y-1,0)}
z=outer(u,v,f)
z is a 30x30 matrix which is fine, but why all its entries are equal to
1? for example, the maximum between u[22]+v[22]-1 and 0 is not 1?? I
don't really know where I went wrong!
thanks,
Dominique
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
